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Abstract 

This is the first Internet course on elementary quantum mechan- 
ics written in Spanish ( "castellano" ) for the benefit of Spanish speak- 
ing students. I thank my eight Mexican students at the Institute of 
Physics, University of Guanajuato, Leon, (IFUG), for the collaboration 
that made this possible. The topics included refer to the postulates of 
quantum mechanics, one-dimensional barriers and wells, angular mo- 
mentum, WKB method, harmonic oscillator, hydrogen atom, quantum 
scattering, and partial waves. 
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0. Introduction 



The energy quanta occured in 1900 in works of Max Planck (Nobel prize, 
1918) on the black body electromagnetic radiation. Planck's "quanta of 
light" have been used by Einstein (Nobel prize, 1921) to explain the pho- 
toelectric effect, but the first "quantization" of a quantity having units of 
action (the angular momentum) belongs to Niels Bohr (Nobel Prize, 1922). 
This opened the road to the universality of quanta, since the action is the 
basic functional to describe any type of motion. However, only in the 1920's 
the formalism of quantum mechanics has been developed in a systematic 
manner and the remarkable works of that decade contributed in a decisive 
way to the rising of quantum mechanics at the level of fundamental theory 
of the universe from the mankind standpoint and one of the most successful 
from the point of view of technology. Moreover, it is quite probable that 
many of the cosmological misteries may he disentangled by means of various 
quantization procedures of the gravitational held leading to our progress in 
understanding the origins of the universe. On the other hand, in recent 
years, there is a strong surge of activity in the information aspect of quan- 
tum mechanics, that has not been very much exploited in the past, aiming 
at a very attractive "quantum computer" technology. 

At the philosophical level, the famous paradoxes of quantum mechanics 
(showing the difEculties of the 'quantum' thinking) are actively pursued ever 
since they have been first posed. Perhaps the famous of them is the EPR 
paradox (Einstein, Podolsky, Rosen, 1935) on the existence of elements of 
physical reality, or in EPR words: "If, without in any way disturbing a 
system, we can predict with certainty (i.e., with probabihty equal to unity) 
the value of a physical quantity, then there exists an element of physical 
reality corresponding to this physical quantity." Another famous paradox 
is that of Schrddinger's cat which is related to the fundamental quantum 
property of entanglement and the way we understand and detect it. What 
one should emphasize is that all these delicate points are the sourse of many 
interesting experiments (such as the so-called "teleportation" of quantum 
states) pushing up the technology. 

Here, we present eight elementary topics in nonrelativistic quantum me- 
chanics from a course in Spanish ("castellano") on quantum mechanics that 
I taught in the Institute of Physics, University of Guanajuato (IFUG), Leon, 
Mexico, during January-June semester of 1998. The responsability of the 
idiom belongs mostly to the eight students, respectively. 

Haret C. Rosu 
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0. Introduccion General 

Los cuantos de energia surgicron cn 1900 como consecucncia de los trahajos 
de Max Planck (prcmio Nohcl 1918) sobre el problema de la radiacion del 
cuerpo negro. Los "cuantos de luz" de Planck fueron usados por Einstein 
(premio Nobel 1921) para explicar el efecto fotoelectrico, pero la primer a 
"cuantiRcacion" de una cantidad que tiene las unidades dc una accion (cl 
moniento angular) sc dcbc a Niels Bohr (prcmio Nobel 1922). Eso abrio cl 
camino de la universalidad de los cuantos ya que la accion es la funcion basica 
para describir cualquier tipo de movimiento. Aun asi, solo los aiios veinte se 
consideran como el inicio del formalismo cuantico que levanto a la mecanica 
cuantica al nivel de teoria fundamental del universo y una de las mas exi- 
tosas en cuanto a la tecnologia. En verdad, es muy probable que muchos de 
los misterios cosmologicos estan por ejemplo detras de las diferentes man- 
eras de cuantiRcar el campo gravitacional y tales avances puedcn contribuir 
al entendimiento de los origincs del universo. Por otro lado, el aspecto in- 
formatico de la mecanica cuantica, que no se aprovecho en el pasado, se esta 
desarrollando de una manera muy activa en los ultimos aiios con el proposito 
de investigar la posibilidad dc la construccion de las llamadas "computado- 
ras cuanticas". En la parte Rlosofica cabe mencionar que en la mecanica 
cuantica hay paradojas famosas que todavia se mantienen en polemica y 
que reEejan las dificultades de la logica que impone la manera de pensar 
cuantica (y/o probabilistica) . Una de las mas cclebrcs es la de Einstein (que 
nunca acepto por completo la MC), Podolsky y Rosen (EPR, 1935) sobre 
si hay o no "elementos verdaderos de la realidad fhica" (segun Einstein la 
MC prohibe la existencia independiente del acto de medicion de los sistemas 
fisicos). Otra de igual celebridad cs la del "gato de Schrddinger" . Lo que 
se debe subrayar es que todos cstos puntos teoricos delicados generan ex- 
perimentos muy interesantes (como son por ejemplo los de la Uamada "tele- 
transportacion" de estados cuanticos) que impulsan a la tecnologia. Lo que 
sigue son algunos temas introductivos en la mecanica cuantica. norclativista 
que sirvieron como base para el curso de macstn'a dc mecanica cuantica I en 
el IFUG durante el semestre Enero-Junio de 1998. Este curso fue impartido 
por mi a los estudiantes enlistados, los cuales se encargaron de los temas 
correspondientes. La responsabihdad del idioma pertenece en gran parte a 
cada uno de los estudiantes. 

Haret C. Rosu 
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1. LOS POSTULADOS DE LA MC 



Los siguientes 6 postulados se pueden considerar como la base de la teona 
y los experimentos de la mecanica cuantica. 

PI.- A cualquier cantidad fisica L le corresponde un operador Hermitiano 
L. 

P2.- A cualquier estado (fisico) estacionario de un sistema fisico le corre- 
sponde una funcion de onda normalizada tp (\\ i/j |p= 1). 

P3.- La cantidad fisica L puede tomar solo los valores propios del operador 
L. 

P4.- Lo que se mide es siempre el valor promedio L de la cantidad L en el 
estado ip, la cual en teoria es el elemento de matriz diagonal 

<f\L\f>=L. 

P5.- Los elementos de matriz de los operadores coordenada cartesiana y 
momento Xi yp^, calculados entre las funciones de onda fyg satisfacen 
a las ccuaciones de Hamilton de la mecanica clasica en la forma: 

i<f\m\9>=-<f\§A9>,i<f\^^\9>=<f\§\9> 

donde H es el operador Hamiltoniano. 

P6.- Los operadores Pi y Xk tienen los siguientes conmutadores: 

[Pi,Pk\ = 
[xi,xl] = 0, 

h = h/2'K = 1.0546 X 10-27 erg.seg. 
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1. - La correspondencia de una cantidad fisica L que tiene un analogo 

clasico L{xi,pk) se hace sustituyendo xi, pk por Xj, pf^. La funcion L 
se supone que se puede desarrollar en seric do potencias. Si la funcion 
no contiene productos XkPk, ei operador L es directamente hermitiano. 
Ejemplo: 

T = {Ehi)/'^rn ^ f = (Ef p')/2m. 

Si L contiene productos XiPi, L no es hermitiano, en tal caso L se susti- 
tuye por A la parte hermitica de L (A es un operador autoadjunto). 
Ejemplo: 

w{xi,pi) = T,iPiXi — ^ w = l/2Ei (Pi^i + XiPi). 
Resulta tambien que el tiempo no es un operador sino un parametro. 

2. - (Probabilidad en el espectro discreto y continuo) Si jpn es funcion 

propia del operador L, entonces: 

L =< n \ L \ n >=< n\ Xn\ n >= A„ < n | n >= Snn^n = ^n- 
Tambien se puede demostrar que = (A„)^. 

Si la funcion (j) no es funcion propia de L se usa el desarroUo en un 
sistema completo de L, entonces: 
Sean las siguientes deRniciones: 

combinando estas dos definiciones obtenemos lo siguiente: 

Con las definiciones pasadas ya podremos calcular los elementos de 
matriz del operador L. Entonces: 
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<4>\L\(f)>= J2n,m Kn'^nK < m \ n >= J2m I «m P >^m, 

lo cual nos dice que el resultado del experimento es Xm con la proba- 

bilidad \ a,n p. 

Si el espectro es discrete: de acuerdo con el postulado 4 eso signiBca 
que I Qm o sea, los coeGcientes del desarroUo en un sistema completo 
dctcrminan las probabilidades de observar el valor propio A„. 
Si el espectro es continuo: usando la siguiente definicion 

<f>{T) = Ja{X)i;{T,X)dX, 
se calcularan los elementos de matriz para el espectro continuo 

L =< (f) \ L \ (f) > 
= JdTja* (A)V'* (r, X)dX J Ma(/x) ^'(t, n)dii 
= J J a*a{ij,)ii J iI}*{t, X)'tjj(tau, fj,)dXdiJ,dT 
— II a* {X)a{fj,)fj,d{X — iJ.)dXdfi 
= Ja*{X)a{X)XdX 
= J I a(A) p XdX. 

En el caso continuo sc dice que \ a (A) p cs la densidad de probabilidad 
de observar el valor de X del espectro continuo. Tambien vale 

L =< (f) \ L \ (f) >. 
3.- Definicion de la derivada con respecto a un operador: 

dFik - F{L+ei)-F(L) 
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4.- (Representacion del momento) Cual es la forma concreta de pi, pi J 
P3, si los argumentos de las funciones de onda son coordenada carte- 
siana xi. 

Vamos a considerar el siguiente conmutador: 

iPij -^i \ ~ Pi-^i Pi 
— Pi^i^i ^iPi^i ~t~ ^iPi^i ^i^iPi 
~ {,Pi^i ^iPi)^i ~f~ ^iipi^i •^iPi} 

— \pi jXi\Xi -\- Xi \pi , Xi\ 

= —ihxi — ihxi = —2ihxi. 

En general se tiene: 

Pi^i Pi — TliTlX'i 

Entonces para todas las funciones anaUticas se tiene lo siguiente: 

Piil^{x) - '4^{x)pi = 

Ahora sea pi(p = f{xi,X2,xs) la accion de pi sobre X2, xa) = 1. 
Entonces: 

Pii/j = —i^-§^ + flip y hay relaciones analogas para x^ y xz- 
Del conmutador \pi,pk] = se obtiene V x / = 0, por lo tanto, 
h = V,F. _ 
La forma mas general de Pi es: 

Pi = —ih-^ + donde F cs cualquicr funcion. La funcion F se 
puede eliminar utilizando una transformacion unitaria W = expfi^. 

Pi = uH-in£ + §:)u 
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= -^^^ 

resultando que pi = —ifi-£^ — ^ P = —ihV. 

5.- (Calculo de la constante de normalizacion) Cualquier funcion de onda 
i>{x) G de variable x se puede escribir como: 

Hx) = J Six - c)md^ 

y considerar la expresion como desarroUo de ip en las funciones propias 
del operador coordenada x6{x — ^) = Ci^ ~ 0- Entonces, | V'(a^) P es 
la densidad de probabilidad de la coordenada en el estado il^{x). De 
aqm resulta la interpretacion de la norma 

II ipix) |p= / I ipix) p dx = 1. 

El sistema descrito por la funcion ip{x) debe encontrarse en algiln lugar 
del eje real. 

Las funciones propias del operador momento son: 
—ih^ = Pii^, integrandola se obtiene '4>{xi) = A expt*''^*, x yp tienen 
espectro continuo y entonces se tiene que hacer la normalizacion con 
la funcion delta. 

Como se obtiene la constante de normalizacion? 

se puede obtener utilizando las siguientes transformadas de Fourier: 

f{k) = J g{x) e^p~^''^ dx, g{x) = f{k)exp^'^^dk. 

Tambien se obtiene de la siguiente manera: 

Sea la funcion de onda no normalizada de la particula libre 

ipx 

(j)p{x) = ^exp n y la formula 

S{x-x') = ^ exp^*^(^-^') dx 

se ve que 
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= A* exp~ft~ Aexp~fr dx 



fOO I A |9 ix(p-p ) 

/f^ I ^ |2 exp n dx 



. ,0 . „0O ) 



= 27r;i I A |2 S{p-p) 
entonces la constante de normalizacion es: 

V2TTn 

Tamhien resulta que las funciones propias del operador momento for- 
man un sistema completo (en el sentido del caso continuo) para las 
funciones L^. 



'^^^'^ ^ vhh ^ ^^^^ ^ 

^'^^^ ^ ^ ^^^^ exp-??^ dx. 
Estas formulas establecen la conexion entre las representaciones x y p. 

6.- Representacion p: La forma explicita de los opcradores pi y x^ se 
puede obtener de las relaciones de conmutacion, pero tamhien usando 
los niicleos 



x{p-il3) = U^xU = 2^ J exp n xexp n dx 



-, -iyx Q i0x 

2^jexp ft {-ifiQ^exp ft ). 
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La integral pasada tiene la forma siguiente: 

M(A,A') = JU^X,x)MU{X' ,x)dx, y usando xf = J x{x,^)f{^)d^. 
Entonces la accion de x sobre a(p) G es: 

xa{p) = J x{p, P)a{l3)dl3 
= ! "i^V^ {-iT^^(iW^)dx)a{l3)dl3 

— ipx fx i0x 

= I exp~s- ^ exp— a{(3)dxd(3 
= exp— ^ exp- a{P)d^d(3 

= ^//exp^ ^a{P)d^dp 
= -inj^6iP-p)df3 = -ih^, 
donde S{P - p) = ^ J exp^^^^ d|. 

El operador momento en la representacion p se caracteriza por el 
nucleo: 

pip, 13) = U^pU 
= 2^1 J exp {-ihg^)exp h dx 

-j —ipx i0x 

= 2il J f3exp h dx = I35(jp - (3) 

resultando que pa{p) = pa{p). 

Lo que pasa con x y p es que aunque son hermiticos sobre todas f(x) 
G no son hermiticos sobre las funciones propias. 
Si pa{p) = Poa{p) y X = x'^ p = p^ , entonces: 
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<a\px\a> — <a\xp\a >= —ih < a \ a > 

Po[< a\x\a> — <a\x\a>] = —ih < a \ a > 

Po[< a\x\a> — <a\x\a>]=0 

La parte izquierda es cero, mientras tanto la derecha esta indefinida, 
lo que es un contradiccion. 

7- (Representaciones de Schrddinger y Heiscnberg) 

Las ecuaciones de movimiento dadas por el postulado 5 tienen varias 
interpretaciones, por el hecho de que en la expresion ^ < / | L | / > 
uno puede considerar la dependencia del tiempo atribuida completa- 
mente a las funciones de onda o completamente a los operadores. 

• Para un operador dependiente del tiempo O = 0{t) tenemos: 

ri. - -dS ^.-dS 
^» ~ dxi ' ~ dpi 

\P, f]=Pf-fP = -if^m: 

[x, f]= xf - fx = -ifi§i: 

se obtienen las ecuaciones de movimiento de Heisenberg: 



Pi = T'\p,H], Xi = ^[x,H]. 

• Si las funciones dependen del tiempo, todavia se puede usar 
Pi = ^\pi,H], porque es consecuencia solo de las relaciones de 
conmutacion y entonces no dependen de la representacion. 



i<f\Pi\9>=-i<f\ \p,H]\g>. 

Si ahora pi y H no dependen del tiempo y teniendo en cuenta su 
hermiticidad se obtiene: 
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= :^{fJ,Hg) + ^{f,Hp,g) 
= ^{pf,Hg) + l(Hf,p,g) 

{% + iHf,p,g) + {p,f,^-iHg)=0 

La ultima rclacion se cumple para cualquier pareja de funciones 
f{x) y g{x) al momento inicial si cada una satisface la ecuacion 

ih^ = Hi,. 

Esta cs la ecuacion de Schrddinger y la descripcion del sistenia 
por operadores independientes del tiempo se llama representacion 
de Schrddinger. 

En las dos representaciones la evolucion temporal del sistema se carac- 
teriza por el operador H, el cual se obtiene de la funcion de Hamilton 
de la mecaxiica clasica. 

Ejemplo: H de una particula en potencial U{xi,X2,X2) es: 
^ ~ 2m ^ (a^i) 3^2, x^), y en la representacion x es: 

H = -^V + U{X1,X2,X3). 

8.- El postulado 5 vale en las representaciones de Schrddinger y de Heisen- 
berg. Por eso, el valor promedio de cualquier observable coincide en 

las dos representaciones, y cntonccs, hay una transformada unitaria 
que pasa de una representacion a otra. Tal transformacion es del tipo 

s' = exp R . Para pasar a la representacion de Schrddinger hay 
que usar la transformada de Heisenberg tp = f con f y L, y para 
pasar a la representacion de Heisenberg se usara la transformacion 
de Schrddinger A = s^Ls con ip y A. Ahora so obtendra la ecuacion 
de Schrddinger: como en la transformacion tp = s^f f la funcion f no 
depende del tiempo, derivaremos la transformacion con respecto al 
tiempo obteniendose lo sig.: 
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por lo tanto, tenemos: 

Enseguida calcularemos la ecuacion de Heisenberg: poniendo la trans- 
formacion de Schrddinger de la siguiente manera sAst = L y derivan- 
dola con respecto al tiempo se obtiene la ecuacion de Heisenberg 

^ = IfAst + sA^ = i#exp^ Ast - fsAexp^ H 

= K-ffsAst - sAs^H) = ^{HL - LH) = L]. 
Por lo tantOjtenemos: 

Tambien la ecuacion de Heisenberg se puede escribir de la sig. manera: 

§ = ^s[H,A]s^. 

A L se le conoce como la integral de movimiento si < ip \ L \ tp >= 
y esta caracterizada por los siguentes conmutadores: 

[H,L] = 0, [H,A] = 0. 

9.- Los estados de un sistema descrito por las funciones propias do H 
se Uanian estados estacionarios del sistema, y al conjunto de valores 
propios correspondientes se les Uaman espectro de energia (espectro 
energetico) del sistema. En tal caso la ecuacion de Schrddinger es : 
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Y su solucion es: tpn{x,t) = exp 's"* (t)n{x). 

• La probabilidad es la siguiente: 

S{x) =1 il^n{x,t) p=| exp^f^ p 
= exp R exp R I (pn[x) 1^ = 1 0„(x) 1^. 
Resultando que la probabilidad es constante en el tiempo. 

• En los estados estacionarios, el valor promedio de cualquier con- 
mutador de tipo [H, A] es cero, donde A es cualquier operador: 

< n I HA - Ah I n >=< n \ HA | n > - < n | AH \ n > 
=< n I En-A \ n> — <n \ AEn \ n > 
= En <n \ A\n> -En <n \ A\n>=0. 

• Teorema del virial en mecanica cuantica.- Si H es el operador 
Hamiltoniano de una particula en un campo U (r) y usando 
A= l/2J2i=i{PiXi — XiPi) se obtiene lo siguiente: 

< V I [i, -ff] I V- >= =< ^ I i-ff - I V > 

= ELi < I PiXiH - HpiXi I V > 

= Ef=i <-il^\[H, Xi]pi + Xi[H,pi] I V >• 

usando varias veces los conmutadores y pi = —ihA, H = T+U{r), 
se tiene entonces: 

<ip\[A,H]\tP>=0 
= -in{2 <^|f|V>-<V'l^- VC/(r) I i/j >). 



16 



Que es el teorema del virial. Si el potencial es U{r) = Uor''\ en- 
tonces tenemos el teorema del virial como en mecanica clasica, 
solo que para valores promedios 

T=^U. 

• Para un Hamiltoniano H = + J7(r) y [f, H] = -^p, y 

calculando los elementos de matriz se tiene: 

{Ek - En) < n \ r \ k >= < n \ p \ k >. 
- (Densidad de corriente de probabilidad) La siguiente integral : 

I I V'n(a;) P dx = 1, 

es la normalizacion de una funcion propia de un espectro discreto 
en la representacion de coordenada, y ocurre como una condicion de 
movimiento en una region Bnita. Por eso, los estados del espectro dis- 
creto se Uaman estados ligados. 

Para las funciones propias del espectro continuo ipx{x) no se puede dar 
de manera directa una interprctacion de prohahilidad. 
Supongamos una funcion dada ^ G L^, ia cual la escribimos como 
combinacion lineal de funciones propias en el continuo: 

4> = I aW'ip\{x)dx. 

Se dice que (j) corresponde a un movimiento infinito. En muchos casos, 
la funcion a(A) es diferente de cero solo en una vecindad de un punto 
A = Ao- -En este caso (f) se le conoce como paquete de onda(s). 
Vamos a calcular el cambio en el tiempo de la probabilidad de encon- 
trar el sistema en el volumen fi. 

P = In \^{x,t) dx = f^ip*{x,t)ilj{x,t)dx. 
Derivando la integral con respecto al tiempo tenemos: 
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Utilizando la ecuacion de Schrddinger del lado derecho de la integral 
se tiene lo siguiente: 



Usando la identidad fVg — gVf = div[{f)grad{g) — {g)grad{f )\ y la 
ecuacion de Schrddinger de la forma siguiente: 



Sustituyendo lo anterior en la integral se obtiene: 

Usando el teorema de la divergencia para transformar la integral de 
volumen en una de superGcie, entonces tenemos lo siguiente: 



f = -/^(W*-V'*vv)dx. 

A la cantidad J{ip) = ^{ipVip* - se le conoce como densidad 

de corriente de probabilidad y de inmediato se obtiene una ecuacion 
de continuidad, 



§ + div{J) = 0. 

• Si iIj{x) = AR(x), donde R{x) es funcion real, entonces: J{i^) = 
0. 
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• Para las funciones propias del impulso ^^{x) = (27r?t)^/2 ^■^P 
ohtiene: 

t/t / 2i'{) \ _ p 

lo cual nos dice que no depende de la coordenada la densidad de 
probabilidad. 

- (Opcrador de transporte espacial) Si H es invariante ante translaciones 
de cualquier vector a, 

H{r + a) = H{r). 

Entonces, hay un T{a) unitario T\d)H{r)T[a) = H{r + a). 
Por la conmutacion de las translaciones 

f (a)f (6) = f (5)f (a) = f (a + 6), 

resulta que T tienc la forma T = exp*'^'^, donde, k = ^. 
En el caso infinitesimal: 

f{6a)Hf{6a) ^ {I + ik6a)H{I - ikSa), 

H{r) + i[k, H\5a = H{r) + {VH)6d. 

Tambien [p, H] = 0, donde p es una integral de movimiento. El sis- 
tema tiene funciones de onda de la forma ip{p,r) = (27rnV/2 ^^P"^ y 

ipa 

transformada unitaria hace que exp n '^(f) = '0(r + a). 
El operador de transporte espacial = exp~n~ es analogo al 
= exp h el operador de transporte temporal. 
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12.- Ejemplo: Si H es invariante con rcspccto a una traslacion discreta (por 
ejemplo en una red cristalina) H{f+ a) = H{f) donde a = J2i (H'lT'i, 
e N y tti son los vectores baricos. 
Entonces: 

H{f)i){f) = E^{r), 

H{r-\- a)V'(f + a) = Eip{f+ a) = H{r)ip{f+ a). 

Resultando que i/j{r) y i/j{r-\- a) son funciones de onda para el mismo 
valor propio de H. La rclacion cntrc %lj{r) y ip{f+ a) sc pucdc buscar 
en la forma ilj{f + a) = c{a)'ip[r) donde c{a) es una matriz gxg (g es 
el grado de degeneracion del nivel E). Dos matrices tipo c, c{a) y c{b) 
conmutan y entonces son diagonalizables en el mismo tiempo. 
Adcmas para los clcmcntos diagonales hay dc tipo 
Cii{a)cii{b) = Cii{a+h), donde i=l,2,....,,g. Lasolucidn de esta ecuacion 
es del tipo Cii{a) = exp*'^'" resulta que = Uk{f) exp*'^" donde k es 

un vector real cualquiera, y la funcion Ukif) es periodica con periodo 
a, entonces: Uk{r + a) = Ukif)- 

La aseveracion de que las funciones propias de un H periodico cristal- 
ino H{r + a) = H{r) se pueden escribir ipkir) = Uk{r)expika con 
Uk{f-\- a) = Uk{r) se llama teorema de Bloch. 

En el caso continuo, Uk debe ser constante, porque la constante es la 
unica funcion periodica para cualquicr a. 

El vector p = Uk se llama cuasi-impulso (analogia con el caso con- 
tinuo). El vector k no esta determinado de manera umvoca, porque se 
le puede juntar cualquier vector g para el cual ga = 2-?™ donde n e N. 
El vector g se puede escribir g = bimi donde son numeros 

enteros y hi estan dados por 

* di{d'jXd'k) 

si i j k son los vectores baricos de la red cristalina. 

Citas: 

1. Acetatos del Prof. H. Rosu. 
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2. E. Farhi, J. Goldstone, S. Gutmann, "How probahility arises in quantum 
mechanics", 

Annals of Physics 192, 368-382 (1989) 

3. N.K. Tyagi cn Am. J. Phys. 31, 624 (1963) da una dcmostracion muy 
corta del principio de incertidumbre de Heisenberg, que dice que la medicion 
simultanea de das operadores hermitianos que no conmutan produce una 
incertidumbre relacionada con el valor del conmutador. 

Notas: 

1. Para "la creacion de la MC...", Werner Heisenberg ha sido galardon- 
ado con el premio Nobel en 1932 (y lo recibio en 1933). El articulo "Zur 
Quantenmechanik. U", Zf. f. Physik 35, 557-615 (1926) (recibido el 16 de 

Noviembre de 1925) de M. Born, W. Heisenberg y P. Jordan, sc Ic conocc 
como "el trabajo de las tres gentes" y esta considerado como el que abrio 
de verdad los grandes horizontes de la MC. 

2. Para "la interpretacion estadistica de la funcion de onda" Max Born 
recibio el premio Nobel en 1954. 
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Problemas 



Problema 1.1; 

Considerar dos operadores A y B los cuales por hipotesis, conmutan. 
Entonces se derivara la relacion: 

exp^ exp^ = exp"^+^ exp^/^["^'^I, (formula de Glauber). 

DeGniremos un operador F(t), una funcion de variable real t, por: 
F{t) = exp^*exp^*. 

tenemos: 

'-^ = Aexp^* exp^* + exp^* B exp^* = {A + exp^* B exp-^*)F(i). 

Ahora aplicando la formula [A,F{B)] = [A, B]F' (B), tenemos que 
[cxp^*, B] = t[A.B] exp^*, por lo tanto: exp^* B = B exp^* +t[A, B] exp^* 
multiplicando ambos lados de la ecuacion pasada por exp^"^* y sustituyendo 
en la primera ecuacion, obtenemos: 

^ = {A + B + t[A,B])F{t). 

Los operadores A y B y [A, B] conmutan por hipotesis. Por lo tanto, 
podemos integrar la ecuacion diferencial como si A + B y [A, B] fueran 
numeros. 

Entonces tenemos: 

F{t) = F(0)exp(^+-B)*+V2[A,B]t2_ 

Poniendo t = 0, se ve que F{0) = 1, y : 
F{t) = exp(^+^)*+V2[A,B]t2 _ 

Entonces poniendo t = 1, obtenemos el resultado deseado. 
Problema 1.2; 
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Se calculara el conmutador [X, Dx] ■ Para hacerlo, tomaremos una funcion 
arbitraria jp{f): 

[X,Dx]i^{r) = {xl - i^x^r) = x^^(r) - l[x^{f)] 
= x^ip{f) - ip{f) - x-^^ipir) = -V'(f). 
Entonces si es valido para toda ijj{r), se puede deducir que: 

[X,Dx] = -l. 

Problema 1.3; 

Se probara que la traza es invariante ante un cambio de base ortonormal 
discreta. 

La suma de los elementos de la diagonal dc la matriz la cual representa 
un operador A cn una base arbitraria no depende dc la base. 
Se derivara esta propiedad para el caso del cambio de una base ortonormal 
dicreta [\ Ui >] a otra base ortonormal dicreta [| tk >]■ Tenemos: 
T,i<Ui\A\ui>=J2i<Ui\ [Efc I tk >< tk WA I Ui > 

(donde se ha usado la relacion dc cerradura para el estado tk). El lado 
derecho de la relacion pasada es igual a: 

<Ui\tk><tk\ A\ui >= Y.i,j <tk\ A\ui><Ui\tk >, 

(es posible el cambio de orden del producto de dos numeros). Entonces, 
podemos reemplazar J2i | ^i^ >< «i | por uno (relacion de cerradura para el 
estado \ Ui >), y se obtiene Gnalmente: 

J2i < Ui\ A \ Ui >= J2k < ^fc I ^ I ^fc >• -for io tanto, se ha mostrado la 
propiedad de invariancia para la traza. 
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2. POTENCIALES BARRERAS y POZOS 

Comportamiento de una Funcion de Onda Esta- 
cionaria iIj{x) 

Regiones de Energia Potencial Constante 

En el caso de un potencial cuadrado, V{x) es una funcion constante 
V{x) = Ven cierta region del espacio. En tal region, la ecuacion de Schrddinger 
puede ser escrita: 



^^(x) + ^{E-V)i,ix)=0 (1) 



Distinguiremos entre varies casos: 

(i) E>V 

Introduzcamos la constante positiva k, definida por 

fc = i^?™fHS (2) 

La solucion de la ecuacion (1 ) puede ser entonces escrita: 

iPix) = Ae'^"" + A'e-'^'' (3) 

dondc A y A' son constantes complejas. 

(ii) E <V 

Esta condicion corresponde a regiones del espacio las cuales estarian 
prohibidas para la particula por las leyes de la mecanica clasica. En este 
caso, introducimos la constante positiva q definida por: 

,= «HS (4) 

y la solucion de (1 ) puede ser escrita: 

ij{x) = Se«^ + 5'e-«^ (5) 

donde B y B' son constantes complejas. 

(iii) E = V 

En este caso especial, ip{x) es una funcion lineal de x. 
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Comportamiento de V'(x) en una discontinuidad de energia 
potencial. 

Podna pensarse que en el punto x = xi, donde el potencial V{x) es 
discontinuo, la funcion de onda il){x) se comportaria extranamente,llegando 
a ser por si misma discontinua, por ejemplo. Este no es el caso: ip{x) y ^ 
son continuas, y es solo la segunda derivada la que es discontinua en x = xi 

Vision general del calculo 

El procedimiento para determinar el estado estacionario en un "poten- 
cial cuadrado" es por lo tanto el siguiente: en todas las regiones donde 
V{x) es constante, escribimos ip{x) en cualquiera de las dos formas (3) o (5) 
scgun sea. aplicable; entonces pegamos estas funciones por requerimientos de 
continuidad de ^{x) y de en los puntos donde V{x) es discontinuo. 

Examinacion de ciertos casos simples 

Llcvemos a cabo cl calculo cuantitativo de los estados estacionarios, 
hecho de acuerdo al metodo descrito arriba. 
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Potencial escalon 



V(x) 



V 



II 



Fig. 2.1 



a. Caso donde E > Vq; reflexion parcial 
Pongamos la ec. (2) como: 



h = 



V2mE 

n 



(6) 



ko = 



V2m(E - Vq) 

n 



(7) 



La solucion de la ec. (1) tiene la forma de la ec. (3) en las regiones 
I{x < 0) yllix > 0); 



ikix I /)/ —ikix 



En la region I la ec. (1 ) toma la forma: 



iIj"{x) + ^^^(x) = ^"{x) + k'^^ix) = 

h 
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y en la region 11: 

2m 



^"(x) -^[Vo- E]ct>{x) = r{x) - q^x) = 



Si nos limitamos al caso de una particula incidente viniendo desde x = —oo, 
debemos elegir A2 = y podemos determinar los radios A[/Ai y A2/A1. 
Las condiciones de pegado entonces dan: 

• ipj = ipjj, en X = : 

Al + A[ = A2 (8) 

• ■(/'/ = "0//; en X = : 

Aiiki - A[iki = A2ik2 (9) 

Sustituyendo Ai y A[ de (8) en (9): 

A, = Mi±M („, 

La igualacion de la constante A2 en (10) y (11) resulta: 

A[ _ki-k2 



Al ki + k2 



(12) 



y un despeje en (11) nos da: 



A2 _ 2ki 

Ti = hTh 



^(x) es la supcrposicion de dos ondas. La primcra (el termino cn Ai) 
corresponde a una particula incidente, con momento p = hki, propagandose 
de izquierda a derecha. La segunda (el termino en A[) corresponde a una 
particula reBejada, con momento —hki, propagandose en la direccion op- 
uesta. Ya que hcmos clcgido A2 = 0, ipn{x) consistc de una sola onda, la 
cual esta asociada con una particula transmitida. (En la siguiente pagina 
se muestra como es posible, usando el conccpto de una corriente de proba- 
bilidad, definir el coeGciente de transmision T y el coeficiente de reflexion 
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R del potencial escalon). Estos coeficientes dan la probabilidad de que la 
particula, Uegando de x = — oo, pase el potencial escalon en x = o se 
regrese. Asi encontramos: 

i?=|4^P (14) 



y, para T: 



T='A^l^ (15) 



Tomando en cuenta a (12) y (13), tenemos: 

(fci + 

Es facil vcriGcar que R+T = 1; cs cicrto que la particula sera transmitida 
o reRejada. Contrariamente a las predicciones de la mecanica clasica, la 
particula incidente tiene una probabilidad no nula de regresarse. 

Finalmente es facil vcriRcar, usando (6) y (7) y (17), que, si E ^ Vo, T c^i 
1: cuando la encrgia de la particula es suGcientemente grande comparada 
con la altura del potencial escalon, la particula salva este escalon como si 
no existiera. 

Considerando la solucion en la region I: 



J = -^(0*V<^-^V^*) (18) 
con Aie^^^^ y su conjugado A\e~'^'^^^: 

j = -^[(Ate-^*^i^)(^izA;ie*^)-(y4ie^'=i^)(-^tzfcie-*^)] 

J = ^\A,f 
m 

Ahora con A'le^'-^'^'-^ y su conjugado A\e''^^''^ resulta: 

m 
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Deseamos en seguida veriGcar la proporcion de corriente que se reBeja 
con respecto a la corriente que incide (mas precisamente, queremos veriGcar 
la probabiUdad de que la particula se regrese): 

j,_\ji4>-)\ \ - _ , 

En forma similar, la proporcion de lo que se transmite con respecto a 
lo que incide (o sea la probabiUdad de que la particula se transmita) es, 
tomando ahora en cuenta la solucion de la region II: 



iMllyl 121 kiAi 

a. Caso donde E <Vq; reflexion total 
En este caso tenemos: 



=^'2 (20) 



^/2mE 

k = — (21) 



= «^ (22) 

En la region I{x < 0), la solucion de la ec. (1) [dada como iIj{x)" + kftp{x) = 
Oj tiene la forma de la ec. (3): 

iPi = Aie''''^ + A'le-^'^i^ (23) 

Y, en la region II {x > 0), la misma ec. (1) [ahora dada como il){x)" — 
q2tp{x) = 0] tiene la forma de la ec. (5): 

^pn = S2e«2^ + S^e-«2=^ (24) 

Para que la solucion permanezca limitada cuando x +oo, es necesario 
que: 

B2 = (25) 
Las condiciones de pegado en x = dan en este caso: 

» ipj = ipjj^ en X = : 

Ai + A[ = B'^ (26) 
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• = ^jj, en X = 



Aiiki - A[iki = (27) 



Sustituyendo Ai y A'^ de (26) en (27): 



La igualacion de la constante B2 en (28) y (29) resulta: 

A'l _ iki +q2 _ ki — iq2 
Ai iki — q2 ki+ iq2 ' 

y un despeje en (29) nos da: 

B'o 2iki 2ki 



(30) 



(31) 



Ai iki — 92 ^1 — iq2 
El coeBciente de reEexion R es entonces: 

i?=l4^P = I^^P = 4^ = l (32) 
Ai 'ki+iq2 kf + q^ 

Como en la mecanica clasica, la parti'cula es siempre reBejada (reBexion 
total). Sin embargo, hay una diferencia importante: debido a la existencia 
de una onda evanescente e"^^^, la particula tiene una probabilidad no nula 
de estar presente en la region del espacio la cual, clasicamente, le ser/a 
prohibida. Esta probabilidad decrece exponencialniente con x y Uega a 
ser despreciable cuando x es mas grande que el "rango" l/q2 de la onda 
evanescente. Notemos tambien que el coeBciente A[/Ai es complejo. Un 
cierto cambio de fase aparece a causa dc la rcBexion, el cual, fisicamonte, es 
debido al hecho de que la particula cs rctardada cuando penetra la region 
X > 0. No hay analogia en la mecanica clasica. 
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Potencial Barrera 



V(x) 
Vo ^_ 



II 



III 



Fig. 2.2 

a. Caso donde E > Vo; resonancias 
Pongamos aqu/ia ec. (2) como: 



h = ^ (33) 



= vM|ZS (34) 

La solucion de la ec. (1) tiene la forma de la ec. (3) en las regiones 
I{x < 0), 7/(0 <x<a) y III{x > a) : 

■4)1 = Aie'^^i^ + yl'ie-^'^i^ 

Si nos limitamos al caso de una particula incidente viniendo desde x = 
— oo, debemos elegir A'^ = 0. 



31 



• ipl = ipii, en X = : 

Ai + A[ = A2 + (35) 



• V/ = i^'iij en X = : 

Aiiki - A[iki = A2ik2 - A2ik2 (36) 

• ipii = '4'ilij en X = a : 

Aae^'^^" + A'2e-'''^'' = A^e^^^" (37) 

• = ipjii, en X = a : 

A2ik2e'''^'' - A'2ik2e-'''^'' = Azikie'^^"" (38) 

Las condiciones de continuidad en x = a dan entonces a A2 y A2 cn funcion 
de A3, y aquellas en x = dan a Ai y A'^ en funcion de A2 y A'2 (y, 
consecuentemente, en funcion de A^). Este proceso es mostrado enseguida: 
Sustituyendo A2 de la ec. (37) en (38): 

A3e^'^"(fc2 + fci) 

^2 = ^j—^a (39) 

Sustituyendo A2 de la ec. (37) en (38): 

^ ^ A,e^^n^2 - fci) 



Sustituyendo Ai de la ec. (35) en (36): 



, _ A2{k2 - ki) - A'2{k2 + ki 



(41) 



Sustituyendo A'^ de la ec. (35) en (36): 



A2{k2 + ki) - A'2{k2 - ki 
= 2k[ ^^^^ 



Ahora, sustituyendo en (41) las ecuaciones (39) y (40): 



A[ = iM_^(sinfc2a)e*M3 (43) 
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Y, finalmente, sustituyendo en (42) las ecuaciones (39) y (40): 



Ai = [cos k2a - i \ , ^ sin k2a]e^^'^'^ (44) 
2k\k2 

A'l/Ai y As/Ai [relaciones que salen de igualar las ecuaciones (43) y (44), 
y del despeje de la ec. (44), respectivamente] nos capacita para calcular el 
coeBciente de reEeccion R y el de transniision T de la barrera, los cuales 
aqm son iguales a: 

Entonces es facil de veriRcar que R + T = 1. 
b. Caso donde E <Vq; efecto tunel 

Ahora tendnamos a las ecuaciones (2) y (4) dispuestas: 

h = ^ (47) 



La solucion de la ec. (1) tiene la forma de la ec. (3) en las regiones 
I{x < 0) 7 III{x > a) y la forma de la ec. (5) en la region 11(0 < x < a): 

Las condiciones de pegado enx = Oyx = a nos capacita para calcular el 
coeficiente de transmision de la barrera. De hecho, no es necesario realizar 
otra vez el calculo: todo lo que debemos hacer es sustituir, en la ecuacion 
obtenida en el primer caso de esta misma seccion, k2 por —iq2. 

Estados Ligados; Pozo de Potencial 
a. Pozo de profundidad finita 
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V(x) 





Vn 









Fig. 2.3: Pozo finite 

En esta parte nos limitaremos solo a tratar el caso < E < Vq (el caso 
E > Vq es exactamente igual al calculado en la seccion precedente, "barrera 
de potencial" . 

Para las regiones exteriores I {x < 0) y III {x > a) usamos la ec. (4): 



^ ^2m{Vo - E) 
Para la region central II {0 < x < a) usamos la ec. (2): 



(49) 



k = 



n 



(50) 



La solucion de la ec. (1) tiene la forma de la ec. (5) en las regiones 
exteriores y la forma de la ec. (3) en la region central: 
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En la region {0 < x < a) la ec. (1 ) toma la forma: 
2mE 

y en las regiones exteriores: 

2m 

'n 

Ya que debe ser limitada en la region I, debemos tener: 



^{x)" + -^'^{x) = iP{x)" + k'^iPix) = (51) 



iPix)" - ^ [Vb - E]4>{x) = ilj{x)" - q^i,{x) = (52) 



Si = (53) 



Las condiciones de pcgado dan: 
■0/ = fpiij en X = : 



Bi = A2 + A'^ (54) 

ip'j = ip'jj, en X = : 

Biq = A-2ik - A'2ik (55) 

^11 = 'tplll, enx = a: 

A^e''"' + A'^e-'^" = Sge'^" + S^e'^" (56) 

i^ii = i^'iiD enx = a: 

A2ike''"' - A'^ike-'^" = B^qe'^" - B^ge'^" (57) 

Sustituyendo la constante A2 y la constante A2 de la ec. (54) en la ec. 
(55) obtenemos, respectivamente: 

^ ^ Bi{q-ik) 
^ -2ik 

_ Bi{q + ik) 

- 2ik ^^^^ 

Sustituyendo la constante A2 y la constante A2 de la ec. (56) en la ec. 
(57) obtenemos, respectivamente: 

B'^e-'^''(ik + q)+B^e'i''(ik-q) + A'2e-'^''{-2ik) = 
2ikA2e^^'' + B'^e-'i''(-ik + q)+B^E'i''{-ik-q) = (59) 
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Igualando B'^ de las ecuaciones (59) y tomando en cuenta las ecuaciones 
(58): 

ft p~i'^ 

— = ^[e'^^iq + ikf - e-^^°(g - ikf] (60) 
B\ Aikq 

Pero ip{x) debc tambicn estar limitada en la region III. Pot lo tanto, es 
necesario que = 0, esto es: 



q + ik e" 



-ika 



Ya que q y k dependen de E, la cc. (1 ) pucdc scr satisfecha para ciertos 
valorcs dc E. Imponiendo un limite sobrc il){x) en todas las regiones del 
espacio as/ se ocasiona la cuantizacion de la energia. Mas precisamente dos 
casos son posibles: 

(i) si: 



q + ik 

Si en esta expresion igualamos en ambos miembros la parte real y la imagi 
naria, nos resulta: 

,ka q 



tan(-) = J (63) 



Poniendo: 



l2mVo 



ko = ^—^ = ^k^ + q^ (64) 
Entonces obtenemos: 

1 1 + tan\^) = = (^)2 (65) 



cos2(^) ' 2 ' P ^ k 

La ec.(63) es as/ equivalente al sistema de ecuaciones: 



tan(y) > (66) 

Los niveles de energia estan determinados por la interseccion de una Imea 
recta teniendo una inclinacion con arcos senoidales (Uneas largas inter- 
rumpidas en la Ggura 2.4). As/ obtenemos un cierto numero de niveles de 
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encrgia, cuyas funciones de onda son pares. Esto es mas claro si sustituimos 
(62) en (58) y (60); Es Mcil veriGcar que B'^ = Bi y que A2 = A2, as/ que 
tp{—x) = ijj{x). 
(ii) si: 



q — ik 
q + ik 

Un calculo del mismo tipo nos Ueva a: 



(67) 



I sin( 



ka. 



k_ 
ko 



,ka. 
tan(— ) < 



(68) 



Los nivclcs dc encrgia cstan entonccs dctcrminados por la intcrscccion 
de la misma Hnea recta como antes con otros arcos senoidales (Imeas cortas 
en la Ggura 2.4). Los niveles asi obtenidos caen entre aquellos encontrados 
en (1). Puede ser facilmente mostrado que las correspondientes funciones de 
onda son impares. 




n/a 



27C/a 



37C/a ko 47t/a 



Fig. 2.4 



a. Pozo de profundidad infinita 
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Asumiendo que V{x) es cero para < x < a e infinito en cualquier parte 
mas. Pongamos: 

4>{x) debe ser cero fuera del intervalo [0, a], y continuo en x = 0, tambien 

como en x = a. 
Ahora para < x < a: 

iIj{x) = Af^^"" + A'e"^*^^ (70) 
Ya que ^(0) = 0, puede ser deducido que A' = —A, lo cual nos Ueva a: 

ip{x) = 2iAsm{kx) (71) 

Ademas 'tp{a) = 0, as/ esto: 

k = ^ (72) 

donde n es un arbitrario entero positivo. Si normalizamos la funcion (71 ), 
tomando (72) en cuenta, entonces obtenemos las funciones de onda esta- 
cionarias: 

f2 . ,mTX, 



V'n(x) = ^-sin(— -) (73) 
con energias: 

La cuantizacion de los niveles de energia es asi, en este case, particularmente 
simple. 
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Problemas 



Problema 2.1: El potencial Delta 

Supongamos que tenemos el potencial de la forma: 

V{x) = -Vodix); Vb > 0; x e ^. 

La correspondiente funcion de onda ip{x) se supone que es suave. 

a) Busca los estados ligados (E <0) los cuales estan localizados en este 
potencial. 

b ) Calcula la dispersion de una onda plana Uegando en este potencial y 
encuentra el coeficiente de reflexion 



2 



R 



\Aef\ 



a;=0 



donde ipref, i^lleg GS la onda refiejada y la de Uegada, respectivamente. 
Sugerencia: Para evaluar el comportamiento de iIj{x) en x=0, Integra la 
ecuacion de Schrddinger en el intervalo (—e, +e) y considera el Umite e 
0. 

Solucion. a) La ecuacion de Schrddinger esta dada por 



Lejos del origen tenemos una ecuacion diferencial de la forma 

2m.E 

Las funciones de onda son por lo tanto de la forma 



^(x) + ^{E + Vod{x))i;{x) = (75) 



V'W = -^V'W- (76) 



tp{x) = Ae~*^ + 5e«^ si x>0 a x<0, (77) 



con q = \J —2mE /U'^ G 3ft. Como debe ser integrable, alii no puede 
haber una parte incrementandose exponencialmente. Ademas la funcion de 
onda debe ser continua en el origen. De aqm, 

^(x) = Ae«^; (x<0), 

ij{x) = Ae-«^; (x>0). (78) 
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Integrando la ecuacion de Schrddinger desde —e a +£, obtenemos 

- ^bP'i^) - - VbV'(O) = E i;{x)dx ^ 2eEiP{0) (79) 

Insertando ahora el resultado (78) y tomando el Umite £ — oo, tenemos 

-^{-qA-qA)-VoA = (80) 

o E = — m(V^^/2^^). Claramente hay solo un eigenvalor de energia. La 
constante de normalizacion es encontrada que es A = y^mVo/^^. 



vx-x) = ^e*'=^ e = ^ (81) 



b) La funcion de onda de una onda plana es descrita por 

2mE 

Se mueve de izquierda a derecha y es reBejada en el potencial. Si B o C es 
la amplitud de la onda reflejada o transmitida, respectivamente, tenemos 

ipix) = Ae'^"" + Be-'^''; (x < 0), 

^{x) = Ce'^''; (x>0). (82) 

Las condiciones de continuidad y la relacion ip'^e) — ijj'{—e) = —fip{0) 
con f = 2mVo/n'^ da 

f . 



A+B = C B 



f + 2ik^ 



2jk 

ik{C-A + B) = -fC c = j^A. (83) 

El coeGciente de reflexion deseado es por lo tanto 

Si el potencial es extremadamente fuerte (Vq ^ oo) R ^ 1, o sea, la onda 
entera es reflejada. 

El coeficiente de transmision es, por otro lado, 

Si el potencial es muy fuerte, (Vq oo) T ^ Q, o sea, la onda transmitida 
se desvanece. 

Obviamente, R + T = 1 como se esperaba. 
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Problema 2.2: Particula en un Pozo de Potencial finito Uni- 
dimensional 

Resuelve la ecuacion unidimensional de Schrodinger para un pozo de poten- 
cial Rnito descrito por el siguiente potencial 



Vix) 



—Vq si \x\ < a 
si \x\ > a 



Consider a solo est ados ligados (E < 0). 



-a 



+a 



E 



-Vo 



Fig. 2.5 

Solucion. a) La funcion de ondapara \x\ < a y \x\ > a. La correspondiente 
ecuacion de Schrodinger esta dada por 

+ V{x)Mx) = EMx) (86) 

2m 



De&nimos 
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y obtcncmos: 

1) si X < -a : ^P'|{x) - g^Vi = 0, V'l = + Sie^^"'; 

2) si — a < X < a : (^) + = 0, ■02 = ^2 cos(A;x) + B2 sm(kx); 

3) si X > a : iIj'^{x) - q^ipz = 0, ^3 = -Bge^^ + Sge"''^. 

bj Formulacion de las condiciones de frontcra. La normalizacion de los 
cstados ligados requiere que la solucion sc haga cero en el infinito. Esto 
significa que Bi = = 0. Ademas, ip{x) debe ser continuamente difer- 
enciable. Todas las soluciones particulares son fijadas de tal forma que ip 
tambien como su primera dcrivada ij/ son suaves en ese valor dc x corre- 
spondiendo a la frontera entre cl interior y el exterior. La scgunda derivada 
ip" contiene el salto requerido por el particular potencial tipo caja de esta 
ecuacion de Schrddinger. Todo esto junto nos Ueva a 

V'i(-a) = V2(-a), ^2(0) = 0^3(0), 

i^'ii-a) = V'2(-«), V'2(«) = V'3(«)- (88) 

c) Las ecuaciones de eigenvalores. De (88) obtenemos cuatro ecuaciones 
lineales y homogeneas para los coeRcientes Ai, A2, B2 y B3: 

Aie"'^"- = A2Cos{ka) - B2sm{ka), 

qAie~'^"' = A2ksm.{ka) + B2kcos{ka), 

B^e~'^°' = A2Cos{ka) + B2S\n{ka), 

-qBse-'i"- = -A2ksm{ka)+B2kcos{ka). (89) 

Por adicion y sustraccion obtenemos un sistema de ecuaciones mas facil de 
resolver: 

(^1 + 53)6"'?" = 2A2Cos{ka) 
q{Ai + Bs)e-i'' = 2A2ksm{ka) 

{Ai - ^3)6^''" = -2B2 sin{ka) 
q{Ai - B3)e-'^'' = 2B2kcos{ka). (90) 

Asumiendo que Ai + B3 ^ y A2 0, Las primeras dos ecuaciones resultan 

q = ktaxi{ka). (91) 

Insertando esta en una de las ultimas dos ecuaciones da 

Ai = S3; B2 = 0. (92) 
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De aqm, como resultado, tenemos unasolucion simetrica con ^(x) = ^{—x). 
Entonces hablamos de una paridad positiva. 

Casi identico calculo nos Ueva para Ai — B3 ^ y para B2 7^ a 



-A;cot(A;a) y ^1 = -S3; ^2 = 0. 



(93) 



La funcion de onda as/ obtenida es antisimetrica, correspondiendo a una 
paridad negativa. 

d) Solucion cualitativa del problema de eigenvalores. La ecuacion que 
conecta q y k, la cual ya hemos obtenido, son condiciones para el eigenvalor 
de energia. Usando la forma aorta 



^ = ka, r] = qa, 
obtenemos de la deGnicion (87) 

2 2 2mVoa^ 



(94) 
(95) 



For otro lado, usando (91) y (93) obtenemos las ccuaciones 

77 = ^tan(0, 77 = -Ccot(0- 

For lo tanto los valores de energia deseados pueden ser obtenidos con- 
struyendo la interscccion de esas dos curvas con el circulo definido por (95), 
en el piano ^-77 (ver figura 2.6). 

11 



T) = ^ tan ^ 




2 2 


2 


^ +11 


=r 






r\=-t, cot ^ 


















■ 



Fig. 2.6 



Al mcnos una solucion existe para valores arbitrarios del paramctro Vq, 
en el caso de paridad positiva, porque la funcion tangente intersecta el ori- 
gen. Para paridad negativa, el radio del circulo necesita ser mas grande que 
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un valor mmimo de tal forma que las das curvas puedan intersectarse. El 
potencial deber tener una cierta profundidad en conexion con un tamaiio 
dado a y una masa dada m, para permitir una solucion con paridad nega- 
tiva. El nilmcro dc nivclcs dc cncrgia sc incrcmcnta con Vq, a y la masa m. 
Para el caso en que mVa? — > oo, las intersecciones son encontradas en 



,. , , 2n-l 

tan(A;a) = oo correspondiendo a ka = — - — tt, 

— cot(fca) = oo correspondiendo a ka = mr, 

donde n = 1, 2, 3, ... 
o combinado: 

k{2a) = mr. 
Para el espectro de energia esto signiBca que 



(96) 



(97) 



(98) 



Ampliando el pozo de potencial y/o la masa de la particula m, la diferencia 
entre dos eigenvalores de energia vecinos decrecera. El estado mas bajo 
(n = I) no csta localizado en —Vq, sino un poco mas arriba. Esta diferencia 
es Uamada la energia de punto cero. 

e) La forma de la funcion de onda es mostrada en la figura 2.7 para la 
solucion discutida. 





Fig. 2,7: Formas de las funciones de onda 
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Problema 2.3: Particula en un Pozo de Potencial infinite Uni- 
dimensional 

Resuelve la ecuacion unidimensional de Schrodinger para una particula que 
SB encuentra en un pozo de potencial infinito descrito como sigue: 

, ^ _ J si x' < X < x' + 2a 

1 oo si x' > X o X > x' + 2a 

Tenemos que la solucion en forma general es 

'4j{x) = A sin{kx) + B cos{kx) (99) 

donde 



k = \l^ (100) 



l2mE 

Como ip debe cumplir que il^{x') = ip{x' + 2a) = 0, se tiene: 



A sm{kx) + B cos{kx) = (101) 
Asm[k{x' + 2a)] + B cos[k{x' + 2a)] = (102) 

Multiplicando (101) par s'm[k{x' + 2a)] y (102) por sin(kx') y, enseguida, 
restando el segundo resultado del primero obtenemos: 

B[ cos(fcx') sm[k{x' + 2a)] - cos[k{x' + 2a)] sm{kx') ] = (103) 

o a traves de una identidad trigonometrica: 

Bsm{2ak) = (104) 

Multiplicando (101) por cos[k{x' + 2a)] y restandole la multiplicacion de 
(102) por cos{kx'), se tiene: 

A[ sm{kx') cos[k{x' + 2a)] - sin[k{x' + 2a)] cos(kx') ] = (105) 

o a traves de la misma identidad trigonometrica: 

Asin[k(-2ak)] = Asin[k(2ak)] = (106) 

Como se descarta la solucion trivial ip = entonces por las ecuaciones 
(104) y (106) se tiene que sin(2aA;) = y esto solo ocurre si 2ak = mr, siendo 
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n un entero. Segun lo anterior k = nn/2a y como k"^ = 2mE/'h? entonces 
tenemos que los eigenvalores estan dados por la expresion: 

E = 107 

La energia esta cuantizada ya que solo se le permiten ciertos valores; para 
cada kn = 711^ I2a le corresponde la energia = [n^/2m][7r/i/2a]^. 
La solucion en la forma general queda como: 

V'„ = asin(^) + 5cos(^). (108) 

Normalizan do: 

fx' -\-1a 

1 = / '^ip*dx = a{A^ + (109) 

Jx' 

lo que nos Ueva a 

A = ±^l/a- (110) 
Sustituyendo este valor de A en (101) se obtiene: 

„ 1 . ,mrx\ 

Sustituyendo este valor de B en (110) se tiene: 

1 ,mrx\ 

Tomando los signos superiores de Ay B, y sustituyendo sus valores en (108) 
se tiene: 

1 7?7r 

Utilizando los signos inferiores de A y B se tiene: 

1 Tiir 
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3. MOMENTO ANGULAR EN LA MC 



Introduccion 

Dc la Mecdnica Cldsica se sabe que, el momento angular 1 de las particulas 
macroscopicas esta dado por 

l = rxp, (1) 

donde r y p son respectivamente el vector de posicion y el momento lineal. 

Sin embargo, en Mecdnica Cudntica, el operador momento angular (OMA), 
en general no es un operador que se exprese solamente por el operador co- 
ordenada xj y el operador momento lineal pk, los cuales, solo actuan sobre 
las funciones propias (FP) de coordenadas. Por lo tanto, es muy impor- 
tante establecer, antes que nada, las rclacioncs de conmutacion para las 
componentes del OMA, es decir, en Mecdnica Cudntica 1 se representa por 
el operador 

1 = -ihr X V, (2) 

cuyas componentes son operadores que satisfacen las siguientes reglas de 
conmutacion 

ademas, cada uno de ellos conmuta con el cuadrado del OMA, esto es 

1^ = 11 + 11 + 11 [kA=^, z = 1,2,3. (4) 

Estas relaciones, ademas de ser validas para el OMA, tambien se cumplen 
para el operador espm (OS), el cual, carece de analogo en la mecdnica cldsica. 

Estas relaciones de conmutacion son basicas para obtener el espectro de 
dichos operadores, asi como sus representaciones diferenciales. 

Momento Angular Orbital 

Para un punto cualquiera de un cspacio espacio Bjo (EP), se puede tener 
una funcion il}{x,y,z), para el cual, consideramos dos sistemas caxtesianos 
S 7 S', donde S' se obtiene de rotar el eje z. 

En general un EP se refiere a un sistema de coordenadas diferentes a S y S'. 
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Ahora bien, coniparemos los valores de ip en dos puntos del EP con las 
misnias coordenadas (x,y,z) enY, yT,' , esto es equivalente a rotar un vector 



/ 






( COS 4> 


— s'mcj) 







?) 


- 


sin (j) 


cos (j) 





V 






{ 





z 



^ix',y',z') = R'ip{x,y,z) (5) 
donde R es la matriz de rotacion en R'^ 

s'mcf) \ j X \ 

\ y \ , (6) 

\ z' J \ z ) \ z j 

entonces 

Ri/j{x, y, z) = ip{x cos (j) — y sincp, x sin + y cos (p, z). (7) 

For otro lado es importante recordar que las funciones de onda (FO) no 
dependen del sistema de coordenadas y que la transformacion a rotaciones 
de las FF se hace por medio de operadores unitarios, luego entonces, para 
establecer la forma del operador unitario (OU) U^{(j)) que Ueva tp a ip' , 
siempre se considera una rotacion infinitesimal d(p, manteniendo solamente 
los terminos lineales en d<p cuando se expande ip' en series de Taylor alrededor 
del punto x 

i^{x',y',z') ^ i:{x + yd(p,xd(p + y,z), 

~ i'i^^y^z) +d(p(^y^ - x^ 

{l-id(pl^)'ip(x,y,z), (8) 

donde hemos empleado la notacior^ 

h = fr'^{xpy-y'p.j,), (9) 

la cual, como se vera mas adelante, corresponde al operador de proyeccion 
en z del momento angular de acuerdo con la definicion en (2) y dividido por fi, 
tal que, para la rotacion del angulo (p finito, se pueda representar como una 
exponencial, es decir 

^{x',y',z) = e^'^^tP{x,y,z), (10) 

donde 

IJ\(P) = e^'-<^. (11) 



^La demostracion de (8) se presenta como el problema 3.1 
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Para reaBrmar el conccpto dc rotacion, considcrcmosla dc una manera mas 

general con la ayuda del vector A que actua sobre if) asumiendo que A^, Ay 
Az tienen la misma forma en E 7 S', es decir, que los valores promedio de 

A calculados en S 7 S' deben ser iguales cuando se ven desde el EF, esto es 
J nf^){Aj + Ayj' + A,k')^*{f^) dr 

= Jr (0 (A J + Ayj + A,k)r (r) df, (12) 

donde 

i! = icoscf) + jsiiLcf), y = zsin^ + jcos ^, k' = k. (13) 
Luego entonces, si combinamos (10), (12) y (13) tenemos 

e^'^'^i^e"*'^'^ = cos (/) - ij, siiK/), 
e'^^'f'Aye-'^^'^ = A^csincf)- Ay cos (p, 
e^'-'^i^e"*^-* = A^. (14) 

Nucvamcntc consideremos rotaciones inEnitcsimalcs y desarroUando las 
partes de la izquicrda en (14) se pueden determinar las relaciones de con- 
mutacion de A,j., Ay y A^ con 1^ 

[lz,Ax] = iAy, [lz,Ay] = —iAx, [lz,Az] = 0, (15) 

7 de manera similar para Ix y ly 

Las condiciones basicas para obtener tales relaciones de conmutacion son 

•k La FP se transforma como en (7) cuando S — >^ S' 

•k Las componentes Ax, Ay, Az tienen la misma forma en T, y T,' 

-k Los vectores de los valores promedio de A en E 7 S' coinciden (son 
iguales) para un observador en el EF. 

Tambien se puede usar otra representacion en la cual la FO ip{x,y,z) 
no cambia cuando S — > E' 7 7o operadores vectoriales se transforman como 
vectores. Para pasar a tal representacion cuando (f) rota alrededor de z se 
usa el operador tJ{^), es decir 
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e'^^W{x,y,z) = ij{x,y,z), (16) 

entonces 

^-iiz4> Xf^^'^'t' = A. (17) 
y utilizando las relaciones dadas en (14) 

= Axcoscp + Ay sin 4> = e-'^^^Axe'^^'f', 
A'y = -Axsm4> + Ay cos (l) = e-'^^'^ Ay e'^^"^, 

= (T'^^'^A^e^^^'''. (18) 

Dado que las trasformacioncs dc la nueva representacion se hacen con 
operadores unitarios, las relaciones de conmutacion no se cambian. 

Observaciones 

★ El operador A? es invariante ante rotaciones, esto es 

★ Resulta que 

[k,A^]=0, (20) 

★ Si el operador Hamiltoniano es de la forma 

6 = ±.f + u(f). (21) 

entonces se mantiene invariante ante rotaciones dadas en cualquier eje 
que pasa por el origen de coordenadas 

[k,H]=f) (22) 

donde las li son integrates de movimiento. 

Definicion 

Si Ai son las componentes de un operador vectorial que actua sobre una 
FO de coordenadas y si hay operadores li que cumplen con las siguientes 
relaciones de conmutaciones. 

[li, Aj] = ieiji;Ai~, [li,lj] = i^ijk^k- (23) 
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Luego cntonces, las li se Uaman componentes del operador momento an- 
gular orbital, y podemos concluir de (20) y (23) que 

[ii,P] = 0. (24) 

Por lo tanto, las trcs componentes asociadas con las componentes de un 
momento angular clasico arbitrario satisfacen las relaciones de conmutacion 
Mas aun, puede mostrarse que el origen de estas relaciones conducen a 
propiedades geometricas de rotacion en un espacio tridimensional. Esto 
es porque adoptamos un punto de vista mas general y dcfiniremos un mo- 
mento angular J como cualquier conjunto de tres observables Jx, Jy y Jz los 
cuales cumplen las relaciones de conmutacion, es decir 

[Ji, Jj] = iSijkJk- (25) 
Entonces introducimos el operador 

i^ = Jl + Jl + Jl (26) 

el (escalar) cuadrado del momento angular J. Este operador es hermitiano, 
dado que Jx, Jy y Jz son hermitianos, y asumiremos que es un observable, 
ahora mostremos que conmuta con las tres componentes de J 

[J2,J]=0. (27) 

Por el hecho dc que 3^ conmuta con cada una de sus componentes, 
entonces hay un sistema completo de FP 

J^V'tm = tV'tm, Ji'^-yn = JJ-i'-yiJ.- (28) 

Los operadores Ji y Jk {i 7^ k) no conmutan, es decir, no tienen FP en 
comiin. 

En lugar de usar las componentes Jx y Jy del momento angular J, es 
mas conveniente introducir las siguientes combinaciones lineales 

J-\- ~ Jx ~t~ "^Jyi J~ ~ Jx iJy' (2*^) 

Los cuales no son hermiticos, tal como lo operadores a y del oscilador 

armonico, solamcntc son adjuntos. 

Ahora, si hacemos las operaciones indicadas en la izquierda concluimos 

que 

[Jz, J±] = ±J±, [J+, J-] = 2Jz, (30) 
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[J\J+] = [J\J.] = [J\J,]=^. (31) 

Jz{J±'>P-ff,) = {J±Jz + [Jz, J±]}'^-1^l = (/i ± 1)^±V'7M- (32) 

Luego entonces las ip^^ son FP de los valores de (autovalores) /x ± 1, 
es decir 

J-ipji^ = Pi^tpj^-i. (33) 

pero 

taJ que, aJ tomar una fase del tipo e*" (con a real) para la funcion ^/^^^ se 
puede hacer real e igual a (3^, esto quiere decir 

J+lp-yix-l = Q;mV'7M' = aMV'7M-i' (^5) 

y por lo tanto 

7 = '(^J'y^,[Jx+Jy+Jz]i^'y|X = ^J''^+a + b, 

^ = {'<P'yn,j'^1p'Yix) = {Jyll^'YH,Jy'<P'yn)>0. (36) 

La norma de cualquier funcion es no negativa, esto implica que 

7 > m', (37) 

para 7 fijo, el valor do fi tiene Hmite superior e inferior (es decir, tiene valores 
en un intervalo finito). 

Sean A y A esos limites (superior e inferior de n) para un 7 dado 

J+V7A = 0, J-V7A = 0. (38) 
Ahora, utilizando las siguientes igualdades operatoriales 

J-J+ = + J, = J2_ J^(J^_l), 

J+J_ = + = J2_ J^(J^ + l)^ (39) 

actuando sobre ip^A J i^-yX se obtiene 

7-A2-A = 0, 
^-A^ + A = 0, 
(A-A + 1)(A + A) = 0. (40) 
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La condicion 

A> A^A = -A = J^7 = J(J+1). (41) 

Para un 7 dado (fijo) la proyeccion del momento ^ toma 2 J + 1 valores 
que difieren por una unidad de J a — J. For eso, la diferencia A — X = 2 J 
debe ser un numero entero y consecuentemente los autovalores Jz tornados 
por m son enteros, esto es 

m = k, k = 0,±l,±2, ... , (42) 

o semienteros 

m = k + ^, k = 0,±l,±2, ... . (43) 

Para un estado dado 7 = J {J +1), es degenerado con grado g = 2J +1 
respecto a los autovalores de la posicion del momenta m (esto es porque 
Ji, Jk conmutan con pero no conmutan entre ellos). 

Por "estado de momento angular J" se entiende en la mayon'a de los 
casos, un estado con 7 = J( J+1) en el cual, el valor maximo de la proyeccion 
es J. Tales estados se dcnotan por o \jm). 

Vamos a encontrar los elcmcntos de matriz de Jx, Jy de la representacion 
en la cual, y Jz son diagonales, luego entonces de (35) y (39) se tiene que 

J—J+'4^jm—l — OlmJ—^jm — C>lm'4^jm—li 

J{J + 1) - (m - if - (m - 1) = a^, 

am = ^(J + m)(J-TO + l). (44) 
Ahora, combinando (44) con (35) se obtiene 

J+i^jm-i = \J{J + m){J -m + l)'ipjm, (45) 
resulta que el elemento de matriz de J+ es 

{jm\J+\jm - 1) = ^J{J + m){J -m + l)5nm, (46) 
y de manera analoga 

{jn\J-\jm) = {J + m){J - m + l)Snm-i, (47) 
por ultimo, de la deHnicion en (29) J"^, J_ se obtiene facilmente 

{jm\Jx\jm-l) = ^^(J + m)(J-m+ 1), 
{jm\Jy\jm-l) = ^^{J + m){J-m + l), (48) 
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A Modo de Conclusion 

a (Propiedades de los Eigenvalores de J y J^) 

j{j + y jz son los eigenvalores de 3 y Jz asociados con los 
eigenvectores \kjm), entonces j ym satisfacen la desigualdad 

~j < m < j. 

(3 (Propiedades del Vector J-\kjm)) 

Sea \kjm) un eigenvector de y jm con los eigenvalores j{j + 1)^^ 7 
mh 

- (i) Si m = - j, J-\kj - j) = 0. 

— (a) Si m > —j, J-\kjm) es un eigenvector no nulo de y Jz 
con los eigenvalores j{j + l)h^ y [m — l)h. 

7 (Propiedades del Vector J+ 1 kjm) ) 

Sea \kjm) un eigenvector de 3^ y Jz con los eigenvalores j{j + l)h y 

■k Si m = j, J^\kjm) = 0. 

-k Si m < j, J+ = \kjm) es un vector no nulo de 3^ y jz con los 
eigenvalores j{j + 1)^^ y (m + l)h 

6 Luego entonces 

Jz\kjm) = mh\kjm).. 

J+\kjm) = rnh^ j{j + 1) — m{m + l)\kjm + 1), 

J-\kjm) = mhyj j{j + 1) — m(m — l)\kjm + 1). 

Aplicacion del Momento Angular Orbital 

Hemos considerado las propiedades del momento angular derivadas unicamente 
de las relaciones de conmutacion, ahora retomaremos el momento angular L 
de una particula sin giro y veremos como dicha teoria desarroUada en la 
seccion anterior se aplica a un caso particular, esto es 

[ii,Pj] = ieijkPk- (49) 
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Por lo tanto, Iz y pj tienen un sistema comun de funciones propias. Por otro 
lado, el Hamiltoniano de una particula libre 




por el hecho dc ser el cuadrado de un operador vectorial tiene el mismo 
sistema de FP que I y Iz- Ademas, esto implica que la particula libre se 
puede encontrar en un estado con E, I, m bien determinados. 



Eigenvalores y Eigenfunciones de y L 

Es mas conveniente trabajar con coordenadas esfericas (o polares), dado que, 
como veremos, varios operadores del momento angular actuan solamente 
sobre los angulos variables 9, cp y no en la variable r. En lugar de caracterizar 
el vector r por sus componentes cartesianas x, y, z Uamaremos el punto 
correspondiente M en el espacio (OM. = r) por sus coordenadas esfericas, 
esto es 

X = rcos(l)sm6, y = r siiK^sin^, z = rcos9, (50) 

con 

r > 0, < 6* < TT, 0<(f><2Tr. 

Sean ^{r,9,(j)) y ^>'(r, 9, 4>) las FO de una particula en S y S' en la cual 
la rotacion infinitesimal esta dada por 6a alrededor de z 



^\r,9,(l>) = ^r,9,<f> + Sa), 

I - 

d(t> 



^{r,9,(j)) + 5a^, (51) 



d bien 

$'(r, 9, <^) = (1 + iiz5a)^{r, 9, cj)). (52) 
Luego entonces, resulta que 

9$ . ^ d 

^ = = (53) 

Para una rotacion infinitesimal en x 

= {l + iUa)^r,9,(l>), (54) 
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pero en tal rotacion 

z' = z + ySa; z' = z + y5a; x' = x 

y de (50) se obtiene 

rcos{6 + d6) = r cos9 + r sin 6 sin (j)Sa, 
r sin (f) sin 6 + d6) = r sin^sin^ + r sin^siiK^ — r cos ^(Jo;, 

es decir 

d9 

sin OdO = sin d sin (pSa ^ — = — sin 4>, 

da 

y 



cos 6 sin (f) d6 + sin cos (f)d(f) = —cos 6 5a, 

coscpsmd— = — cosy — cos sin 0—, 
da da 



ahora, sustituyendo (57) en (56) 

d(f> 



da 



= — cot cos (p, 



(55) 

(56) 
(57) 

(58) 
(59) 



tal que al sustituir (56) y (58) en (51 ) y comparando las partes de la derecha 
en (51 ) se obtiene 



d d 
It. = i{ sin 0— + cot cos S— 



En el case de la rotacion cn y, cl resultado es similar, tal que 

{ d d 
ly = i\ sin 0— + cot 6 cos (j)— 
\ oU 0(p 

Usando Ix, ly tambien se puede obtener l±, P, esto es 

/ d d 
l± = exp(ibz0) I ±^ + zcot " 



(60) 



(61) 



d(t>) ' 



n 



1^1+ + P + h, 
1 ^2 



+ 



1 52 



sin^ed(t>'^ sin^ede 



sm( 



(62) 



de (62) se ve que P es identico hasta una constante al operador de Laplace 
en la parte angular, esto es 

1 92 - 



1 d 



dr V dr 



+ 



1 d 

'^ae 



sin( 



dej sin"^ edcj)'^ 



+ 



(63) 
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Funciones Propias de Iz 



L$rr,, = m$ = -i- 



d(j) ' 



= -^e'^t (64) 



Condiciones de Hermiticidad de Iz 

fj rhg dcp = (^j^J g%f + f*gi2n) - /*ff(0). (65) 

Entonces Iz es hermitico en la class de funciones para las cuales 

/*5(27r) = rg{0). (66) 

Las funciones propias de de Ijn pertenecen a (0, 27r) y cumplen con 
(66), asi como para cualquier funcion que se pueda desarroUar en $^((^), 
esto es 

k 

Fi4>) = J2''ke"'t k = 0,±l,±2, ... , 

k 

G{4>) = Y^he""^, fe = ±l/2,±3/2,±5/2 ... , (67) 

solo m enteros 6 m semi-enteros, pero no para comhinaciones de F{(j)), G{(j)). 

Las elecciones apropiadas de m estan basadas en el experimento de FP 
comunes de Iz y r. 

Armonicos Esfericos 

En la representacion {?}, las eigenfunciones asociadas con los eigenvalores 
lil + l)n^ de I? y mh de h son las soluciones de la ec. diferencial parcial 

~i-^xb{r,e,ct)) = mn^{r,e,4>). (68) 

Considerando que los resultados generales obtenidos son aplicables al 
momento angular, sabemos que I es un entero o un semientero y que para 
l,m fijos, pueden tomarse solamente los valores 

-1,-1 + 1, ... ,1-1,1. 
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En (68), r no aparece en el operador diferencial, as/ que puede conside- 
rarse como un parametro y tomar en cuenta solo la dependencia en 6, <p de 
■0. Luego entonces podemos denotar por YimiO^cj)) como una eigenfuncion 
comun de y 1^ la cual corresponde a los eigenvalores de l{l + \)h^,mTi, 
esto es 

L''Yi^{e,c/)) = l{l + l)h^Yim{0,^), 
hYUOA) = mhYUe,^). (69) 

Para ser completamente rigurosos, tenemos que introducir un I'ndice adi- 
cional con el objeto de distinguir entre varias soluciones de (69), las cuales 
correspondan al mismo par de valores I, m. En efecto, como se vera mas 
adelante, estas ecs. tienen solamente una solucion (en un factor constante) 
para cada par de valores permitidos de I, m; esto es porque los submdices 
I, m son suRcientes. 

La ec. (69) did a 6, (p dependencia de las eigenfunciones de y h- 
Una de las soluciones de YimiO, 4>) de estas ecs. ban sido encontradas de la 
siguiente manere^ 

i;Ur,e,^) = f{r)i;Ue,^), (70) 

donde f{r) es una funcion de r la cual aparece como una constante de 
integracion para las ecuaciones diferenciales parciales de (68). El hecho de 
que f{r) sea arbitraria muestra que y no forman un conjunto completo 
de observables comune^ en el espacio o funciones de f (o de r,6,4>). 

Con el objeto de normalizar ipimii^y ^) </*)! conveniente normalizar Yi^iO, 4>) 
y f(r) separadamente (como se muestra). Entonces, debemos tomar un 
diferencial de angulo solido 

r del) r sm9\i;im{e,^)\^de = 1, (71) 
Jo Jo 

POO 

/ r^\fir)\^dr = 1. (72) 
Jo 

^Demostracion en el problema 3.4 

^Por definicion, el operador hermi'tico A es un observable si este sistema ortogonal de 
vectores forma una base en el espacio de estados 

*Cada estado cuantico de la parti'cula es caracterizado por un estado vectorial 
perteneciente a un espacio abstracto Sr 
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Valores de I, m 

a : I, m deben ser enteros 

Usando Iz = j-§^, podemos escrihir (69) como sigue 

^^Yim{e,ct)) = mnYi^{e,(i>), (73) 

la cual muestra que 

YU0,(p) = FUe,cl))e'^t (74) 

Si permitimos que < (p < 2n, entonces podemos cubrir todo el espacio 
ya que la funcion debe ser continua en todas partes, tal que 

Yim{e,cp = 0) = Yim{e, (I) = 27r), (75) 

Jo que implica 

e'"""" = 1. (76) 

Scgun se vio, m es un entero o un semientero, la aplicacion al momento 
angular orbital, muestra que m debe ser un entero. (e^^"^^ sera igual —1 si 
m fuera semientero). 

(3: Todo valor entero (positivo o cero) de I puede ser encontrado escogiendo 
un valor entero de I, se sabe de la teoria general que Yim(9, (p) debe cumplirse, 
esto es 

L+Yim{e,4>) = ^, (77) 
la cual, al combinar = tie^'^ y (62) 

(^■^-lcote^Fu{e) = 0. (78) 

Esta ec. de primer orden puede ser integrada inmediatamente si notamos que 

d(siii 9) 

cotede=— — (79) 
smt' 

su solucion general es 

Fii = ci{smey, (80) 

donde ci es una constante de normalizacion. 

Consecuentemente, para cualquier valor positivo o cero de I, existe una 
funcion Yii{6,(f)) la cual es igual (con un factor constante) 

Y"{e,(f>) = ci{sm9ye'^'f'. (81) 
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A travcs dc la accion repetida de L_ , construimos (p), ... , Yi^ {0, (p) , 

. . . , Yi_i{6, (j)). Luego entonces, vemos que la correspondencia para el par de 
eigenvalores l{l + l)h, mh (donde I es un entero positivo arbitrario o cero y m 
es otro entero tal que I < m < I ), dc (78) una y solamente una cigcnfuncion 
Yim{G,<P)> puede ser ambiguamente calculada de (78). A las eigenfunciones 
Yim{(^:4') se les conocen como armonicos esfericos. 



Propiedades de los Armonicos Esfericos 

a Relaciones de Recurrencia 

Segun los resultados generales podemos tener 



l±Yim{eA) = n^l{l + l)-m{m±l)Yim±i{e,ct>). (82) 

Usando la expresion (62) para l± y el hecho de que Yi^(d, (f)) es el producto 
de una funcion dependiente solo de 6 y e^"^, obtenemos 

e^''*' (^-^-mcote^ YU9,^) = ^Jl{l + 1) - m{m±l)Yi^±i{e,(l)) (83) 
P Ortonormalizacidn y Relacion de Cerradura 

La Ec. (68) determina solamente los armonicos esfericos con un factor cons- 
tante. Ahora eligiremos este factor tal como ortonormalizacidn de Yijn{0, (p) 
(como funcion de variable angular 9, (j)) 

(84) 

Jo Jo 

Ademas, cualquier funcion de 9, (p pueden ser expresadas en terminos de los 
armonicos esfericos, esto es 

oo I 

f{0,(l>)=Yl E QmYlmieA), (85) 
1=0 m=-l 

donde ^ 

cim = r dcp r Sine deYil,{e,^)f{e,cP). (86) 
Jo Jo 

Los armonicos esfericos constituyen una base ortonormal en el espacio 
en de funcioncs de 9, (p. Este hecho se expresa por la relacion de cerradura 

oo I 

J2Y.^lmid,<t>)yi*mid'A) = S{cOSd-COSe')5icP,cP), 



1=0 m=l 

1 



sin^ 



5{e-e')5{4>A)- (87) 
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Es 5{cos 6 — cos 9') y no 6{9 — 9') los cuales entran en el lado derecho de la 
relacion de cerradura porque la integral sobre la variable 9 se efectua usando 
el elemento diferencial sin 9 d9 = —d{cos9). 

Operador de Paridad V 

El comportamicnto de V en tres dimensiones es esencialmente igual que en 
una dimension, es decir, al aplicarlo sobre una funcion ^{x, y, z) de coorde- 
nadas cartesianas solo le cambio el signo, esto es 

Vipix, y, z) = ip{-x, -y, -z). (88) 

V tiene las propiedades de un operador hermitico, ademas es un operador 
unit aria y dc proycccion. 

El operador es un operador identidad 

{r\V\r') = (r|-r') = 6{r + -r'), 

r^\r) =r{r\r) =V\-r) = \r), (89) 

entonces 

V^ = i, (90) 

cuyos valores propios son V = ±00. Ademas se tiene que las FP se Uaman 

pares si V = 00 e inipares siV = —00. En mecanica cuantica no relativista, 
el operador H en un sistema cerrado es invariante ante transformaciones 
unitarias directas 

VHV = V^HV = H. (91) 

Entonces H conmuta con V y consecuentemente la paridad del estado es 
una integral de movimiento. Tambien se cumple para I 

[v,k] = o, [r,i±]. (92) 

Si H es par y imo dc sus eigenestados |$„) cl cual tiene paridad dcHnida 
(V\^n))> colincar a {ipn), se ha encontrado y puede inferirse que el eigen- 
valor correspondiente es degenerado con un grado de degeneration n^, dado 
que V conmuta con H, (7^|$„)) es un eigenvector de H con el mismo eigen- 
valor como |<&n))- Si Tp es FP deV, I y Iz de (92) resulta que las paridades 
de los estados diferentes solo en coinciden. Queda asi determinado la 
paridad de una particula de momento angular I. 



61 



En coordenadas esfericas, para este operador se considera la siguiente 
sustitucion 

r^r, e^TT-e (f>^Tr + (f>. (93) 

Consecuentemente, si usamos una base estandar para el espacio de fun- 
ciones de onda de una particula sin giro, la parte radial de la funcion base 
ipklmi'f^) 6s modificada por el operado paridad. La transformacion solo se da 
en los axmonicos esfericos, como se vera. 

De (93) podemos notar que 

sine sine, cose ^ - cos0 e'"*"^ ^ ^_^■^ni^^m<p ^ ^94^ 

bajo estas condiciones, la funcion Yii{9,(f)) es transformada en 

Yu{cl>-e,7r + cl)) = i-iyYu(e,4>), (95) 

de (95) podemos ver que la paridad de los armonicos esfericos va como (— 1)' . 
For otro lado 

d d d d 

Relacionando (95) y (96) mostramos que l± permanece sin cambio (lo 
cual implica que los operadores l± son pares). Consecuentemente podemos 
calcular Yi^niO, (f)), 

Yi^{(l>-e,7r + 4>) = i-iyYim{e,4>). (97) 

For lo tanto, los armonicos esfericos son funciones cuya paridad esta bien 
definidad e independiente de m, par si I es par e impar si I es impar. 



El Operador Espm 

Algunas paxticulas, ademas de su momento angular tienen un momento 
propio, el cual, sc Ic conoce como espm y denominaremos como S. Este 
operador no esta rclacionado con rotaciones normales espaciales, aim asi, 
cumple con las mismas relaciones de conmutacion que tienen el momento 
angular, esto es 

[Si,Sj]=iSijkSk- (98) 
Asi como, las siguientes propiedades 

(1). Para el espm, valen todas las formulas de (23) a (48) del momento 
angular, las cuales son similares a (98) 



62 



(2) . El espectro de la proyeccion del espm, es una secuencia de numeros 

enteros 6 semienteros que diReren por una unidad. 

(3) . Los valores propios de son 

SV = 5(S + 1)V'.. (99) 

(4) . Para un S dado, la componente Sz solo puede tomar 25 + 1 valores, 

de -S a +S. 

(5) . Las FP de las particulas con espm, ademas de depender de f 6 p, 

dependen de una variable discreta (propia del espm ) a, la cual denota 
la proyeccion del espm en z. 

(6) . La FP il){f, a) de una particula con espm se puede desarroUar en FP 

con proyecciones dadas del espm Sz, esto es 

s 

^(r,a)= Yl V'a(r)x(^), (100) 

a=-S 

donde ipai^ es la parte orbital y x(<t) es la parte espinorial. 

(7) . Las funciones de espm (espinores) x(o'i) son ortogonales para cualquier 

par dc di ^ (Tfc. Las funciones V'o-(r)x(o') se Jes conoce como las com- 
ponentes de las FO de una particula con espm. 

(8) . La funcion if^aif) se llama parte orbital de la FO 6 solo orbital. 

(9) La normalizacion se hace como sigue 

E IIV'.(r-)|| = l. (101) 

(7=-S 

Las relacioncs dc conmutacion permiten establecer la forma concrcta dc 
los operadorcs (matrices) de espm que actiian en el espacio de las FP del 
operador proyeccion del espm. 

Muchas particulas elementales tales como el electron, el neutron, el 
proton, etc. tienen espm 1/2, por eso la proyeccion toma solo dos valores, es 
decir = ±1/2 (en unidades h). 
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For otro lado, las matrices Sx,^y, Sz en el espacio de las FP de 5^, Sz son 



Definicion de las Matrices de Pauli 

Las matrices 

(Ji = 2Si (103) 

se llaman matrices de Pauli, las cuales son matrices hermitianas, tienen la 
misma ec. caracteristica 

- 1 = 0, (104) 
par consiguiente, las eigenvalores de ax, cry y cTz son 

\ = ±1. (105) 

For lo tanto, son consistentes con el hecho de que Sx, Sy y Sz sean iguales ail. 
Ademas 

I 

En el caso para el cual un sistema con espm tiene simetria esferica 
(esferico simetrico) 

Mr,+k), Mr,-l), (107) 

son diferentes soluciones por la proyeccion Sz- 

El valor de la probabilidad de una u otra de las proyecciones esta deter- 
minada por el cuadrado de ||'0i,2|P de tal modo que 

||Vi|P + ||V2|P = l. (108) 
Como las FP de Sz tiene dos componentes, entonces 

XI = ( J ) , X2 = ( J ) , (109) 
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tal que, las FP de una particula de espm 1/2 sc puede escribir como 

V' = V'ixi + te= (| j . (110) 

A continuacion las orbitas van a ser sustituidas por numeros dado que 
solamente la parte del espm es importante. 



Las Transformaciones a las Rotaciones 

Sea la FO de un sistema con espm en E. Encontremos la probabilidad de la 

proyeccion del espm en una direccion arbitraria en el espacio tridimensional 
(3D) que la toma como eje z' de Tl . Como ya se vio, se tienen dos metodos 
para su solucion 

a) ^ no cambia cuando T, T,' y cl opcrador A cambia como un vector. 
Debemos encontrar las FP de la proyeccion S'^ y desarroUamos en 
esas FP. Los cuadradados de los modulos de los coeGcientes dan el 
resultado. 

5; = cos (j) + Sysin(f) = e-'^'^S^e'^'^, 
S'y = - sin ,P + Sy cos (p = e''^'^ Syc'^'^ , 

S', = -S, = e'^'f'S,, (111) 

con rotaciones inGnitesimales y de las relaciones de conmutacion se 
puede encontrar 

L = S„ (112) 

donde L es el generador. 

P) La segunda representacion es: 

S no se cambia a Ja S — > S' 7 Jas componentes de ip se cambian. La 
transformacion a esta representacion se hace con una transformacion 
unitaria 



V^S'V = A, 




(113) 
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de (111) y (113) resulta que 

V^^-iSz't'Se'^-'I'V = S, 

de (114) se obtiene 

Usando la forma concreta de Sz y las propiedades de las matrices de 
Pauli se obtiene la forma concreta , tal que 

Un Resultado de Euler 

Cualquier sistema de referenda S' de orientacion arbitraria con respecto a 
S puede ser alcanzado con solo tres rotaciones, primero alrededor del eje z, 
en seguida una rotacion del angulo 9 sabre el nuevo eje de coordenadas x' y 
por ultimo el angulo ipa en el nuevo eje z' . 

Los parametros (if, 9, ipa) se les llama angulos de Euler 

VH^,e,^a) = Vl{4^a)Vi,{9)V^{v). (117) 

Las matrices son del tipo de (116), en cuanto a es de la forma 

^^(^) = f '"'I '""IV (118) 

de tal modo que 

e 2 COS I le 2 sin| \ , . 

. . -i£±la I ■ (119) 

le 2 sm I 6 2 cos | / 

Entonces por la rotacion de S, las componentes de la funcion espinorial 
se cambian como sigue 

I :] •P + 'i'a ,: Va-y 9 

ipi = tpie 2 cos - + 1^/^26 2 sm-, 

,; . , ; V-^a . 6 , ,■ y+V-a 9 /-, 

tp2 = i'fpie 2 sin- + '02e 2 cos-. (120) 
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De (120) se puede ver que para una rotacion en E3 le corresponde una 
transformacion lineal en E2 — espacio Euclidiano bidimensional (2D) — las 
dos componcntes de la funcion espinorial. La rotacion en no implica una 
rotacion en E2, la cual signiGca 

i^'liP') = {^\^) = + ^;V'2. (121) 

De (119) se encuentra que (121) no se cumple, sin embargo, hay una 
invariancia en las transformaciones (119) en el espacio E2 de las funciones 
espinoriales, el cual es 

= '01$2 - '02^1- (122) 

Las transformaciones lineales que dejan invariantes tales formas bilin- 
eales se llaman binarias. 

Una transformacion fisica con dos componentes para la cual una rotacion 
del sistema de coordenadas es una transformacion binaria se llama espm de 
primer orden o solamente espm. 



Funciones de Onda Espinoriales de un Sistema con 2 Fermiones 

Las funciones propias de is^ — i = 1, 2 — tienen la forma siguiente 

iH^(l). (123, 

Una variable — o mejor dicho, un operador — corriente en un sistema de 
dos fermiones es el espm total 

S=iS+2S (124) 

Las funciones propias espinoriales de 1^ Sz son los kets \S,a), las cuales 
son combinaciones lineales de las FP de iSz, esto es 

Las funciones de (125) son ortonormalizadas. En En el estado | + +) es 
Sz = 1 y al mismo tiempo es funcion propia del operador 

S=is^ + 2(is){2s)+2S^- (126) 
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Como se puede ver de 

= \ + +) = ^\ + +) +2{lS^ -2 3^+1 Sy -2 Sy+lS^ -2 s,)\ + +)(127) 

= I + +) = 2| + +) = 1(1 + 1)1 + +). (128) 
Si se introduce el operador 

S- =1 s_ +2 (129) 

se obtiene que 

[S-,S'^] = 0. (130) 

Entonces (S'_)'^|l, 1) se puede escribir en funcion de las FP del operador S^, 
esto es 

5_|1, 1) = 5_| + +) = V2\ + -) + V2\- +). (131) 

Resulta que Sz = en el estado S-\l, 1). For otro lado, de la condicion de 
normalizacion tenemos 

|1,0) = ^(| + -) + !-+)) (132) 
5_|1,0) = !--) + I --) = a|l,-l). (133) 

De la coordenada de normalizacion 

|1,-1) = |-,-). (134) 

Solo hay una combinacion lineal independiente mas de funciones de (125) 
diferentes de |1, 1), |1,0) y |1, —1), esto es 

^4 = ^(1 + -)-!-+)), (135) 

Szip4 = 0, S'^ipi. (136) 

Por consiguiente 

V'4 = |0,0). (137) 

ijj4 describe el estado de un sistema de dos fermiones con el espm total igual 
a cero. Este tipo de estado se llama singlet. Consecuentemente el estado de 
dos fermiones de espm total igual a uno se llama triplet y tiene un grado de 
degeneracion g = 3. 
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Momento Angular Total 

Se define como la suma del momento angular orbital mas el spin, esto es 

J = i+S, (138) 

donde, I y S como hemos visto, actuan en espacios diferentes, pero el 
cuadrado de l y S conmutan con J, es decir 

[J^,Jj]=ieijkJk, [Ji,P] = 0, [Ji,S^]=0, (139) 

(139) que r y S*^ tienen un sistema de FP con J^, y Jz- 

Encontramos el espectro de las proyecciones de Jz para un fermion. el 
estado de proyeccion de maximo de Jz se puede escribir como 

^ = i^ii(^l^=\l,l,+) (140) 

= (/ + i)v^,^ j = / + i (141) 

Si introducimos el operador J_ definido por 

J_ = L + S_=L+(5 l). (142) 



De la normalizacion a = a/( J + M)(J — M + 1) se obtiene 
J.^u ^ J j =^|z,Z-i,+) + |/,Z-i,-), 



(143) 



tal que el valor de la proyeccion de j- en j^il) sera 

jz = {l-l) + \ = l^\, (144) 

resulta que j_ disminuye por una unidad a Jz- 
En el caso general tenemos que 

f_ = t + kl'l-^S^, (145) 

se observa que (145) se obtiene del desarroUo binomial si se considera que 
s?. y todas las potencias superiores de s son cero. 

= t\l,l,+) + kt-%l,-). (146) 
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Sabemos que 

tal que al usarla se obtiene 

i-l^' +) = /IHl^' l-k,+) + l-k + l,-), (147) 

con la notacion m = I — k 

= ^l^^^\l,^,+) + ^[^^^il-m)\l,m+l,-). (148) 

Los valores propios de la proyeccion del momento angular total es la secuen- 
cia de numeros que diGeren por la unidad desde j = I + \ hasta j = I — \. 
Todos estos estados pertenecen a la misma funcion propia de J como \l,l,+) 
porque [Ji, J^] = C 

J^\l,l,+) = iP + 2lS + S^)\l,l,+), 

= [Z(Z + l) + 2Zi + |]|/,/,+) (149) 

dondej{j + l) = {l + ^){l + l). 

En la dcrecha dc (149) una contribucion diferente de cero da solamente 
j = IzSz- Entonces las FP obtenidas corresponden a j = I + rrij = m + ^. 
Las FP son de forma 

El numero total de estados lineales independientes es 

N = {2l + l){2s + l) = Al + 2. (151) 

El sistema de FP constituido de tal manera contiene 2j + 1 + 21 + 1 
estados 

\l-hm-l) = m, +) - y^±^|l, m + 1, -). (152) 

Si dos subsistemas estan interaccionando de tal mancra que cl momento 
angular de cada Ji se conserva, entonces las FP del operador momento 
angular total 

i = ii + j2, (153) 
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sc pucdc cncontrar como lo hicimos anteriormente. Para valores propios de 
h y h hay 

(2ii + l)(2i2 + l), 

FP ortonormalizadas dc la proyeccion del momento angular total Jz- La FP 
que corresponde al valor maxlmo de la proyeccion Jz es deck, 

Mj =jl +j2, 

se puede construir de manera unica y por lo tanto J = ji + j-i cs cl valor 
maximo del momento angular total del sistema. Aplicando el operador J = 
Ji + j2 de manera repetida a la funcion 

\ji +32,h + h,h +32} = \ juh) ■ \j2,j2). (154) 

se obtienen todas las 2{ji + 22) + 1 funciones ortogonales de la FP de J = 
ji + J2 con varios M 

-(J1+J2) < M< (ji+j2). 
For ejemplo las FP para M = ji + j2 — 1 es: 



/ Ji I 32 

|jl+j2, jl+j2-l, Jl,i2) = \ . , . jl-l,j2,j2) + W ■ , ■ jl,j2,i2-l) 

\jjl+j2 V-?l+-?2 

(155) 

Aplicando en seguida vaiias veces el operador J_ se puede obtener las 
2{ji + i2 — 1) — 1 funciones de J = ji + ^2 — 1- 
Se puede demostrar que 

\ji -j2\< J < ji +32 

tal que 

max J 

J2 (2J + 1) = (2Ji + 1)(2J2 + 1), (156) 

min J 

Consecuentemente 

\J,M,ji,j2)= {3i,m,j2\J+)\ji,mi,j2,m2) (157) 

mi+m2=M 

Citas: 1. Acetatos del Prof. H. Rosu 
Referenda bibliografica: 

1. H.A. Buchdahl, "Remark concerning the eigenvalues of orbital angular 
momentum", 

Am. J. Phys. 30, 829-831 (1962) 
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Problemas 
Problema No. 3.1 

Mostrar que si ip' = Rip, entonces el operador R se puede representar como 

un operador exponencial 

Solucion 

Para mostrarlo, expandemos ip' (r) en serie de Taylor en el punto x' = x + dx 
y considerando solo las primeras potencias tenemos 

d 

i^ix' ,y' ,z') = ip{x,y,z) + {x' - x)—tp{x',y',z 

d 

+{z'-z)-^iPix',y',z') 



ahora considerando el hecho de que 



x' = X - yd(f>, y' = y + d(j), z' = z, 
tal que, esto reduce la serie de tres dimensiones a solamente dos 

il){r') = 'ijj{r) + {x — yd(/) — x) — \- {y + xd(j) — y)- 



dx 

= iP{r)-dct)y^^-d(l>x^^, 



dy' 



1 



d__ d_ 
dy ^ dx 



considerando que ilz = (2;^ — y-§^ entonces tenemos que R = 
entonces R puede escribirse como una exponencial 



— dcj) 



dy 



R = e 



ilzd4> 



Problema No. 3.2 

Mostrar que de las expresiones dadas en (14) se puede Uegar a (15) 
Solucion 
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Nuevamentc, consideremos solo los terminos lineales en la expansion de la 
serie de Taylor y dado que tenemos rotaciones infinitesiniales, entonces 

/l.d^ = 1 + ildcp + ^{iU(Pf + ... , 

entonces tenemos que 

(1 + ilzd(j))Ax{l — ilzdcj)) = Ax — A^dcp, 
{Ax + iizd(j)Ax){l - ilzdcj)) = Ax- Axdcf), 
Ax - Axilzdcf) + ilzd(j)Ax + lzd(f)Axlzd(f) = Ax- Axd(f), 
i(lzAx — Axlz)d(j) = —Ayd(f). 

Luego entonces, concluimos que 

[^Z) Ax] = lAy. 

entonces tenemos que 

(1 + ilzd(f))Ay{l — ilzdcj)) = Axdcj) — Ay, 
{Ay + ilzdcj)Ay){l — ilzdcj)) = A^dcj) — Ay, 
Ay — Ayilzdcj) + ilzdcj) Ay + lzdcj)Aylzd<j) = Axdcj) — Ay, 
i{lzAy — Aylz)d(j> = —Axdcj). 

Problema No. 3.3 

Determine la precesion del espm en un campo magnetico homogeneo. 
Solucion 

Si cl cuerpo cargado se mueve en un campo magnetico homogeneo circular 
alrededor de la direccion del campo con una frecuencia 

u = 2ul = • 

mc 

Aqui, la carga del electron es —e. Esto se sigue del hecho de que la fuerza 
de Lorentz equilibra la fuerza centrifuga. 

eBv 2 

= myu , 

c 

entonces 

mc 



73 



Por lo tanto 



eB 



UJL = 



2mc 



la cual, se conoce con el nombre de frecuencia de Larmor. 
Problema No. 3.4 

Resolver la ec. de Laplace usando coordenadas esfericas 
Solucion 

Asumiendo que podemos tener 



U{r,e,ct>) = R{r)Q{e,<i)) 



de este modo vemos que 



r d'^ 
R{r) dr"^ 



[rR{r)\ = ^L^e = l{l + 1) 
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4. El Metodo WKB 



Para estar en condiciones de estudiar los efectos de potenciales mas rea- 
listas, que los correspondientes a barreras y pozos de potencial, es necesario 
cncontrar mctodos que permitan resolver la ecuacion de Schrddinger para 
dichos potenciales. 

En general no es posible construir soluciones exactas para tales casos, 
y lo que se Lace, es recurrir a metodos de aproximacion que proporcionen 
una solucion suficientemente huena y simple, como para poder estudiar el 
comportamineto del sistema con ella. 

Metodos como estos hay muchos, pero nos concentraremos en el metodo 
desarroUado simultaneamcnte por G. Wentzel, M. A. Kramers y L. Brillouin 
en 1926; y de cuyos apellidos deriva el acronimo WKB. 

Es importante mencionar que el metodo WKB, solo es aplicahle a la 
ecuacion de Schrddinger 1-dimensional. 

Para resolver la ecuacion de Schrddinger 



2m dy'^ 



+ u{y)'ip = Elf) 




supongamos que el potencial tiene la forma: 




(2) 



y hacemos los cambios de variables: 



e = 



(3) 



E 



(4) 



y 



(5) 



X 



a 



de la ecuacion (5) obtenemos: 



d dy d 
dx dx dy 



d 



(6) 



SP__d_^ d_ 
dx^ dx ^ dy 




(7) 
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J la ecuacion de Schrddinger se escribe: 



2^ ip ^ 



(8) 



multiplicdndola por —1/^^ y deGniendo r{x) = rj — f{x), es posible escribirla 
en la forma: 

9^ + ^r{x)i; = (9) 
para resolver esta ultima, proponemos la siguiente solucion: 



ip{x) = exp 



q{x)dx 



(10) 



En general: q{x)dx = Q{x)\l = Q{x) - Q{a) 9 ^ = = q{x). 



Esto de acuerdo con el teorema fundamental del calculo. 
Por lo que: 



dx'^ dx\dx J dx\^^^^^ 



I 



q{x)dx 



Factorizando ijj tenemos: 

_ 

dx"^ 



- / q{x)d2 

S Ja 



dq{x) 
dx 



^ ^ e dx 



exp 



I 



q{x)dx 



(11) 



Olvidandonos de la dependencia en x, la ecuacion de Schrddinger se 
escribe ahora: 

1 o i da 1 



^1 



En general V' 7^ por lo que: 



^dx e 



i^^+r-q' = 

dx 



(12) 



(13) 



que es una ecuacion diferencial lineal tipo Riccati, cuya solucion se busca 
como una serie de potencias de ^; suponiendo que ^ es muy pequeHa. 
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Dicha serie proponemos que tiene la forma: 

oo 

q{x) = 5](-iO"g„(x) (14) 
n=0 

Sustituimos esta en la Riccati para obtener: 

oo 1 oo oo 

^ + rix) - Y: i-^(r% E (-^0^3. = (15) 

n=0 "'^ n=0 u=0 

Rearreglando los terminos de la ecuacion (15) tenemos: 

oo 7 oo oo 

E(-1)"(^0"+'-^ + r(x) - E E(-^0'^+'^^m5. = (16) 

n=0 /i=0 1^=0 

Las series dobles cumplen con: 

oo oo oo n 

E E ^1^'^ ~ E E ^'m,n-m 
ix=0 i'=0 n=0 m=0 

Donde: ji = n — m u = m 
De esta forma: 

oo 7 oo n 

E(-l)"(^0"+'^ + rix) -EE (-^^"-"^^"'^mgn-m = (17) 

n=0 n=0 m=0 

Veamos por separado unos cuantos terminos de cada una de las series 
contenidas en la ecuacion (17): 



v(-i)-(zO"+i^ = + - +... (18) 

' ^ Hif rlif rlnc 

n=0 



E E i-iO'"Qmqn-m = Qo - i^CqoQi + ■■■ (19) 
n=0 m=0 

Para que ambas series contengan en su primer termino, de sus respectivos 
desarroUos a i^, debemos escribirlas: 

oo 7 oo n 

E(-1)""'(^0"^ + r{x) - - E E i-^Crimqn - m = 

n=l n=l m=0 
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De la cual obtenemos: 



E 

n=l L 



He^^-E(-^0^m9n-. 



m=0 



+ 



r(x) - ql 



= (20) 



Para que se satisfaga la igualdad anterior, debemos exigir que: 

fix) -ql = =^ qo = ±\fr(x) 

dx 
dqn- 



(21) 



-(-iO"^ - E i-iO^'qmqn-m = 
m=0 



dx 



- = - E ^'"5"- 



n > 1 



(22) 



m=0 



A la cual Uamaremos relacion de recurrencia. Recordando que deGnimos 
r{x) = 1] — f{x), 11=^ ^ fix) = ^ ; obtenemos con ayuda de la 
ecuacion (21) que; 



qo = ± Jr? - fix) 



I E _ u 



l2miE-u) 
2muo 



(23) 



Esta ultima es la conexion clasica para el momento de una particula de 
energia E en el potencial u, en unidades de ^/2muo. For ello: 



qo = pix) = ^Jt]- fix) 



(que no es operador) Si aproximamos hasta segundo orden, tenemos lo 
siguiente: 

qix) =qo- i^qi - fq2 
y empleando la ecuacion de recurrencia calculamos qi y q2: 



dqo o 

— = -2qoqi 
ax 



2 go 2 ax 



qi 



1 d 
2dx 



(In \pix)\) 



dqi 
dx 



-2qoq2 - qi 



q2 



2qo 



(24) 
(25) 
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De la ecuacion (24), nos percatamos de que qi es la pendiente con el signo 
cambiado de In \qo\; cuando qo es muy pequeno, gi <C =^ —^Qi ^ 0, 
y en consecuencia la serie diverge. Lo que nos Ueva a exigir la siguiente 
condicion WKB; 

kol > I - ^Qi\ = ^kil 

La condicion WKB no se satisface para puntos Xk tales que: 



Qo(,Xk) =p{xk) = 



recordando que qo = p 



2m{E—u) r ■ ' J. ■ J 

2mu^^ ? ccuacion anterior nos conduce a: 



E = u{xk) 



(26) 



En fisica clasica puntos Xk que satisfacen la ecuacion (26), son Uamados 
puntos de retorno; ya que en ellos la particula invierte el sentido de su 
movimiento. 

En base a lo anterior, podemos decir que qo es una solucion clasica del 
problema, y que qi y q2 son respectivamente, la primer y segunda correc- 
ciones cuanticas del problema. 

Para obtener las funciones de onda, solo consideraremos la solucion 
clasica, y la primer correccion cuantica del problema; y las ssutituimos en 
nuestra propuesta para ip: 



i/j = exp 



i 




- j q{x)dx 


= exp 


S Ja 





V' = exp y qodx^ ■ 



i. Ja 

qidx 



exp 



Para el segundo factor tenemos: 



exp 



qidx 



exp 



exp 



i/;A(ln|pM|)cfa 



- 2 (111 \Pi^)\) 



A 



con A una constante. Para el primer factor tenemos: 



exp 



- J qodxj = exp ^ J P{^)dx 
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y ip puede ser escrita como: 



exp 



p{x)dx 



(27) 



y se Uaman las soluciones WKB de la ecuacion de Schrodinger uno- 
dimensional. 

La solucion general WKB en la region para la cual se cumple la condicion 
WKB, se escribe: 

i) = a+V^ + a-i^" (28) 

Como ya se dijo, no hay solucion WKB en los puntos de retorno; lo que 
nos lleva a cuestionarnos como es que ip{x < x^) pasa a ip{x > x^), y para 
esto se hace necesaria la introduccion de las formulas de conexion. 



Las Formulas De Conexion 

Ya se dijo que las soluciones WKB, son singulares en los puntos clasicos 

dc retorno; no obstante cstas soluciones son validas a la izquicrda, y a la 
derecha dc un punto clasico dc retorno x^. Y par ello nos cuestionamos 
como es que tp{x < Xk) pasa a ■iIj{x > Xk); es decir, debemos encontrar las 
formulas de conexion. 

De la teoria de ecuaciones diferenciales ordinarias, y con apoyo del 
analisis de funciones de variable compleja, puede demostrarse que las formulas 
de conexion existen y que son las siguientes: 



ipi{x) = — rexp (- [ J-r{x)dx 



donde ipiix) solo tiene un comportamiento exponencial decreciente para 
X < Xk- Lo que significa nuestra primer formula de conexion, es que una 
funcion il^{x), que a la izquierda de un punto clasico de retorno se comporte 
como una exponencial decreciente, pasa a la derecha del punto clasico de 
retorno como un coseno de fase (l> = j, y con el doble de la amplitud de la 
exponencial. 
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Ahora, en el caso de una funcion il){x) mas general; es decir, una funcion 
que tenga un comportamiento exponencial creciente y decreciente; la formula 
de conexion correspondiente es: 

sin ((p + — T- exp ( [ J-r(x)dx] 

< 3-^ cos ( r Jr{x)dx + (30) 

[r(x)]4 V-'^^fc / 

siempre que (p no tenga un valor muy cercano a —j; la razon de ello es que 
si (f) = — J, el seno se anula. Esta segunda formula de conexion, significa 
que una funcion que se comporta como un coseno a la derecha de un punto 
cMsico de retorno, pasa a la izquierda de el como una exponencial creciente 
con amplitud modulada por un seno. 

Para ver los detalles de como son obtenidas estas formulas de conexion, 
dehe consultarse el libro: Mathematical Methods of Physics by J. Mathews 
& R.L. Walker. 



Estimacion Del Error Introducido en la Aproximacion WKB 

Hemos encontrado la solucion a la ecuacion de Schrddinger en cualquier 
region donde se satisfaga la condicion WKB; no obstante, las soluciones 
WKB divergen en los puntos clasicos de retorno como ya se ha senalado. 
Analizarcmos un tanto superficialmente esta problematica a fin de proponer 
las Uamadas formulas de conexion en una vecindad proxima a los puntos 
clasicos de retorno. 

Supongamos que x = x^, es un punto clasico de retorno; es decir, es un 
punto tal que se cumple: qo{xk) = p{xk) = ^ E = u{xk). Ahora 
bien, a la izquierda de x^; es decir para puntos del espacio 1-dimensional 
tales que x < xj-, supongamos que E < u{x), de modo que en esta region la 
solucion WKB es: 



tl){x) = 



a 



~ u{x)-E 



exp 




+ 



81 



+ 



u{x)-E '[ 4 
uo 



b 




1 u{x) -E 

1 Jx ^ Uo 



dx 



) 



(31) 



de igaal forma a la derecha de Xk, es decir para puntos del espacio 1- 
dimensioanl tales que x > x^ supondremos que E > u{x), en consecuencia 
la solucion WKB en esta region es: 



Si ip(x) es una funcion real, lo sera tanto a la derecha como a la izquierda 
de Xk, a esto le Uamaremos "reality condition", y establece que si a, 6 G 3fJ, 
entonces c = d* . 

Nucstro problema es conectar las aproximaciones a cada lado de x^ a 
modo de que se re&eran a la misma solucion exacta; esto es encontrar c 
y d si conocemos a y b, y viceversa. Para hacer dicha conexion, debemos 
utilizar una solucion aproximada, la cual sea valida a lo largo de un camino 
que conecte las regiones a ambos lados de x^, donde las soluciones WKB 
sean validas tambien. 

Lo mas comun es recurrir a un metodo propuesto por Zwann y Kemble, 
el cual consiste en evadir el eje real en las cercanias de x^, mcdiante el 
recorrido de un camino que encierre a Xk en el piano complejo; a lo largo de 
este camino las soluciones WKB seguiran siendo validas. En esta exposicion 
recurriremos a dicho metodo, pero solo con la Bnalidad de obtener un medio 
de cstiniar errores en la aproximacion WKB. 

La cstimacion de errores es importante, a causa de que se desca obtener 
soluciones aproximadas, en un amplio intervalo de puntos del espacio 1- 
dimensional; y se debe estar preocupado en si el error se acumula, y si 
poster iormente traera consigo corrimientos dc fase. 

Para esto deGnimos las funciones WKB asociadas como: 




+ 



+ 




(32) 




(33) 
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a estas las consideraremos como funciones de variable compleja respecto de 
X, y emplearemos cortes de tipo rama para eludir discontinuidades en los 
ceros de r(x) = Estas funciones satisfacen una ecuacion diferencial, 

que puede obtenerse diferenciandolas respecto a x, para tener: 



w4 



I Ir 
i 4 r 



r 1 r" 



r 

16 \ r 



nombramos a: 

, ^ Ir" 5 

six) = 

^ ^ 4 r 16 

entonces las W± son soluciones exactas de: 

1 



w'l + 



,^r{x) + s{x) 



W± = 



(34) 
(35) 

(36) 



mientras que solo satisfacen aproximadamente a la ecuacion de Schrddinger; 
la cual es regular en x = Xk, mientras que la ecuacion satisfecha por las 
funciones WKB asociadas es singular en dicho punto. 

Procedamos a deRnir funciones a±{x) tales que cumplan con las dos 
relaciones siguientes: 



tp{x) = a-i-{x)W+{x) + a-{x)W-{x) 

^'(x) = a+{x)W'^{x) + a-{x)W'_{x) 



(37) 

(38) 



donde tp{x) es una solucion a la ecuacion de Schrddinger. Resolviendo las 
ecuaciones anteriores para las a±; tenemos: 



W+W'_ - W'^W. 



a_ = — 



siendo el denominador de estas el Wronskiano de y W-; no es dificil 
demostrar que este toma el valor — |i, as/ que las a± se simpHHcan a: 



a- 



(39) 
(40) 
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Tomando la derivada respecto a x de las ecuaciones (39) y (40), tenemos: 

(41) 



dx 2" 

dentro del parentesis el primer y cuarto termino se anulan; recordemos que: 



iP" + ^r{x)il; = & + 



podemos escribir la ecuacion (41) como: 

da± ^ 



^r{x) + s{x) 



W± = 



dx 



2 

= T^isix)ip{x)W^{x) 



y en base a las ecuaciones (33) y (37); 



da± 
dx 



^ [r(x)]2 



a± + q;=p exp ^ =F^^ 



r{x)dx 



(42) 



(43) 



Las ecuaciones (42) y (43) son usadas para estimar el error que se 
comete en la aproximacion WKB para un punto particular del espacio 1- 
dimensional. 

La razon de que se considere a como una estimacion del error que 
se comete en la aproximacion WKB, es que cn las ecuaciones (31) y (32) las 
constantes a, b y c, d respectivamente, tan solo nos dan soluciones ijj apro- 
ximadas; mientras que las funciones a± al introducirlas en las ecuaciones 
(37) y (38), nos proporcionan soluciones ip exactas; y al tomar su derivada 
obtenemos la pendiente de la recta tangente a ellas, y esta nos dice cuanto 
es que se desvian las a± de las constantes a, b, c y d. 

Nota: Los articulos WKB originales son: 

G. Wentzel, "Eine Verallgemeinerung der Wellenmechanik" , 
Zeitschrift fiir Physik 38, 518-529 (1926) [recibido 18 Junio 1926] 



L. Brillouin, "La mecanique ondulatoire de Schrddinger; une methode generale 
de resolution par approximations succcssivcs" , 

Compte Rendue Acad. Sci. (Paris) 183, 24-26 (1926) [recibido 5 Julio 1926] 
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H.A. Kramers, "Wellenmcchanik und halbzahlige Quantisierung" , 
Zf. Physik 39, 828-840 (1926) [rccibido 9 Sept. 1926] 

H. Jeffreys, "On certain approx. solutions of linear diff. eqs. of the second 
order", 

Proc. Lond. Math. Soc. 23, 428-436 (1925) 



85 



Problema 4.1 



Problemas 



Veamos un ejemplo dc como sc usa cl mctodo WKB cn mccanica cuantica: 
Consideremos unaparticula de energia E que se mueve en un potencial u{x), 
la correspondiente ecuacion estacionaria de Schrddinger es: 

g + |j[^-.(x)]^ = (44) 
y considerenaos que u{x) tiene la forma que se muestra en la figura 4.1. 



u(x) 



Fig. 4.1 



Como podemos ver: 



, , 2m . ^ , , , I es positiva para a < x <b 
r(x) = \E — u(x)\ < , 

I es negativa para X < a,x > b 



Si iIj{x) corresponde a puntos tales que x < a, al pasar al intervalo 
a < X < b, nuestra formula de conexion es la ecuacion (29) y nos dice que: 



4 



[E-u 

donde A es una constante arbitraria, 




i'ix)^— ^cos(/ \l-^{E-u)dx--\ (45) 
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Cuando iIj{x) corresponde a x > b, al pasar al intervalo a < x < b 
similar mente: 



ip{x) ^ — - 



B 



[E-u] 



cos 




— rrl-C/ — u)ax — - 



(46) 



donde B es una constante arbitraria. La razon de que nuestra formula de 
concxion sea nuevamente la ecuacion (29), es que cuando la particula Ucga 
al segundo punto clasico de retorno, x = b, esta invierte la direccion de su 
movimiento, y entonces es como si hubiera venido de derecha a izquierda; 
lo que equivale a ver la primer situacion de izquierda a derecha en el punto 
X = a, cn un espejo. 

Estas das expresiones independientemente de las constantes A y B, 
deben de ser las mismas; as/ que; 



cos 





(47) 



Recordando que: 



/A + B\ fA-B 
cos A + cos B = 2 cos ( — - — I cos 



la ecuacion (47) se escribe: 
2 cos 




iE-u)dx-^ + J^{E-u)dx-^ 



cos 



I j 1 2m TT 1 2m vr 



= (48) 



lo que implica que los argumentos de estos cosenos sean un multiplo entero 
de |; el argumento del segundo coseno no nos Ueva a algun resultado intere- 
sante, por lo que solo prestaremos atencion al argumento del primer coseno, 
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el cual por el contrario si nos Ueva a un resultado de gran importancia; 
entonces: 




Este resultado es muy similar a las reglas de cuantizacion de Bohr - 
Sommcrfcld. 

Recordemos que cl postulado de Bohr establece que el momento an- 
gular de un electron, que se mueve en una orbita permitida en torno a 
un nucleo atomico, esta cuantizado y su valor es igual a: L = nh, n = 
1,2,3,.... Y recordemos tambien que las reglas dc cuantizacion dc Wil- 
son - Sommerfeld, establecen que toda coordenada de un sistema fisico que 
varie periodicamente en el tiempo debera satisfacer la condicion cuantica: 
§ Pqdq = righ; donde q es una coordenada periodica, pg es el momento asoci- 
ado con ella, Ug es un numero entero y h es la constante de Planck. Y como 
podemos ver el resultado obtenido de la aproximacion WKB es muy similar 
a estos dos. 
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Problema 4.2 



Estimemos el error que se comete en la solucion WKB, en un punto 
xi 7^ Xk, con Xk un punto clasico dc rctorno; para la ccuacion difcrcncial 
y" + xy = 0. La solucion de este problema en fi'sica, es importante para el 
estudio de campos uniformes; tales come el campo gravitacional 6 el campo 
electrico uniforme debido a placas planas con carga electrica. 
Solucion: 

Para esta ecuacion diferencial tenemos que: 

5 

^ = 1, r{x) =x & s{x) = - J^x''^ 

r{x) = X solamente tiene un cero en x^ = 0, as/ que para x 3> 0; 

W± = x^s exp j ^/xdx^ = x"* exp ^zb^ixz^ (50) 

Derivando las W± una y dos veces respecto a x, nos damos cuenta de que 
satisfacen la siguiente ecuacion diferencial: 

+ (x - ^x-^)W± = (51) 

La solucion exacta y{x) a esta ecuacion diferencial, la escribiremos conio 
una combinacion lineal de las W±, tal y como se indico en la seccion corre- 
spondiente a la estimacion de error en la aproximacion WKB; si recordamos 
la combinacion lineal se propuso de la forma: 

y{x) = a+(x)W+(x) + a_(x)W_(x) 

Para x muy grandes una solucion general de nuestra ecuacion diferencial, 
esta descrita por la aproximacion WKB como: 

y{x) = Ax'i cos (^^^ + cuando x —>■ oo (52) 

de modo que a_|_ — > ^e^^ y ol- ^ 4^"*'^ para x — > oo. Deseamos calcular 
el error en esta solucion WKB; el cual es medido por la desviacion de a+ y 
a- respecto de las constantes A. Para esto utilizamos la ecuacion: 



da± ^ s(x) 



dx 2 -^r(x) 



a± + Q!=p exp (^2i J ^J r{x)dx 
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Y efectuando las sustituciones correspondientes: 



da± 
dx 



=F- X 



_ i_ 

X 2 



2 2 ^13 



(53) 



sabemos que Aa±, representa los cambios que sufren las a± cuando x va 
desde xi hasta oo, y estos cambios se calculan mediante: 



Aa± 
A/2 



2 
A 



da± 
dx 



dx 



ibi — e 
32 



9 3 /"OO 



us I 

J XI 



X 2 exp ( =F^-x2 



dx 



(54) 



El segundo termino dentro del parentesis es menos importante que el primero, 
esto se debe a que la exponencial compleja, oscila entre 1 y —1 y x~^ < x~^ . 
de este modo: 

Aa± , 5 



A/2 
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X, 



(55) 



y como podemos ver el error que se introduce es realmente pequeno, esto 

_3 

porque igualmente la exponencial compleja oscila entre —lyljyx-^"^ sera 
tambien pequeno. 



Problema 4.3 



iPorque la ecuacion diferencial que satisfacen las funciones WKB asoci- 
adas, difiere de la ecuacion de Schrodingcr que es satisfecba por las funciones 
WKB; en la inclusion de la funcion s{x), si las funciones WKB y las asoci- 
adas WKB tienen la misma forma? 



Justificacion: 



Recordemos que en el proceso de obtencion de las soluciones WKB, nos 

encontramos con una ecuacion diferencial tipo Riccati; para la cual pro- 
pusimos una solucion en forma de serie de potencias de —i^, dicha serie 
es q{x) = J2'^=oi~^0"''ln{x)- Pero recordemos tambien que esta la aproxi- 
mamos solo hasta segundo orden, por lo que nuestras funciones if)^ (x) sat- 
isfacen la ecuacion dc Schrddinger solo aproximadamentc. Por otra parte se 
proponen las funciones WKB asociadas W±, como funciones que tienen la 
misma forma que las funciones ip'^; ypara obtener la ecuacion diferencial que 
estas satisfacen, simplemente las derivamos; y en consecuencia esta ecuacion 
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diferencial es satisfecha exactamente por ellas, y como vemos se introduce 
de manera natural la funcion s(x); y hace su aparicion para indicarnos que 
tanto se "desvian" las funciones ip^ de la solucion exacts a la ecuacion de 
Schrddinger 1 - dimensional. 
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5. EL OSCILADOR ARMONICO 
Solucion de la ecuacion de Schrodinger 

El oscilador armonico es quiza el modelo mas usado en la Fisica, y su 
utilidad va desde los campos de la Fisica clasica hasta la Electrodinamica 
cuantica. 

De la mecanica clasica sabemos que muchos potenciales complicados, pueden 
ser aproximados en la vecindad de sus puntos de equilibrio por un oscilador 
armonico 

V{x)^^V"{a){x-af (1) 

Esto en el caso unidimensional. Para este caso tenemos que la funcion 
hamiltoniana clasica, de una particula con masa m, oscilando con frecuencia 
Lo, toma la siguiente forma : 

H = ^ + -mu?x^ (2) 
2m 2 ^ ^ 

y el correspondiente hamiltoniano cuantico en el espacio de con&guraciones 
es : 

Dado que el potencial es independiente del tiempo, lo que determina las 
eigenfunciones y sus correspondientes eigenvalores En, es la ecuacion de 
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Schrodinger independiente del tiempo : 

H^n = En^n (5) 

Considerando el hamiltoniano para el oscilador armonico , se tiene que 
la ecuacion de Schrodinger correspondiente es : 



Hemos suprimido los submdices de E y por comodidad. Definiremos 
ahora las siguientes cantidades: 



Con estas deGniciones, nuestra ecuacion de Schrodinger es: 



+ [r - A^x^]* = (9) 



A esta ultima ecuacion se le conoce como "ecuacion diferencial de Weber" 
Haremos enseguida la transformacion: 

y = Ax^ (10) 

En general, en un cambio de variable, suponiendo que hacemos el cambio 
de la variable x a la variable y , se tiene que los operadores diferenciales 
toman la forma siguiente: 

d_^dyd_ 
dx dx dy 
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dx"^ dx dx dy dx"^ dy dx dy"^ 
Aplicando esto a la transformacion propuesta obtenemos la siguiente 
ecuacion diferencial en la variable y : 

o bien, definiendo : 

e p E 

= TZV = 7i = 1— (14) 

Obtenemos entonces la ecuacion: 



Id^ 1 ,^ ^ 



Pasaremos enseguida a resolver esta ecuacion, haciendo primeramente 
el analisis asintotico en el Umite y ^ oo , para hacer esto reescribimos la 
ecuacion anterior como sigue : 

I d-^ , K 1,^ ^ 
+ 77--7- + [77- -7)^ = (16) 



2y dy 2y i' 

Observamos que , en el Umite y — > 00, esta ecuacion se comporta as/; 

d^^^-oo 1 



dy^ 4 
Esta ecuacion tiene como solucion 



7*00 = (17) 



*oo(y) =^exp| + Bexp^ (18) 
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Desechamos exp | dado que esta diverge en el Umite y ^ oo, nos quedamos 
entonces con la exponencial decreciente. Podemos sugerir entonces que ^ 
tiene la siguiente forma: 

^'(y) = exp (19) 
Sustituyendo esto en la ecuacion diferencial para y ( ec. 15) se tiene: 



Lo que hemos obtenido es la ecuacion hipergeometrica confiuente f\ : 

La solucion general a esta ecuacion es : 

y{z)=A iFi{a;c,z) +B z^^^ iFi(a - c + 1; 2 - c, z) (22) 

Con la funcion hipergeometrica confiuente : 

,F,ia;c,z) = f:^-^ (23) 

Comparando entonces nuestra ecuacion , con la ecuacion hipergeometrica 
confiuente , se observa que la solucion general a nuestra ecuacion es : 

1 X 3 

ij{y) = A ^Fi{a;-,y) + B y\ ^Fi{a + -;-,y) (24) 



donde 



-(f-i) (25) 



^tambien conocida como la ecuacion diferencial de Kummer 
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Si dejamos estas soluciones asi como estan, entonces la condicion de nor- 
malizacion de la funcion de onda no se cumple, pues del comportamiento 
asintotico de la funcion hipergeometrica conRuente ^ se sigue que ( con- 
siderando linicamente el comportamiento exponencial, dado que es el dom- 
inante) : 

— y y 1 

^{y) = e~'4){y) const, e^y'^"^ (26) 



Esto nos lleva a una divergencia en la integral de normalizacion, lo cual 
es fisicamente inaceptable. Lo que se hace entonces, es imponer la condicion 
de terminacion de la serie f\ , esto es, la serie se corta y surge entonces un 
polinomio de grado n. 

Observamos entonces que, el hecho de pedir que la integral de normalizacion 
sea finita (para tener significado fi'sico en terminos de probabilidades), nos 
lleva al truncamiento de la serie, y esto a su vez es lo que da lugar a la 
cuantizacion de la energia. 
Consideremos enseguida los dos posibles casos : 
1) a = -n yB = 

- - - = n (27) 
2 4 ^ ^ 



® El comportamiento asintotico para j a; |^ cxd es: 



^La condicion de truncamiento de la serie para la funcion hipergeometrica confluente 
iFi{a; c, z) es a = — n, con n un entero no negativo ( esto es, incluye el cero ). 
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Con las eigenfunciones: 

—Xx 



=L>„exp— ^ iFi(-n;-,Aar2) (28) 



y la energia: 



= nuj{2n + ^) (29) 



2) a+\ = -n yA = Q 



i-i = "4 (30) 
Teniendo las eigenfunciones siguientes: 

-Xx^ 3 
*„(x) =D„exp^— X iFi(-n;-,Aa;2) (31) 

y para la energia se tiene: 

En = nu[{2n + 1) + ^] (32) 

Los polinomios dados anteriormente por las funciones hipergeometricas, 
son conocidos como polinomios de Hermite, y estan deGnidos en terminos 
de la funcion hipergeometrica como sigue : 

H2niv) = (-1)"^ iFi(-n;i r/2) (33) 
n! 2 

H2n-i{v) = i-ir^^^^^ r, iFi(-n; In') (34) 

Podemos finalmente combinar los resultados obtenidos ( pues unos nos 
dan los valores pares y otros los impares ) en una sola expresion para los 
eigenvalores y las eigenfunciones , obteniendose : 

-Xx"^ 

*„(x) = Dn exp -^Hn{\fXx) (35) 
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En = {n + -)hu n = 0,l,2... (36) 

El espectro dc cncrgi'a del oscilador armonico es equidistante, esto es, 
existe la misma diferencia Hlo entre cualesquiera dos estados . Otra obser- 
vacion que podemos hacer, es acerca del mmimo valor de energia que toma 
el oscilador; lo sorprendente es que este es distinto de cero; esto es un resul- 
tado puramente mecanico cuantico, a este valor se le conoce como energi'a 
de punto cero y el hecho de que sea distinta de cero , es una caracteristica 
de todos los potenciales ligantes (aquellos que confinan a las particulas) . 



La constante de normalizacion puede ser calculada, y tiene el valor: 

Con lo cual obtenemos Gnalmente las eigenfunciones del oscilador armonico 
unidimensional, normalizadas : 

/A 1 1 — At^ 

-^)^ exp^ Hr,{V\x) (38) 

Operadores de creacion (a)) y aniquilacion (a) 

Existe otra forma de tratar el oscilador armonico de una forma distinta a 
la convencional de resolver la ecuacion de Schrddinger, esta otra manera es 
Uamada el metodo algebraico o metodo de operadores, este es un poderoso 
metodo el cual es aplicado tambien en otra clase de problemas mecanico 
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cuanticos. 

DeRniremos dos operadores no hermitianos a y : 

Imuj , iv . , . 

V 2ft mu 

a} = J — (x-—) (40 

Estos operadores son conocidos como operador de aniquilacion y op- 
erador de creacion, respectivamente (las razones para estos nombres se 
veran despues ). 

Vamos a calcular ahora el conmutador de estos dos operadores: 

i muj , ip ip , 1 

[a, at] = + — >^ - — ] = 2^(-^[^,P] +^b,^]) = 1 (41) 

Donde hemos usado el conmutador: 

[x,p]=in (42) 

Esto es, tenemos que los operadores de creacion y aniquilacion satisfacen 
la relacion de conmutacion : 

[a, at] = 1 (43) 
Vamos a deGnir tamhien el Uamado operador de numero N como: 

N = a^a (44) 

Este operador es hermitiano como lo podemos demostrar facilmente 

usando {AB)^ = B^A^ : 

iVt = (ata)t = at(at)t =a^a = N (45) 
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Ahora bien, considerando que : 

observamos que el hamiltoniano est a dado de una man era simple en 
terminos del operador de numero : 

H = noo{N + (47) 

El operador de numero recibe su nombre debido a que sus eigenvalores 
son justo el mdice de la funcion de onda sobre la que opera, esto es: 

N I n >= n\n> (48) 

Donde hemos usado la notacion: 

I *n > = \n> (49) 

Aplicando este hecho a (47) tenemos : 

H \n>=fiu}{n + ^)\n> (50) 

Pero sabemos de la ecuacion de Schrodinger que H \ n >= E \ n > de 
lo cual se sigue que los valores de la energia estan dador por : 

En = huj{n + ^) (51) 

El cual es identico (como debia ser ) con el resultado (36). 
Vamos enseguida a mostrar el por que de los nombres que se le dan a los 
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operadores a y . Para hacer esto comenzaremos calculando dos conmu- 
tadores: 

[N, a] = [a' a, a] = a' [a, a] + [a^, a]a = —a (52) 
Lo anterior se sigue de[a,a] =07 (AS) .Similarmente calculamos: 

[N,a^ = [a^a,a^] = a^[a,a^] + [a^,a^]a = (53) 

Con estos dos conmutadores podemos escribir: 

n > 

= {of + a^N) I n > (54) 
= a\l + n) \ n >= (n + l)a^ | n > 

Con un procedimiento similar se obtiene tambien: 

N{a I n >) = {[N, a] + aN) \ n >= {n - l)a \ n > (55) 

La expresion (54) implica que el ket \ n > es un eigenket del operador 
de numero , donde el eigenvalor se ha incrementado por uno, esto es , se 
ha creado un cuanto de energia al actuar a} sobre el ket, de ahi su nombre 
de operador de creacion. Comentarios siguiendo la misma Hnea de 
razonamiento son validos para el operador a, lo cual le da el nombre de 
operador de aniquilacion ( un cuanto de energia es disminmdo al actuar 
este operador ). 

La ecuacion (54) tambien implica que el ket \ n > y el ket \ n + 1 > son 
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proporcionales, podemos escribir esta relacion as/; 

I n >= c \ n + 1 > (56) 

Donde c es una constante que hay que determinar. Considerando ademas 
que : 

(a^ I n >)^ =< n\a = c*<n + l\ (57) 
Podemos entonces realizar el siguiente calculo: 

< n I a(a^ \ n >) = c* < n + 1 \ {c \ n + 1 >) (58) 

< n I aa^ \ n >= c*c <n + l\ n + l> (59) 
< n I aa^ \ n >=| c |^ (60) 
Pero de Ja relacion de conmutacion para los operadores a y a} : 

[a, a^] = aa^ — a^a = aa^ — N = 1 (61) 

Esto es : 

aa^ = N + 1 (62) 

Sustituyendo en (60); 

<n \ N + l\ n >=< n \ n > + < n \ N \ n >= n + 1 =\ c f (63) 
Pidiendo que c sea real y positiva ( por convencion ), obtenemos su valor: 

c = Vn + 1 (64) 
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Con lo cual se tiene la relacion: 

\n>= Vn + l | n + 1 > (65) 

Siguiendo el mismo camino se Uega a una relacion para el operador de 
aniquilacion : 

a I n >= -y/n | n — 1 > (66) 

Vamos ahora a mostrar que los valores de n deben ser enteros no neg- 
ativos. Para esto, acudiremos al requerimiento de positividad de la norma, 
aplicado cn especial al vector de estado a \ n>. Este requerimiento nos dice 
que el producto interior de este vector con su adjunto({a \ n >)^ =< n \ a^) 
debe ser mayor o igual que cero : 

{<n\a^)-{a\n>)>0 (67) 

Pero lo anterior no es mas que : 

<n\a^a\n >=< n \ N \ n>=n>0 (68) 

Por lo tanto n nunca puede ser negativo. Y tiene que ser entero pues 
si no lo fuera al aplicar en repetidas ocasiones el operador de aniquilacion 
nos Uevaria a valores negativos de n, lo cual esta en contraposicion con lo 
anterior. 

Es posible expresar el estado n {\ n >) en terminos del estado base (| >) 
usando el operador de creacion, veamos como hacerlo: 
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1 >= aT I > 



2>^[^]|1>^[ 



(at) 



t^2 



■]|o> 



at fat)3 



(69) 
(70) 

(71) 



fat)" 

n >= [ V=-] I > (72) 



Podemos tambien aplicar este metodo para encontar las eigenfunciones 
en el espacio de conBguraciones. Para hacer esto, partiremos del estado 
base: 

a I >= (73) 

En la representacion x tenemos: 



a*oW = ^^(.+ £;)*oW=0 (74) 

Recordando la forma que toma el operador momento en la representacion x, 
podemos Uegar a una ecuacion diferencial para la funcion de onda del estado 
base; introduciremos tambien la deBnicion siguiente xq = J , con esto : 



{x + xl^)^Q = Q (75) 

Esta ecuacion se puede resolver facilmente, resolviendola y normalizandola 
( su integral de — oo a oo debe ser la unidad ), Uegamos a la funcion de onda 
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del estado base: 



M^) = ( 



1 




(76) 



Las demas eigenfunciones, esto es, las eigenfunciones para los estados ex- 
citados del oscilador armonico , se pueden obtener usando el operador de 
creacion, el procedimiento es el siguiente: 



Siguiendo con este procedimiento, por induccion se puede mostrar que: 



Evolucion temporal del oscilador 

En esta seccion vamos a ilustrar con el oscilador armonico una manera en la 
cual se trabaja con la representacion de Heisenberg, esto es, dejaremos que 
los estados esten Rjos en el tiempo y haremos evolucionar a los operadores 
en el. Veremos a los operadores como funciones del tiempo, especiGcamente, 
encontraremos como es que evolucionan los operadores posicion , momento, 
a y en el tiempo, para el caso del oscilador armonico . Las ecuaciones de 
movimiento de Heisenberg para pyx son : 




(77) 



(78) 




(79) 



dp d 
-dt = 



(80) 
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De aqm se sigue que las ecuaciones de movimiento para x y p en el caso 
del oscilador armonico son: 

^ = -mco'^x (82) 
dt ^ ' 

f = I (83) 
dt m ^ ^ 

Se tiene un par de ecuaciones acopladas , estas son equivalentes a un 

par de ecuaciones para los operadores de creacion y aniquilacion, salvo que 

estas dos ultimas no estan acopladas, veamoslas explicitamente : 

da jmoj d ,^ ip . ,^ ^, 



da jmco ,dx i dp, , , 

" y^^'dt'^^W ^ ^ 

Sustituyendo (82) y (83) en (85) : 

da frmv p . . 

Similarmente se puede obtener una ecuacion diferencial para el operador de 
creacion, la cual no esta acoplada: 

da^ + 

— = luja^ (87) 
dt ^ ' 

Las ecuaciones diferenciales que hemos encontrado para la evolucion tem- 
poral de los operadores de creacion y aniquilacion , pueden ser integradas 
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inmediatamente, dandonos la evolucion expHcita de estos operadores en el 
tiempo : 

a{t) = a{0)e-''^^ (88) 
a^{t) = at(0)e^'"* (89) 

Podemos observar de estos resultados y de las ecuaciones (44) y (47) , 
que tanto el hamiltoniano como el operador de numero , no dependen del 

ticmpo, tal y como podnamos csperar. 

Con los dos resultados anteriores , podemos encontrar los operadores de 
posicion y momento como funcion del tiempo, pues ellos estan dados en 
terminos de los operadores de creacion y aniquilacion: 



P = ^y^(at-a) (91) 

Sustituyendo los operadores de creacion y aniquilacion se obtiene: 

x(t) = x(0) cos Lot + --^ sin ujt (92) 

ma; 



p{t) = — mujx (0) sin Lot + p(0) cos tot (93) 

La evolucion temporal de los operadores de posicion y momento es la 
misma que las ecuaciones clasicas de movimiento. 

Hemos Gnalizado esta seccion, mostrando la forma expHcita en que evolu- 
cionan cuatro operadores en el caso del oscilador armonico , reAejando de 
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esta manera una forma de trabajar en la poco mencionada representacion 
de Heisenberg. 

El oscilador armonico tridimensional 

Al iniciar nuestro estudio cuantico del oscilador armonico , haciamos co- 
mentarios acerca del por que la importacia del oscilador armonico . Si 
hicieramos un analogo tridimensional, consider ariamos entonces un desar- 
roUo de Taylor en tres variable^ reteniendo terminos solo hasta segundo 
orden, lo que tenemos es una forma cuadratica (en el caso mas general), el 
problema de resolver para esta aproximacion no es sencillo, es decir, para el 
caso : 



Afortunadamente hay varios sistemas que se ajustan bien a la simetria 
esferica, esto es, para el caso: 



Esto ultimo, equivale a decir que las parciales segundas ( no cruzadas 

) toman todas el mismo valor ( en el caso anterior representado por K), y 
®Es posible expresar el desarroUo de Taylor como un operador exponencial como sigue: 



V{x, y, z) = ax^ + by^ + cz'^ + dxy + exz + fyz 



(94) 




(95) 




Esto es el desarroUo de Taylor en tres variables alrededor de To. 
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podnamos agregar que esta es una huena aproximacion en el caso en que los 
valores de las parciales cruzadas sean pequenas comparadas con las parciales 
segundas no cruzadas. 

Cuando se satisfacen los requerimientos anteriores, y tenemos un potencial 
como el dado por (95), entonces tenemos el denominado oscilador armonico 
tridimensional esfericamente simetrico. 
El hamiltoniano para este caso es de la forma: 

H = V + r (96) 

Donde el laplaciano esta dado en coordenadas esfericas yresla variable 
esferica convencional. 

Tenemos entonces que el potencial es independiente del tiempo, por tanto 
la energia se va a conservar; ademas dada la simetria esferica , el momento 
angular se conservara tambien, se tienen por tanto dos cantidades conser- 
vadas, y puesto que a cada cantidad conservada le corresponde un numero 
cuantico, podemos adelantar que nuestras funciones de onda dependeran de 
dos numeros cuanticos (aunque en este caso, como veremos surge otro ). 
Esto es , necesitamos solucionar la ecuacion : 

H^nl = Enl^nl (97) 

El laplaciano en coordenadas esfericas es : 
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2 d"^ 2 d L2 
v' = ^ + -^-T23 (98) 



Qj.2 J. Qj. ^2^2 

Esto se sigue del hecho que : 



L' = -^'[-^^(sin^^) + -^#^] (99) 
Las eigenfunciones de son los armonicos esfericos, se tiene: 



PYimie, if) = -hH{i + i)Yim{e, ^) (loo) 

Podemos observar que el hecho de que los armonicos esfericos Ueven el 
niimero cuantico m , produce la intromision del mismo en la funcion de 
onda, es decir tendremos '^nlm- 

Para separar la ecuacion diferencial se propone la sustitucion: 

*nim(r, 9, ^) = ^Il^Yimie, ^) (101) 

Esto al sustituirlo en la ecuacion de Schrddinger, nos va a separar la 
parte espacial de la parte angular; la parte angular , son las eigenfunciones 
del operador momento angular ( al cuadrado ), en la parte espacial llegamos 
a la ecuacion : 



<, + (^-^r^-^)fl„,M = o (102) 
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Usando las definiciones (7) y (8) , la ecuacion anterior toma exactamente 
la misma forma que (9), excepto por el termino del momenta angular, este 
termino usualmente es Uamado la barrera de momento angular. 



+ - - = (103) 

Para resolver esta ecuacion , partiremos del analisis asintotico de la 
misma. Si consideramos primero el Um r oo, observamos que el termino 
del momento angular es despreciable, de manera que el comportamiento 
asintotico en este limite es identico a (9) con lo cual obtenemos: 

Rni{r) ~ exp — ^ — en lim r ^ oo (104) 

Si observamos ahora el comportamiento cerca de cero, vemos que el 
comportamiento dominante esta dado por el termino de momento angular, 
es decir, la ecuacion diferencial (102) se convierte en este limite en : 

Kl - ^^T^^ni = (105) 

Esta es una ecuacion diferencial tipo Euler , solucionandola encon- 
tramos dos soluciones independientes: 

Rnl{r) ~ o r"' en lim r ^ (106) 



®Una ecuacion tipo Euler es : 

a;"j/'-"' (x) + x"^^ y'^"'^^'^ (x) + ■ • • + xy'{x) + y(x) = 
La cual tiene soluciones del tipo x", se sustituye y se encuentra un polinomio para a. 
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Lo anterior nos Ueva a proponer la sustitucion : 

Rnl{r) = r'+i exp ^</>(r) (107) 
Podriamos hacer tambien la sustitucion: 

Rni (r) = r-' exp -^v{r) (108) 

Sin embargo esto Ueva a las mismas soluciones que (107) ( mostrar esto 
es un buen ejercicio). Sustituyendo (107) en (103) , se obtiene una ecuacion 
diferencial para (f) : 

41' + 2(^ - \t)^ - \\{2l + 3) - k^U = (109) 
r 

Haciendo ahora la sustitucion de la variable w = Xr^, obtenemos: 

+ {1 + 1- -i^i^ + l)- = (110) 

Donde hemos introducido k- = ^ = Tenemos nuevamente una ecuacion 
diferencial tipo hipergeometrica confiuente la cual tiene por soluciones ( 
vease {21) y (22)): 

ct^{r) = A ,F,[^(l+^-n);l+^,Xr']+B r-(2'+^ ,i.,[i(_/+i_^); , Ar^] 

(111) 

La segunda solucion particular no puede ser normalizada , pues diverge 
fuertemente en cero, de manera que tomamos B = y se tiene : 

<t>{r) = A ,F^[l{l + \-n)-l + \,\T^] (112) 
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Empleando los mismos argumentos que para el oscilador unidimensional, es 
deck, pedimos que las soluciones sean regulares en el inRnito, nos Ueva a 
la condicion de truncamiento de la serie, lo cual nos Ueva nuevamente a la 
cuantizacion de la energia ; imponiendo la condicion de truncamiento: 

lil + l-^) = -n (113) 

Esto es , poniendo expUcitamente k, obtenemos el espectro de energia : 

Eni = nco{2n + 1 + 1) (114) 

Podemos observar que se tiene una energia de punto cero igual a |^ 
para el oscilador armonico tridimensional esfericamente simetrico. 
Las eigenfunciones son (no normalizadas ): 

^nlm{r,e,ip)=r^e^ iF,{-n;l + ^, Xr^) Yim{e,ip) (115) 



113 



Problemas 

Problema 5.1 

Encuentre los eigenvalores y eigenfunciones del oscilador armonico 
en el espacio de momentos 

El hamiltoniano mecanico cuantico para el oscilador armonico esta dado 
por: 

H = — mu X 

2m 2 

Ahora bien , en el espacio de momentos, los operadores x y p toman la 
siguiente forma : 

P 

X — ^ in,Tr 
op 

Por tanto el hamiltoniano mecanico cuantico para el oscilador armonico en 
el espacio de los momentos es : 

TT P ^ 2f-2 " 

H = mco n -—^ 

2m 2 ap'^ 

Tenemos entonces que resolver el problema de eigenvalores ( esto es, encon- 
trar las eigenfunciones y los eigenvalores) dado por (5) , lo cual nos lleva, 
con el hamiltoniano anterior, a la siguiente ecuacion diferencial : 

cip mn m'^Ti u'^ 
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Se puede observar que la ecuacion diferencial obtenida, es identica, salvo 
constantes, con la ecuacion diferencial que obtuvimos en el espacio de con- 
figuraciones ( ec. (6) ). Solo que para ejempUGcar otra forma de resolverla, 
no seguiremos exactamente el niismo camino que se siguio para obtener la 
solucion de aquella. 

DeRniremos dos parametros, de manera similar a como hicimos en (7) y (8); 

fc2 ^ A = ^ (117) 

Con estas definiciones, arribamos exactamente a la ecuacion diferencial (9), 
se sigue por tanto que la solucion buscada ( despues de realizar el analisis 
asintotico ) es de la forma: 

^{y) = e--^my) (118) 

Dondey estadadapory = Xp^ yX definidaen (117). Sustituiremos entonces 
(118) en (116), solo que regresando (118) a la variable p. Haciendo esta 
sustitucion se obtiene una ecuacion diferencial para (f) : 

^-2Ap^ + (*^-AWp) = (119) 

Haremos Gnalmente el cambio de variable u = VXp, la ecuacion anterior se 
transforma en la ecuacion diferencial de Hermite: 

d^6(u) „ d6(u) , , , , , 

Con n un entero no negativo, y donde hemos hecho: 

^ - 1 = 2n 
A 
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De aqm y de las definiciones dadas en (117) se sigue que los eigenvalores de 
la energia estan dados por: 



Y las soluciones a (120) estan dadas por los polinomios de Hermite, de 
manera que (f){u) = Hn{u), con lo cual se tiene que las eigenfunciones (no 

normalizadas) estan dadas por: 



Problema 5.2 

Demuestre que los polinomios de Hermite pueden ser expresados 
con la siguiente representacion integral: 



La representacion anterior de los polinomios de Hermite es una poco 
usual, pcro que es muy util en ciertos casos. Lo que vamos a hacer para 
demostrar la igualdad, es desarroUar la integral que se presenta y mostrar 
que lo obtenido es identico con la representacion en serie de los polinomios 
de Hermite, la cual esta dada por: 



En = nu{n + -) 



^{p) = Ae-^P'^HniVXp) 




(121) 



n 




(122) 
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Donde el simbolo [c] signiEca, el mayor entero menor o igual que c. 
Lo primero que haremos es desarroUar el binomio que esta dentro de la 
integral usando el teorema del binomio: 

n 1 

{x + y)" = E . '-^^x^y"' 

(n — mjimi 

Usando esto, el binomio dentro de la integral, tiene el desarroUo siguiente: 

{x + iyr = E 7 (123) 

in — m)!m! 

m=0 ^ ' 

Sustituyendo esto en la integral: 

De la forma del integrando podemos ver que la integral es distinta de 

cero cuando m es par, pues de lo contrario el integrando seria impar y la 
integral se anula. Haremos por tanto el cambio m = 2k; con este cambio se 
tiene: 

[-] 

T^tl^^^Wm^'^"-'' 2f ,-e-»'* (125) 

Con el cambio de variable u = la integral se convierte en una funcion 

gamma, haciendo el cambio de variable: 
f-1 

on 12 1 I poo , 

V^^^{n-2k)\{2k)\ Jo ^ ' 

La integral es precisamente r(A; + ^) , la cual puede ser expresada en 
forma de factoriales ( para k entero, desde luego) : 
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Sustituyendo este valor en la sumatoria y usando el hecho de que i = (—1) 

[HI 

El cual es identico con (122), con lo cual se completa la demostracion. 
Problema 5.3 

Muestre que los eigenestados del oscilador armonico satisfacen la 
relacion de incertidumbre 

Debemos mostrax que para cualesquier eigenestado '^n se satisface: 

<{Ap)\Axf> >^ (128) 
Donde la notacion <> significa promedio. 

Vamos a calcular por separado < (Ap)^ > 7 < (Ax)^ >, donde cada una de 
estas expresiones es: 

< (Ap)^ >=< (p- < p >)^ >=< p'^-2p <p> + <p >^>=< p^ > - <p>^ 

< (Ax)^ >=< (x- < X >)^ >=< x^-2x < X > + < X >^>=< x^ > - < X >^ 

Vamos primeramente a mostrar que tanto el promedio de x, como el de 
p se anulan. Consideremos primeramente el promedio de x: 

/oo 
x[*„(x)]^dx 
-00 
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Esta integral se anula pues ['^n{x)]'^ es una funcion par, esto se puede 
ver considerando que la paridad esta dada por la parte polinomial ( pues la 
exponencial involucrada es una funcion par ). Los polinomios de Hermite 
tienen paridad de&nida, y se tienen solo dos casos, n es par o es impar. Si n es 
par se sigue de inmediato que lo es. Si n es impar entonces tenemos 

que Hn{—x) = (— l)"iJ„(x), e inmediatamente se ve que esta funcion al 
cuadrado es par tambien ( esto es, cualquier polinomio par o impar, elevado 
al cuadrado es par ). Hemos mostrado que es una funcion par para 

n cualquiera, por tanto al multiplicarla por x se vuelve impar, de manera 
que la integral se anula. Tenemos entonces como resultado: 

<x>=0 (129) 

Los mismos argumentos son validos para el promedio de p, si calculamos 
este en el espacio de momentos con las funciones encontradas en el problema 
1, pues la forma funcional es la misma. De manera que: 

<p>=0 (130) 

Calculemos ahora el promedio de x^. Para hacer esto usaremos el teo- 
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rema del virial |^. Observemos primeramente que : 

< V >= ]^muP' <x^ > . 

De manera que es posible relacionar el promedio de con el promedio del 
potencial ( y poder usar el teorema del virial). 



< >= <V > . (131) 



Necesitamos considerar tambien el promedio de la energia : 

< H >=< T > + <V > 
Usando el teorema del virial ( para n = 2 ) se obtiene: 

< H >=2 <V > (132) 

Con lo cual se obtiene: 

< >= ^ = ^ , (133) 



<x'^ >= —{n + l) (134) 



El teorema del virial en mecanica cuantica nos dice que: 

2 < T >=< r ■ v1^(r) > 
Para un potencial de la forma V = \x" se satisface: 

2<T>=n<V > 
Donde T representa la energia cinetica y V la energia potencial. 
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De forma similar calculamos el promedio de p^, expUcitamente: 

< >= 2m, < >= 2m < T >= m < H >= mfiojin + -) (135) 
2m 2 

Con (133) y (135) se tiene: 

< {ApfiAxf >= (n + ^fh^ (136) 

De este resultado inmediatamente se puede ver que los eigenestados sat- 
isfacen la relacion de incertidumbre, con el mmimo valor precisamente para 
el estado base ( n = ). 

Problema 5.4 

Obtenganse los elementos de matriz de los operadores a, a\ x y p 

Encontraremos primero los elementos de matriz de los operadores de 
creacion y aniquilacion, ya que estos nos ayudaran a encontrar los elementos 
de matriz para los otros dos operadores. 
Usaxemos las relaciones (65) y (66), con las cuales se tiene: 

< m \ a \ n >= \/n < m \ n — 1 >= \/n^m,n-i (1^'^) 

Similarmente para el operador de creacion se tiene: 

< m \ \ n >= Vri^^ < m \ n + 1 >= Vn + lSm,n+i ■ (1^8) 

Pasaremos enseguida a calcular los elementos de matriz del operador de 
posicion. Para hacer esto, expresaremos el operador de posicion en terminos 
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de los operadores de creacion y aniquilacion. Usando las de&niciones (39) y 
(40), se comprueba inmediatamente que el operador de posicion esta dado 
por: 

' (a + at) . (139) 



V 2muj 

Usando esto, los elementos de matriz del operador x pueden ser inmedi 
atamente calculados: 



< m \ X \ n > = < m \ \ -r-—(a + a^) I n > 

V 2mu 



n 



■[VnSm,n-l + Vn+ldm,n+l] (140) 



V 2mu 

Siguiendo el mismo procedimiento podemos calcular los elementos de matriz 
del operador momenta, considerando que p esta dado en terminos de los 
operadores de creacion y aniquilacion, de la forma que a continuacion se 
muestra: 

P = i\ 

Usando esto se tiene: 



'^(a^-a) (141) 



' mhio 



<m\p\n>=i^ —^[Vn + lSm,n+i - VnSm,n-i] (142) 

Hemos encontrado entonces , los elementos de matriz de los cuatro oper- 
adores y se puede observar la sencillez con la cual son calculados los mismos 
para los operadores de posicion y momento con la ayuda de los operadores 
de creacion y aniquilacion. Finalmente podemos hacer la observacion acerca 
de la no diagonalidad de los elementos de matriz que hemos encontrado, lo 
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cual era de esperarse por el hecho de que la representacion que estamos 
usando es la del operador de numero, y ninguno de los cuatro operadores 
conmuta con el. 

Problema 5.5 

Encuentrense los valores esperados de x"^ y para el oscilador 
armonico unidimensional y lisense estos para encontrar los valores 
esperados de la energia cinetica y la energia potencial. Comparese 
este ultimo resultado con el teorema del virial 

Se encontrara primeramente el valor esperado de x^. Para hacerlo re- 
curriremos a la expresion (139), de la cual se sigue: 

= {a^ + (at)2 + a^a + aa^) (143) 

Recuerdese que los operadores de creacion y aniquilacion no conmutan entre 
SI. Con lo anterior podemos calcular el valor esperado de : 

<x^> = <n\x^ \ n> 

n 



n{n - l)dn,n-2 + J {n + l){n + 2)Sn,n+2 



2mL0 

+ n + (ra + 1) ] (144) 

Esto es, el valor esperado de x^ esta dado por: 

< x^ >=< n\ x'^ \n>= -^—(2n + 1) (145) 

2moj 

Para calcular el valor esperado de solo necesitamos expresar este operador 
en terminos de los operadores de creacion y aniquilacion, lo cual lo podemos 
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realizar a partir de (141), obteniendose: 



f = -^(a^ + (at)2 - aat - a^a) (146) 



Para el valor esperado del cuadrado del momento se obtiene: 

<p^ >=<n\p^ \n>= ——{2n + l) (147) 



Con este ultimo resultado podemos encontrar el valor esperado de la energia 
cinetica : 

< T >=< i— >= — < p2 > (148) 
Usando el (147), se obtiene el valor esperado de la energia cinetica: 

< f >=< n I f I n >= ^(2n + 1) . (149) 
Nos falta obtener ahora el valor esperado de la energia potencial: 

< V >=< ^mJ^x^ >= ^muj^ <x^ > . (150) 

Con el resultado de (145), inmediatamente se sigue el valor esperado para 
la energia potencial: 

< y >=< n I T> I n >= ^(2n + 1) . (151) 

Observamos que los valores esperados de la energia cinetica y la energia 
potencial coinciden para toda n, esto es, para cualquier estado de energia, lo 
cual esta en correspondencia con el teorema del virial el cual establece que 
para un potencial cuadratico, como el del oscilador armonico , los valores 
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esperados de la energi'a cinetica y la energi'a potencial dehen coincidir, y mas 
aun estos deben ser iguales a un medio del valor esperado de la energia total 
del sistema, lo cual efectivamente se satisface. 

Problema 5.6 

Una particula cargada (carga q) se mueve en la direccion z, en la 
presencia de un campo magnetico uniforme en la misma direccion ( 
B = Bk) . Comparando el hamiltoniano para este sistema con el del 
oscilador armonico unidimensional, muestre que los eigenvalores 
de la energia pueden ser inmediatamente escritos: 



donde Uk es el eigenvalor continuo del operador Pz Y n un entero 
no negativo. 

El hamiltoniano para una particula de carga electrica q, la cual se mueve 
en presencia de un campo electromagnetico esta dado por: 



donde A es el potencial vectorial que genera el campo magnetico y (j) es el 
potencial escalar que genera el campo electrico. 

En nuestro problema no tenemos campos electricos presentes, de manera que 
el potencial escalar cf) es igual a cero. Nuestro hamiltoniano toma entonces 





(152) 
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la forma: 



H = ^--^{p.A + A-p) + J-i^fA^ . (153) 

Esto dado que el potencial vectorial y el momento no conmutan, pues el 
potencial vectorial es funcion de las coordenadas. 

Dado que la particula se desplaza solo en la direccion z, solo tenemos mo- 
mento lineal asociado a esta coordenada, esto es : 

P={o,o,pz) (154) 

Y por otro lado, el potencial vectorial A que genera el campo magnetico en 
la direccion z es: 

A = {-By, 0,0). (155) 

Esto ultimo se puede comprobar sabiendo que B = \/ x A. 

Bajo estas circunstancias, el segundo termino del hamiltoniano en (153) se 

anula, de manera que el hamiltoniano a considerar es: 

Observamos que el hamiltoniano es la suma de una parte de particula libre 
y otra de oscilador armdnico , identihcando lo = podemos escribir la 
energia asociada a cada contribucion inmediatamente, con lo cual se obtiene: 

n'^k^ \qB\h. 1, 

Obteniendo as/ el resultado buscado. 
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6. EL ATOMO DE HIDROGENO 



Se estudia el atomo de hidrogeno resolviendo la ecuacion de Schrddinger 
independiente del tiempo con un potencial debido a dos particulas car- 
gadas como lo son el electron y el proton, con el Laplaciano en coordenadas 
esfericas, mediante separacion de variables, dando una interpretacion fisica 
de la funcion de onda como una solucion de la ecuacion de Schrddinger 
para el atomo dc hidrogeno, ademas de las interpretaciones dc los mimeros 
cuanticos y de las densidades de probabilidad. 



INTRODUCCION A LA MECANIC A CUANTIC A 

Como nuestro interes es el de describir el atomo de hidrogeno, el cual esta a 
una escala muy pequeha, se hard mediante el uso de la mecanica cuantica, la 
cual trata las relaciones entre magnitudes observables, pero el principio de 
incertidumbre altera radicalmente la dehnicion de "magnitud observable" en 
el campo atomico. De acuerdo con el principio de incertidumbre, la posicion 
y el momento de una particula no se pueden medir simultaneamente con 
precision. Las cantidades cuyas relaciones busca la mecanica cuantica son 
probabilidades. En vez de ahrmar, por ejemplo, que el radio de la orbita 
del electron en un estado fundamental del atomo de hidrogeno es siempre 
exactamente 5.3 x 10~^^ m, la mecanica cuantica ahrma que este es el radio 
mas probable; si reahzamos un experimento adecuado, la mayor parte de las 
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pruebas daran un valor distinto, mas grande o mas pequeho, pero el valor 
mas probable sera aproximadamente 5.3 x 10"^^ m. 



Como ya se sabe, la cantidad con que esta relacionada la mecanica cuantica 
es la funcion de onda de una particula. Aunque ^ no tiene interpretacion 
fisica, el cuadrado de su valor absoluto | ^' P calculado para un punto 
y en un instante determinado es proporcional a la cantidad de encontrar 
experimentalmente a la particula ahi y en ese instante. El problema de la 
mecanica cuantica es determinar ^ para una particula cuando su libertad 
de movimiento esta limitada por la accion de fuerzas externas. 

Antes de considerar el calculo real de debemos establecer algunos 
requisitos que siempre se deben cumplir. En primer lugar, ya que | ^' p es 
proporcional a la probabilidad P de encontrar a la particula descrita por 
la integral de | * p sobre todo el espacio debe ser Rnita, ya que la particula 
esta en alguna parte. Si tenemos que 



la particula no existe y la integral evidentemente no puede ser oo y tener 



seguido para deGnirla, y asi la unica posibilidad dada es que su integral sea 
una cantidad Gnita para que ^ describa apropiadamente una particula real. 



ECUACION DE ONDA 





cierto signiHcado; \ ^ 



no puede ser negativa o compleja a causa del camino 
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Generalmente es conveniente tener | ^' p igual a la probabilidad P de encon- 
trar la particula descrita por ^, en lugar de ser simplemente proporcional a 
P. Para que | * P sea igual a P se tiene que cumplir la relacion 

/oo 
I * P = 1 (2) 
-oo 



ya que 



oo 



PdV = 1 (3) 



— oo 



es la aRrmacion matematica de que la particula existe en algun lugar en todo 
momcnto. Una funcion que obedezca a la ec. 2 se dice que esta normalizada. 
Ademas de ser normalizable, ^ debe tener un solo valor, ya que P debe tener 
un valor unico en un tiempo y en un lugar deter minados. Otra condicion 
que ^ debe obedecer es que ella y sus derivadas parciales ^) ^ scan 
continuas en cualquier lugar. 

La ecuacion de Schrddinger, que es la ecuacion fundamental de la mecanica 
cuantica, cn el mismo sentido que la segunda ley del movimiento es la 
ecuacion fundamental de la mecanica newtoniana, es una ecuacion de onda 
en la variable ^. Antes de abordar la ecuacion de Schrddinger repasemos la 
ecuacion de onda general 

que gobierna a una onda cuya cantidad variable es y que se propaga en la 
direccion de x con la velocidad v. En el caso de una onda en una cuerda 
tensa, y es el desplazamiento de la cuerda medido desde el eje x; en el caso 
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de una onda sonora, y es la diferencia de presion; en el caso de una onda 
luminosa y es la magnitud del campo electrico o la del campo magnetico. 

Las soluciones de la ecuacion de onda pueden ser de varios tipos, como 
consecuencia de la variedad de ondas que puede haber (un pulso unico en 
desplazamiento, un tren de ondas de ampUtud y longitud de onda constantes, 
un tren de ondas superpuestas de amplitudes y longitudes de onda identicas, 
un tren de ondas superpuestas de amplitudes y longitudes de onda diferentes, 
una onda estacionaria en una cuerda B.ja por ambos extremes, etc.). Todas 
las soluciones deben ser de la forma 



donde F es cualquier funcion que pueda ser diferenciada. Las soluciones 
F{t — x/v) representan ondas que viajan en el sentido +x, y las soluciones 
F{t + x/v) representan ondas que viajan en el sentido —x. Aqm nos interesa 
el equivalente ondulatorio de una particula "libre", es decir, una particula 
que no este bajo la inBuencia de ninguna fuerza y que, por lo tanto, viaja en 
una trayectoria recta a velocidad constante. Este equivalente corresponde a 
la solucion general de la ec. 4 para ondas armonicas no amortiguadas ( es 
decir, de amplitud constante A), monocromaticas (de frecuencia angular lo 
constante) en la direccion +x, 



En esta formula, y es una cantidad compleja, con parte real e imaginaria. 




(5) 




(6) 
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Como 

e"'^ = cosd - isenO (7) 

la ec. 6 se puede escrihir en la forma 

y{x, t) = Acosuj{t — x/v) — iAsenu{t — x/v) (8) 

Unicamente la parte real de la ec. 7 tiene significado en el caso de 
ondas en una cuerda en tension, donde y representa el desplazamiento de la 
cuerda con respecto a su posicion normal, en este caso la parte imaginaria 
se descarta porque no se puede aplicar. 

ECUACION DE SCHRODINGER 

En mecanica cuantica, la funcion de onda ^ corresponde a la variable de 
onda y del movimiento ondulatorio general. Sin embargo, a diferencia de 
y, no es una cantidad mensurable en si misma y puede, por tanto, ser com- 
pleja. Por esta razon supondremos que ^' esta especificada en la direccion 
X por 

*(a;,t) = y4e-^'^(*-^/") (9) 

Cuando se sustituye en esta formula u por 2itv y v por \v, obtenemos 

= Ae-^^^^'^*-^/^) (10) 

que es conveniente, ya que sabemos que v y \ estan en funcion de la energia 
total E y del momento p de la particula descrita por Ya que 

E = hiy = 27rni^ (11) 
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y 




2Trn 



p 



(12) 



tenemos 



{i/h){Et-px) 



(13) 



La ec. 13 es una descripcion matematica de la onda equivalente a una 
particula libre, de energia total E y momento p, que se mueve en la direccion 
y sentido +x, del mismo modo que la ec. 6 es la descripcion matematica de 
un desplazamiento de onda armonica que se mueve libremente a lo largo de 
una cuerda en tension. 

La expresion de la funcion de onda ^, dada por la ec. 13, es correcta 
solamente para particulas que se mueven libremente, pero estamos mas in- 
teresados en situaciones donde el movimiento de una particula esta sujeto 
a varias restricciones, como el caso de un electron ligado a un atomo por 
el campo electrico de su nucleo. Lo que debemos hacer ahora es obtener la 
ecuacion diferencial fundamental para ^, la que se puede resolver en una 
situacion especiRca. 

Comenzamos por la diferenciacion de la ec. 13 dos veces con respecto a 

X, 



y una vez respecto a t 




(15) 
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A velocidades pequenas comparadas con la de la luz, la energia total E 
de una particula es la suma de su energia cinetica /2m y de su energia 
potencial V, donde V es una funcion general de la posicion x y del tiempo 
t: 

Multiplicando ambos miembros de esta ecuacion por la funcion de onda 

= ^ + (17) 
2m 



De las ecs. 14 y 15 vemos que 



y 



E* = 4f (18) 



p=* = -«^|| (19) 
Sustituyendo estas expresiones de y en la ec. 17 obtenemos 

La ec. 20 es la ecuacion de Schrddinger dependiente del tiempo, donde la 
energia potencial V es una funcion de x,y,z,t. En tres dimensiones, la 
ecuacion de Schrddinger dependiente del tiempo es 

ih— - -— + ^ + (21) 

dt 2m I dx^ dy"^ dz^ J 

Una vez conocida V, se puede resolver la ecuacion de Schrddinger para 
la funcion de onda * de la particula, cuya densidad de probabilidad \ ^ f se 
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puede determinar para x, y, z, t. En muchas situaciones, la energia potencial 
de una particula no depende explicitamente del tiempo; las fuerzas que 
actuan sobre ella y, por lo tanto, V, vanan solamente con la posicion de 
la particula. Cuando esto se cumple, la ecuacion de Schrddinger se puede 
simplificar eliminando todo lo referente a t. Notemos que se puede escribir 
la funcion de onda unidimensional de una particula litre 

= ip{x)e-^'^/^^* (22) 

Esto es, "^{x, t) es el producto de una funcion dependiente del tiempo e~('^/^)* 
y una funcion dependiente de la posicion iIj{x, t). Sucede que las variaciones 
con el tiempo de todas las funciones de particulas, sohre las que actuan 
fuerzas estacionarias, tienen la misma forma que las de una particula libre. 
Sustituyendo la ^' de la ec. 21 en la ecuacion de Schrddinger dependiente 
del tiempo, encontramos que 

E^e-^'^/"^' = -^e-(^^/^)*^ + V^e-^'^/"^' (23) 
2m ox-' 

y, as/ dividiendo ambos miembros entre el factor exponencial comun, 

g| + ^(£-V)« = (24) 

que es la ecuacion de Schrddinger en estado estacionario. En tres dimen- 
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siones es 



En general, la ecuacion de Schrddinger en estado estacionario se puede 
resolver unicamente para algunos valores de la energia E. Lo que queremos 
decir con esto no se reRere a las diRcultades matematicas que se pueden pre- 
sentar, sino a algo mas fundamental. ''Resolver" la ecuacion de Schrddinger 
para un sistema dado signiUca obtener una funcion de onda ijj que no solo 
obedezca a la ecuacion y a las condiciones en la frontera que existan, sino 
que tambien cumplan las condiciones de una funcion de onda aceptable, es 
decir, que la funcion y su derivada scan continuas finitas y univaluadas. De 
esta manera, la cuantizacion de energia aparece en la mecanica ondulatoria 
como un elemento natural de la teoria. Asi la cuantizacion de la energia en 
el mundo fisico se ha revelado como un fenomeno universal caracteristico de 
todos los sistemas estables. 

ECUACION DE SCHRODINGER PARA EL ATOMO 
DE HIDROGENO 

A continuacion aplicaremos la ecuacion de Schrddinger al atomo de hidrogeno 
el cual esta formado por un proton, particula con carga electrica +e, y un 
electron, que tiene carga -e y que es 1,836 veces mas ligero que el proton. 

Ahora, si la interaccion entre dos particulas es de tipo u{\ fi — \) , el 
problema de movimiento de tales particulas en mecanica cuantica y tambien 



135 



en mecanica clasica se reduce al movimiento de una sola particula en el 
campo de simetna esferica, entonces tenemos el siguiente Lagrangiano: 

L = ^miT^ + ^7712/1 - u(| ri - r2 |) (26) 

Introduciendo las siguientes expresiones: 

f=ri-r2 (27) 



y 



^^mrn+m^ (28) 
mi + m2 

por lo tanto el Lagrangiano nos queda: 

L = ^MB? + - u{r) (29) 

donde 

M = mi + m2 (30) 

y 

= (31) 

mi + 777-2 

Por otro lado la introduccion del impulso se hace con las formulas de 
Lagrange 

dL 

P = ^ = MR (32) 
dR 



y 



^ dL 

p = — 7^ = mr (33) 
df 
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Jo que permite escrihir la funcion clasica de Hamilton 

Entonces se puede obtener el operador Hamiltoniano del problema cor- 
respondiente cuantico con conmutadores de tipo 

[Pi,Pk] = -ih6ik (35) 

y 

\Pi,Pk] = -iTi^ik (36) 
por lo tanto el operador Hamiltoniano es de la forma 

Este operador Hamiltoniano es la parte fundamental de la ecuacion de 

Schrddinger puesta en la forma 

Hil) = E^l; (38) 

lo cual es una forma muy practica de escribirla, pero lo mas importante 
hasta ahora escrito en esta seccion es que se ha tratado al sistema formado 
por el proton y el electron como un sistema clasico con particulas de masa 
no despreciable, como lo demuestran las ecs. 24-29, ya que no se estan 
tomando en cuenta velocidades cercanas a la de la luz, por este motivo se 
puede aplicar perfectamente la ecuacion de Schrddinger con resultados muy 
satisfactorios. 
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La ecuacion de Schrddinger para el electron en tres dimensiones, que es 
la que debemos emplear para el atomo de hidrogeno, es la ec. 21. Utilizare- 
mos esta ecuacion de Schrddinger independiente del tiempo debido a que el 
potencial V depende solamente de r y no del tiempo. 

La energia potencial V, a causa de la energia potencial electrostatica de 
una carga -e a una distancia r de otra carga +e, es 

V = --—- (39) 

Puesto que V es una funcion de r en vez de serlo de x,y,z, no pode- 
mos sustituir la ec. 39 directamente en la ec. 21. Hay dos posibilidades: 

expresar V en funcion de las coordenadas cartesianas x,y,z sustituyendo a 



r por ^/x"^ + y"^ + z^, o expresar la ecuacion de Schrddinger en funcion de 
las coordenadas polares esfericas r, 9, 4>. Haciendo esto ultimo debido a la 
simetria de la situacion fisica, el problema se simplifica considerablemente. 
Por lo tanto, en coordenadas polares esfericas, la ecuacion de Schrddinger 

es 

J.2 y Qf J r^senO 89 \ 89 J r'^sen'^9 8(p'^ 

(40) 

Sustituyendo la energia potencial V de la ec. 39 y multiplicando toda la 
ecuacion por r^sen'^9, se obtiene 

(41) 
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Esta ecuacion, es la ecuacion diferencial parcial de la funcion de onda ^{rO, (p) 
del electron en un atomo de hidrogeno. Junto con las diversas condiciones 
que ijj{r, 0, 4>) debe cumplir (por ejemplo, il^{r, 9, 0) tiene un solo valor para 
cada punto r,9,(p), esta ecuacion especifica totalmente el comportamiento 
del electron. Para ver cual es este comportamiento, resolveremos la ec. 41 
para ip{r,6,(p) e interpretaremos los resultados obtenidos. 

SEPARACION DE VARIABLES EN LA ECUACION 
DE SCHRODINGER 

Lo verdaderamente valioso de escribir la ecuacion de Schrddinger en coor- 
denadas esfericas para el problema del atomo de hidrogeno esta en que de 
esta forma se puede separar facilmente en tres ecuaciones independientes, 
cada una de ellas con una sola coordenada. El procedimiento consiste en 
buscar las soluciones en que la funcion de onda tpirjO, (j)) tiene la forma de 
un producto de tres funciones diferentes: R{r), que depende solamente de 
r; Q{0) que depende solamente de 0; y $(0) que solo depende de Esto 
es, suponemos que 

V'(r, e, (/>) = i?(r)e(e)$((/>) (42) 

La funcion R(r) describe la variacion dc la funcion dc onda ip del electron 
a lo largo de un radio vector desde el nucleo, siendo 9 y (j) constantes. La 
variacion de ip con el angulo cenital a lo largo de un meridiano de una 

esfera centrada sobre el nucleo esta descrita por la funcion Q(9) para r y (j) 
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constantes. Finalmente, la funcion <&(^) describe como vana i/j con el angulo 
azimutal (f) a lo largo de un paralelo de una esfera centrada sobre el nucleo, 
siendo r y 9 constantes. 

La ec. 42 se puede escribir mas facilmente como = de donde 

vemos que 

Al sustitulr las ecs. 43-45 en la ec. 41, que es la ecuacion de Schrodinger 
para el atomo de hidrogeno, y se divide la ecuacion total entre i?G$ se tiene 
que 

(46) 

El tercer termino de esta ecuacion solo es funcion del angulo (p, mientras que 
los otros dos son funcion der y 9. Volviendo a escribir la ecuacion anterior, 
tenemos 

sen'^9 B f 2dR\ senO B f BOX 2mr^sen'^e / \ _ 15^$ 

~^B^y 'B^)^~e~B9[ ~B9)^ l^W^+ J ~~^W 

(47) 

Esta ecuacion solamente puede ser correcta si sus dos miembros son iguales 
a la misma constante, ya que son funciones de variables diferentes. A esta 



140 



constante es conveniente Uamarla mf . La ecuacion diferencial para la funcion 
$ es 



1 



Si se sustituye mf en el segundo miembro de la ec. 47 , se divide la ecuacion 
resultante entre sen^O y se reagrupan terminos, se tiene 



(48) 



mf ^_d_( dQ\ 

sen^e Qsene 89 [ dO J ^ ^ 



Rdr \ dr J ' \47reor' 
Se tiene otra vez una ecuacion en que aparecen variables diferentes en cada 
miembro, requiriendose que ambas sean iguales a la misma constante. A 
esta constante se le Uamara, par razones que veremos mas adelante, l{l + 1). 
Las ecuaciones para las funciones Q{0) y R{r) son 



( senO— ] =l{l + 1) 



sen'^e QsenO dd \ d0 J 



y 



Id/ 2dR\ 2mr2 / 
Rd^ V 'd^J ^ ~W 



+ E 



y47reor' 

Las ecs. 48,50 y 51 se escriben normalmente como 



1(1 + 1) 



]_d_ 

j,2 



,dR 
dr 



+ 



2m / 



.2 ^ 

+ E 



m 



e = 



sen^ 
1(1 + 1) 



\^47reor' ' J 



R = 



(50) 
(51) 

(52) 
(53) 
(54) 



Cada una de estas ecuaciones es una ecuacion diferencial ordinaria de 
una funcion con una sola variable. Con ello se ha conseguido simplificar la 
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ecuacion de Schrddinger para el atomo de hidrogeno que, al principio, era 
una ecuacion diferencial parcial de una funcion ijj de tres variables. 

LOS NUMEROS CUANTICOS 
0.1 Solucion Pcira La Pcirte Azimutal 

La ec. 52 se resuelve facilmente para encontrar que su solucion es 

$(<^) = A<^e*-'^ (55) 

donde es la constante de integracion. Una de las condiciones establecidas 
previamente que debe cumplir una funcion de onda (y por lo tanto que 
es una componente de la funcion completa Tp) cs que tenga un unico valor 
para cada punto del espacio. Por ejemplo se observa que (j) y (p + 27: se 
identifican en el mismo piano meridiano. Por tanto, debe ser cierto que 
$(^) = $(^ + 27r), o bien, que Ae'"'^'^ = Ae'""^'''+^''\ lo que solamente puede 
ser cuando mi sea o un numero entero positivo o negativo (±1, ±2, ±3, ...). 
La constante mi se conoce como el numero cuantico magnetico del atomo 
de hidrogeno y gobierna a la direccion del momcnto angular L. El numero 
cuantico magnetico mi esta determinado por el numero cuantico orbital I 
que a su vez determina la magnitud del momento angular del electron. 

La interpretacion del numero cuantico orbital I no cs tan evidcntc. Ex- 
aminemos la ec. 54, que corresponde a la parte radial R{r) de la funcion 
de onda ip. Esta ecuacion esta relacionada unicamente con el aspecto radial 
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del movimiento de los electrones, es decir, con el movimiento de aproxi- 
macion y alejamiento de los mismos al nucleo; sin embargo, esta presente en 
ella la energia total del electron E. Esta energia incluye la energia cinetica 
del electron en su movimiento orbital que no tiene nada que ver con el 
movimiento radial. Esta contradiccion se puede eliminar con el siguiente 
razonamiento: la energia cinetica T del electron tiene dos partes, Tradial 
debido a su movimiento de aproximacion y alejamiento del nucleo, y Tgrrhital 
debida a su movimiento alrededor de el. La energia potencial V del electron 
es la energia electrostatica dada por la ec. 39. For io tanto, la energia total 
del electron es 



Sustituyendo esta expresion de E en la ec. 54 obtenemos, despues de rea- 
grupar terminos. 



Si los dos ultimos terminos entre corchetes de esta ecuacion se anulan entre si 
tenemos lo que necesitabamos: una ecuacion diferencial para R{r) constitu- 
ida exclusivamente por funciones del radio vector. Por lo tanto, necesitamos 
que 




(56) 




R = (57) 




(58) 



La energia cinetica orbital del electron es 




(59) 
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Puesto que el momento angular L del electron es 

L = mvorutair (60) 
podemos expresar la energia cinetica orbital 



2 



(61) 



Par la tanto, en la ec. 58 tenemos 

_ h'^iji + 1) 

2mr^ 2mr^ 

lo que nos da 



(62) 



L = ^l{l + l)n (63) 

La interpretacion de este resultado es que, puesto que el numero cuantico 
orbital I esta limitado a los valores I = 0, 1,2, (n — 1), el electron puede 
tener solamente los momentos angulares L que se especiRcan mediante la ec. 
63. Al igual que la energia total E, el momento angular se conserva y esta 
cuantizado. El termino h = h/2'K = 1.054 x 10~^^ J-s es la unidad natural 
del momento angular. 

En el movimiento planetario macroscopico, una vez mas, el numero 
cuantico que describe el momento angular es tan grande que la separacion 
en estados discretos del momento angular no se puede observar experi- 
mentalmente. Por ejemplo, un electron (o para este caso, cualquier otro 
cuerpo) cuyo numero cuantico orbital sea 2, tiene un momento angular 
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L = 2.6 X 10 J-s. Pot el contrario, el momento angular orbital de la 
tierra es /2.7 x 10^° J-s/ 

Se acostumbra designar a los estados del momento angular con la letra 
s para I = 0, con la letra p cuando I = 1, y asi sucesivamente. Este original 
codigo se origino en la clasiGcacion empirica de los espectros en las Uamadas 
series que recibieron los nombres de definida, principal, difusa y fundamen- 
tal, nombres que se les did desde antes de que se desarroUara la teoria del 
atomo. Asi un estado s es el que no tiene momento angular, un estado p 
tiene el momento angular V^h, etc. 

La combinacion del numero cuantico total con la letra que representa al 
momento angular orbital proporciona una notacion apropiada, y que es muy 
comiin para los estados atomicos. En esta notacion, por ejemplo un estado 
en el que n = 2, I = es un estado 2s y uno en el que n = 4, I = 2 es un 
estado Ad. 

Por otro lado para la interpretacion del numero cuantico magnetico, ten- 
emos que, al igual que el momento lineal, el momento angular es un vector, 
de modo que para describirlo se requiere que se especiGque su direccion, su 
sentido y su magnitud. (El vector L es perpendicular al piano en el que 
tiene lugar el movimiento de rotacion, y su direccion y sentido estan dados 
por la regla de la mano derecha: cuando los dedos apuntan en la direccion 
del movimiento, el pulgar tiene la direccion y el sentido de L.) 

iQue signi&cado posible pueden tener una direccion y un sentido en el 
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espacio para un atomo de hidrogeno ? La respuesta es sencilla si pensamos 
que un electron que gira alrededor de un nucleo es un diminuto circuito 
que, como dipolo magnetico, tiene tambien un campo magnetico. En conse- 
cuencia, un electron atomico que posee momento angular interactua con un 
campo magnetico externa B. El niimero cuantico magnetico mi especiRca la 
direccion de L, determinando la componente de L en la direccion del campo. 
Este fenomeno se conoce comunmente con el nomhre de cuantizacion espa- 
cial. 

Si hacemos que la direccion del campo magnetico sea paralela al eje z, 
la componente de L en esta direccion es 

= mih (64) 

Los valorcs posibles de m,i para un valor dado de I, van desdc +/ hasta —I, 
pasando por 0, de modo que las posibles orientaciones del vector momento 
angular L en un campo magnetico son 21 + 1. Cuando I = 0, L^ puede 
tencr solamente el valor cero; cuando I = 1, Lz puede ser h, 0, 6 — fi; cuando 
I = 2, Lz puede ser 2Ti, U, 0, —U, 6 —2h, y asi sucesivamente. Aclaremos que 
L nunca puede estar alineado exactamente (paralela o antiparalelamente) 
con B, ya que Lz es siempre mas pequeno que la magnitud + l)h del 
momento angular total. 

La cuantizacion espacial del momento angular orbital del atomo de hidrogeno 
se muestra en la fig. 6.1. 
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2h/2jl 
h/27t 



-h/2ll 
-2h/2jl 

Fig. 6.1: Cuantizacion del espaclo del mom. angular. 

Debemos considerar al atomo caracterizado por un cierto valor de mi 
como preparado para tomar una determinada orientacion de su momento 
angular L, relativo a un campo magnetico externa en el caso de encontrarse 
en el. 

En ausencia de un campo magnetico externo, la dircccion del eje z es 
completamente arbitraria. Por tanto, debe ser cierto que la componente de 
L en cualquier direccion que escojamos es m;^; el significado de un campo 
magnetico externo es que proporciona una direccion de referenda impor- 
tante experimentalmente. Un campo magnetico no es la unica direccion 
de referenda posible. Por ejemplo, la Imea entre los dos atomos H en la 
molecula de hidrogeno H2 tiene tanto significado experimental como la di- 
reccion de un campo magnetico y, a lo largo de esta Imea, las componentes 
de los momentos angulares de los atomos de H estan determinados por sus 
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m=o Estados 1= 2; L= 6 h/2jl 



valores mi. 

iPor que esta cuantizada unicamente la componente de L? La respuesta 
se relaciona estrechamente con el hecho de que L nunca puede apuntar a 
cualquier direccion z especifica; en lugar de ello describe un cono en el 
espacio, de manera que su proyeccion es miTi. La razon de este fenomeno 
es el principio de incertidumbre: si L estuviera Bjo en el espacio, de manera 
que Lx, Ly y L^ tuvieran valores de&nidos, el electron estaria conBnado en 
un piano deGnido. For ejemplo, si L estuviera en la direccion z, el electron 
tendria que estar en el piano xy todo el tiempo (fig. 6.2a). 

z 

z 

mh/2n 




Fig. 6.2: El principio de incertidumbre prohibe que el vector 
L del momento angular tenga una direccion definida 
en el espacio. 

Esto unicamente puede ocurrirsi la componente del momento del electron 
Pz en la direccion z es in&nitamente incierta, la que, par supuesto, es imposi- 
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ble si es parte de un atomo de hidrogeno. Sin embargo, como en realidad 
unicamente una componente de L junto con su magnitud L tiene valores 
deGnidos y \ L \>\ \, el electron no esta limitado a un piano unico (fig. 
6.2b), y si asi fuera, habria una fundada incertidumbre en la coordenada 
z del electron. La direccion de L cambia constantemente (Gg. 6.3) y as/ 
los valores promedio de y Ly son 0, aunque tenga siempre el valor 
especi&co rriih. 





L ^ (h/2lt) 


^ 




Z- r 


^^^^^^^^\ 































1=2 



Fig. 6.3: El vector L de mom. angular tlene preceslon constante 
en torno al eje z. 

La solucion para $ tambien debe cumplir con la condicion de normal- 
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izacion la cual esta dada por la ec. 2, entonces para $ tenemos 

r2TT 



/ I $ p = 1 (65) 
Jo 



al sustituir $ se tiene 

l2 



4# = 1 (66) 



con lo cual se tiene que = 1/ \/27r y por lo tanto $ ya normalizada esta 
dado por 

m = :^e'^^^ (67) 
Solucion Para La Parte Polar 

La ecuacion diferencial 53 para &{9) tiene una solucion mas complicada y 
esta dada por los polinomios asociados de Legendre 

(68) 

estas funciones cumplen con la relacion de ortogonalidad 

Ahora, @{9), que es la solucion para la ec. 53, esta dada por los polinomios 
de Legendre normalizados, esto es, si 

Q{e) = AeP^' {case) (70) 

entonces la condicion de normalizacion esta dada por 

Al[P[''{cose)fdcose = 1 (71) 
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por lo tanto la constante de normalizacidn para la parte polar es 




2l + l{l-miy. 
2 {l + mi)\ 



(72) 



y por consiguiente, la funcion @{9) ya normalizada es 




(73) 



Para nuestro propdsito, lo mas importante de estas funciones es que, 
como ya se dijo anteriormente, existen solamente cuando la constante I es 
un numero entero igual o mayor que | |, que es el valor absoluto de m,i. 
Esta exigencia se puede expresar como una condicion de mi en la forma 



Armonicos Esfericos 

Las soluciones para la parte azimutal y polar se puedcn juntar para formar 
los Uamados armonicos esfericos, estos dependen de <p y 9 y de alguna forma 
hacen mas facil el manejo de la funcion de onda completa ■0(r, 0, 0). Los 
armonicos esfericos estan dados de la siguiente manera: 



El factor (—1)"*' que se ha introducido sin ningun problema debido a que la 
ecuacion de Schrddinger es lineal y homogenea y ademas de que particular- 
mente es conveniente para el estudio del momento angular. Este factor es 
un factor fase Uamado fase Condon-Shortley y el efecto es para introducir 
una alternancia de signo. 



mi = 0,±1,±2,...,±Z 



(74) 




(75) 
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Solucion Para La Parte Radial 

La solucion de la ecuacion final, ec. 54, para la parte radial R{r) de la 
funcion de onda del atomo de hidrogeno tambien es complicada, y viene 
dada por los polinomios asociados de Laguerre. La ec. 54 solo se puede 
resolver cuando E es positivo o tiene uno de los valores negativos En (lo que 
signiGca que el electron esta unido al atomo), dados por 



donde n es un numero entero y se conoce como niimero cuantico principal y 
describe la cuantizacion de la energia del electron en el atomo de hidrogeno. 
Esta ecuacion es la que obtuvo Bohr para los niveles de energia del atomo 
de hidrogeno. 

Otra condicion que se debe cumplir para resolver la ec. 54, es que n, 
conocido como numero cuantico principal, sea igual o mayor que l + l. Esto 
se puede expresar como una condicion para I en la forma 




(76) 




(77) 



La ec. 54 tambien se puede poner de la siguiente forma 




^r-l{l + l) R = 



(78) 



y su solucion esta dada por los polinomios asociados de Laguerre los cuales 



cumplen con la siguiente condicion de normalizacion 




(79) 
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Por lo tanto la solucion para la ec. 78, que corresponde a la parte radial, es: 



donde p = 2r/nao y ao = f? jme^ . 

Ahora que ya tenemos las soluciones de cada una de las ecuaciones que 
solo dependen de una variable, ya podemos construir nuestra funcion de onda 
para cada estado del electron en el atomo de hidrogeno, esto es si tenemos 
que ijj{r, 9, (p) = R{r)Q{9)^{<^), entonces la funcion de onda completa es 



l 2l + l{l-mi)\ {n-l-l)l f 2 
nr,9,'P)- y 4^ {l + mi)\2n[{n + l)\r\^o) 

{ar^e-'^'I'^LlXl («^)^r (cos^)e™"^ (81) 

donde a = 2/nao. 

Utilizando los armonicos esfericos nuestra solucion queda de la siguiente 
manera 

^('■•^•^ = -/S^ (A)='^„).e-"/^LS/(")l'r'(M) 

(82) 

Esta es la solucion a la ecuacion de Schrddinger para el atomo de hidrogeno, 
la cual describe cada uno de los estados del electron. Esta funcion de onda 
por si sola no tiene interpretacion fisica como se dijo anteriormente, pero 
el cuadrado de su valor absoluto | V' P calculado para un punto y en un 
instante determinado es proporcional a la probabilidad de encontrar exper- 
imentalmente al electron ahi y en ese instante. 
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LA DENSIDAD DE PROBABILIDAD ELECTRONICA 



En el modelo de Bohr del atomo de hidrogeno, el electron gira alrededor del 
nilcleo con una trayectoria circular. Si se realizara un experimento adecuado, 
se vena que el electron estaria siempre a una distancia del nucleo r = n^ao 
(donde n es el numero cuantico de la orbita y uq = 0.53 Aes el radio de la 
orbita mas proxima al nucleo) y en el piano ecuatorial 9 = 90, mientras que 
el angulo azimutal (j) varia con el tiempo. 

La teoria cuantica del atomo de hidrogeno modihca las conclusiones del 
modelo de Bohr en dos aspectos. En primer lugar, no se pueden dar valores 
correctos der,6, (p, sino unicamente probabilidades relativas de encontrar al 
electron en un lugar dado. Esta imprecision es, por supuesto, una conse- 
cuencia de la naturaleza ondulatoria del electron. En segundo lugar, no se 
puede pensar que el electron se mueve alrededor del nucleo en un sentido 
convencional, ya que la densidad de probabilidad \ es independiente del 
tiempo y puede variar considerablemente de un lugar a otro. 

La funcion de onda del electron ip en un atomo de hidrogeno viene dada 
por = i?0$ donde R = Rni{r) describe como varia con r cuando los 
niimeros cuanticos orbital y total tienen los valores n y I; @ = Qimi{Q) 
describe a su vez la variacion de ip con 9 cuando los niimeros cuanticos 
magnetico y orbital tienen los valores I y mi; y ^ = ^mii^P) Que propor- 
ciona la variacion de i/j con (j) cuando el numero cuantico magnetico es mi. 
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Entonces, la densidad de probabilidad \ i/j {"^ se puede escribir como 

I ^ |2 = | |2| e |2| $ |2 (-g3^ 

donde se comprende que si la funcion es compleja, hay que tener en cuenta 
que su cuadrado se debe sustituir por el producto de ella y su conjugada 
compleja. 

La densidad de probabilidad \ ^ \'^, que mide la posibilidad de encontrar 
al electron con un angulo azimutal (j) dado, es una constante que no depende 
para nada de Por lo tanto, la densidad de probabilidad del electron es 
simetrica respecto al eje de las z, independientemente del estado cuantico, de 
manera que el electron tiene igual oportunidad de encontrarse en un angulo 
(j) como en otro. 

La parte radial R de la funcion de onda, en contraste con no solamente 
var/a con r, sino que lo hace de una manera diferente para cada combinacion 
de numeros cuanticos n y I. La fig. 6.4 muestra graficas de R en funcion de 
r para los estados Is, 2s, y 2p del atomo de hidrogeno. Evidcntcmcntc, R es 
mebcimo al ser r = (esto es, en el nucleo mismo) para todos los estados s, 
mientras que es cero en r = para todos los estados que poseen momento 
angular. 
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Fig. 6.4: Variacion con la distancia al nucleo 
de las funciones de onda radiales 
Is, 2s, 2p; ao° 0.53 angs. (radio Bohr) 




5 10 15 r(a ) 



Fig. 6.5: Probabilidad de encontrar un electron 
del atomo de hidrogeno entre r y r+dr 

con respecto al nucleo para Is, 2s, 2p. 



La densidad de probabilidad del electron en el punto r, 6, (p es propor- 
cional a \ ijj pero la probabilidad real de encontrarlo en el elemento de 
volumen infinitesimal dV es \ ip p dV . Ahora, en coordenadas polares 
esfericas 

dV = r^senedrd9d(l) (84) 
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de manera que, como O y # son funciones normalizadas, la prohahilidad 
numerica real P{r)dr de encontrar al electron en el atomo de hidrogeno, a 
una distancia comprendida entre r y r + dr del nucleo, es 



Esta ecuacion esta representada en la fig. 6.5 para los mismos estados cuyas 
funciones radiales R aparecen en la fig. 6.4; en principio, las curvas son 
completamente diferentes. Observamos de inmediato que P no es maximo 
en el nucleo para los estados s, como lo es R, sino que tiene su maximo 
a una distancia finita de el. El valor mas probable de r para un electron 
Is es exactamente qq, que es el radio de la orbita del electron en estado 
fundamental en el modelo de Bohr. Sin embargo, el valor medio de r para 
un electron Is es 1.5ao, lo cual parece enigmatico a primera vista, ya que 
los niveles de energia son los mismos en mecanica cuantica y en el modelo 
atomico de Bohr. Esta aparente discrepancia se elimina cuando se tiene en 
cuenta que la energia del electron depende del/r y no directamente de r, y 
el valor medio del/r para un electron Is es precisamente l/ao- 

La funcion Q varia con el angulo polar 9 para todos los numeros cuanticos 
I y mi, excepto para I = mi = 0, que son estados s. La densidad de proba- 
bilidad \ p para un estado s, es una constante (1/2), lo que significa que, 
como I $ p es tambien constante, la densidad de probabilidad electronica 




= \R\'^ dr 



(85) 
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I V P tiene el mismo valor, para un valor de r dado, en todas las direc- 
ciones. Los electrones en otros estados tienen preferencias angulares que 
algunas veces Uegan a ser muy complicadas. Esto se puede observar en la 
fig. 6.5, donde se muestran, para varios estados atomicos, las densidades de 
probabilidad electronica en funcion de r y 9. (El termino que se representa 
es \ ijj , y no \ ijj dV.) Puesto que \ ijj es independiente de 4>, una 
representacion tridimensional de \ tp se obtiene haciendo girar una repre- 
sentacion particular alrededor de un eje vertical. Al hacerlo, las densidades 
de probabilidad para los estados s tienen evidentemente simetria esferica, 
mientras que las otras no la tienen. Los tipos de lobulos pronunciados, 
caracteristicos de muchos de los estados, tienen importancia en qmmica ya 
que estos modelos ayudan a determinar la manera como interactuan en las 
moleculas los atomos adyacentes. 

Notas: 

1. E. Schrddinger consigio el premio Nobel en 1933 (junto con Dirac) 
por "el descubrimiento de nuevas formas productivas de la teoria atomica". 
Schrddinger escribio una serie notable de cuatro articulos intitulada "Quan- 
tisierung als Eigenwertproblem" (I-IV, recibidos por Annalen der Physik el 
27 de Enero, el 23 de Febrero, el 10 de Mayo y el 21 de Junio de 1926). 
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Problemas 

Problema 6.1 - Obtener la ecuaciones para las orbitas estables y para 
los niveles de energia del electron en el atomo de hidrogeno. 

Respuesta; Tenemos que la longitud de onda del electron esta dada 

por 

mv 

mientras que por otro lado al igualar la fuerza electrica con la fuerza centripeta, 
esto es 

2 1 2 

mv 1 e 



r 47reo 

nosotros obtenemos que la velocidad del electron esta dada por 



V 



•\/47reo mr 

entonces la longitud de onda del electron es 



h j 47reor 

A 



I \ m 

Ahora, si damos el valor 5.3 x 10~^^m al radio r de la orbita electronica, 
vemos que la longitud de onda del electron es A = 33 x 10~^^m. Esta 
longitud de onda tiene exactamente el mismo valor que la circunferencia de 
la orbita del electron, 27rr = 33 x 10~^^m. Como se puede ver, la orbita de 
un electron en un atomo de hidrogeno corresponde asi a una onda completa 
cerrada sobre si misma. Esto se puede comparar con las vibraciones de un 
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anillo de alambre, si las longitudes de onda estan en un niimero entero de 
veces en su circunferencia este podra seguir vibrando indefinidamente, pero 
si un numero no entero de longitudes de onda tiene lugar sobre el anillo se 
producira una interferencia negativa a medida que las ondas se desplacen 
en torno a el, y las vibraciones desapareceran rapidamente. Con esto se 
puede afirmar que un electron puede girar indefinidamente alrededor de un 
nucleo sin irradiar energia con tal que su orbita contenga un numero entero 
de longitudes de la onda de De Broglie. Con esto tenemos que la condicion 
de estabilidad es 



donde r„ designa el radio de la orbita que contiene n longitudes de onda. 
Al sustituir A tenemos 



Para los niveles de energia tenemos que E = T + V, al sustituir las 
energias potencial y cineticas obtenemos 



nX = 2TTr, 



n 




por lo tanto las orbitas estables del electron son 




Trme^ 



J. 2 

E = -mv — 
2 
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o lo que es igual 

En = — - 



87reor-„ 

al sustituir el valor de r„ en esta ultima ecuacion obtenemos 



En = 



ill' ' ( 1 



lo que nos da los niveles de energia. 



Problema 6.2 - El teorema de Unsold dice que, para cualquier valor 

del numcro cuantico orbital I, las dcnsidadcs dc probabilidad, sumadas para 
todos los estados posibles, desde mi = —I hasta mi = +1 da una constante 

independiente de los angulos 6 o (f) esto es 

+1 

mi=-l 

Este teorema signiGca que todo atomo o ion con subcapa cerrada tiene 
una distribucion simetrica esferica de carga eletrica. Comprobar el teorema 
de Unsold para 1 = 0, l = lyl = 2. 

Respuesta; Tenemos que para I = 0, Gqo = l/\/2 y = I/V^tt por 
lo tanto del teorema de Unsold vemos que 

I eo,o n $0 p= I 

Para 1 = 1 tenemos que 

+1 

E I ®lrnt I'l l' = l ©1,-1 I'l ^-1 I' + I ©1,0 I'l $0 I' + I ©1,1 I'l $1 I' 

mi=-l 

161 



por otro lado las funciones de onda son: = {■s/^/2)sen9, #_i = 

{l/^/2^)e-^'t>, 61,0 = {^/^/2)cose, $0 = V^, 61,1 = (V3/2)sen0, $1 = 
(l/\/27r)e*'^ que al sustituirlas en la ecuacion anterior nos queda 

mi=-l 

lo que demuestra que tambien es una constante. 
Para I = 2 tenemos que 

+2 

mi=-2 



62,-2 I'l $-2 I'l 62,-1 I'l $-1 I' + I 62,0 I'l $0 I' + I 62,1 I'l $1 I' + I 62,2 



y las funciones de onda son: 62,-2 = (\/l5/4)sen^6', $_2 = (l/-\/27r)e 
62,-1 = (^/2)sen0cos0, $-1 = {l/V2^)e-^'t', 62,0 = {VT0/A){3cos'^e - 
1), $0 = l/V2n, ©2,1 = (\/T5/2)sen6'cos6', ^>i = {1/V2^)e'^, 62,2 = 
(\/l5/4)sen^^, $2 = (l/A/27r)e^*'^, sustituyendo en la ecuacion anterior nos 
queda 

+2 



_5_ 
47r 



mj=-2 

con lo que queda demostrado el teorema de Unsold. 



Problema 6.3 - La probabilidad de encontrar un electron atomico cuya 
funcion de onda radial sea R{r), fuera de una esfera de radio ro centrada en 
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el niicleo, es 

POO 

/ I R{r) |2 r^dr 

Jro 

La funcion de onda Rioir) corrcsponde al estado fundamental de un atomo 
de hidrogeno, y oq es el radio de la orbita de Bohr correspondiente a este 
estado. Calcular la probabiUdad de encontrar un electron en estado fun- 
damental en un atomo de hidrogeno a una distancia del nucleo mayor que 
Oo- 

Respuesta; Tenemos que la funcion de onda radial que corresponde al 
estado fundamental es 

Rioir) = -4j2^-''^"° 
% 

sustituyendo en la integral nos queda 



poo 

/ l^( 

Jao 

o lo que es igual 

poo 4 

/ I R{r) |2 r^dr = 

Jao Oq 



r) 1^ r^dr = —3 



2 3 " 

^ 2 -2r-/ao _ ^„g-2r-/ao _ ^ -2r/ao 

2 2 4 



00 



esto nos da como resultado al evaluar 



fOO 

/ \R(.r) 

Jao 



P r^dr = 



que es la probabihdad de encontrar al electron. 
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7. DISPERSION EN LA MC 



Introduccion 

Para la teona cuantica de dispersion nos ayudaremos de los resultados ya 
conocidos de la dispersion en campos de fuerzas centrales, y asumiremos 
ciertas situaciones que simplificaran los calculos, si bien no nos alejaran 
demasiado del problema "real". Sabemos que en el analisis experimental de 
una colision podemos obtener datos que nos ayuden a entender el estado 
de la materia "bianco", o bien, la interaccion entre el haz incidente y el 
"bianco". Las hipotesis que asumiremos son: 

i) Asumiremos que las particulas no tienen espm, lo cual no implica que 
este no sea importante en la dispersion. 

a) Nos ocuparemos solo de la dispersion elastica, en la cual no consider- 
amos la posible estructura interna de las particulas. 

Hi) Asumiremos que el bianco es lo su&cientemente delgado como para 
no permitir la dispersion multiple. 

iv) Asumiremos que las interacciones son descritas por un potencial que 
depende solo de la posicion relativa de las particulas. 

Asi, usaremos los resultados ya conocidos de la teoria de dispersion de 
la mecanica clasica, ahi se define: 

do_ ^ I{e,ip) 

dn lo ^ ' 
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donde da es el elemento de angulo solido, Iq es el numero de parti'culas 
incidentes por unidad de area, e Idfl el numero de particulas dispersadas en 
el elemento de angulo solido. 
Con esto, Uegamos a que: 



Si deseamos conocer en terminos de mecanica cuantica lo que sucede, 
debemos estudiar la evolucion en el tiempo de un paquete de ondas. Sea 
ahora Fi el Aujo de particulas del haz incidente, es decir, el numero de 
particulas por unidad de tiempo que atraviesan una superRcie unitaria per- 
pendicular al eje. Colocamos un detector lejos de la region de inRuencia del 
potencial, y que subtiende un angulo solido d^l; con esto, podemos contar 
el numero de particulas dn dispersadas por unidad de tiempo en dCl en la 
direccion {9, if). 



da 



_P_,dp 
sinO^de 



(2) 
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dn es proporcional a dfi y a Fj. Llamaremos a{6,(p) al coeficiente de 
proporcionalidad entre dn y FidO,: 

dn = a{e, ip)FidVL , (3) 

esta es la Uamada seccion diferencial transversal. 

El numero de particulas por unidad de tiempo que alcanzan el detector 
es igual al numero de particulas que cruzarian una superBcie a{6, ip)dQ, colo- 
cada perpendicular al eje del haz. La seccion total de dispersion se define 
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por: 

a = J a{e, (p)dn (4) 

Ya que podemos orientar los ejes de coordenadas a nuestra eleccion, lo 
haremos de tal forma que el eje del haz incidente de particulas coincida con, 
digamos, el eje z (por simplicidad en los calculos, donde usaremos coorde- 
nadas esfericas). 

En la region negativa del eje, para t negativo grande, la particula sera 
practicamente libre; no se vera afectada por V{r), y su estado se puede 
representar por ondas planas, por lo tanto la funcion de onda debe contener 
terminos de la forma e^^^ , donde k es la constante que aparece en la ecuacion 
de Helmholtz. Por analogia con optica, la forma de la onda dispersada es: 

pikr 

fir) = — (5) 



en efecto, pues: 



y 



(V^ + e)e''''' / (6) 



„ikr 

{V^ + k'^)^ = (7) 



para r > ro, donde rg es cualquier numero positivo. 

Asumamos que el movimiento de la particula esta descrito por el hamil- 
toniano: 
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H = ^ + V = Ho + V 



(8) 



V es diferente de cero solo para una pequena vecindad del origen. Sabe- 
mos la evolucion y descripcion de un paquete de ondas ent = 0; 



donde, como ya se sabe, ip es una funcion de ancho Ak, centrada alrede- 
dor de ko- Asumimos tambien que ko es paralela a tq pero en direccion 
opuesta. Para averiguar que sucede con el paquete de onda-que es lo que 
nos interesa-cuando en un tiempo posterior el paquete ha chocado con el 
bianco y ha sido dispersado por este, podemos valernos de la expansion de 
4>{r,0) en terminos de las eigenfunciones tpni'"^), de H. Asi, podemos ex- 
pandir: ip{r, 0) = J2n CnV'n(r)- De esta forma, el paquete de onda al tiempo 
t es: 



Esta es una eigenfuncion del operador Hq y no de H, pero podemos 
sustituir estas eigenfunciones por eigenfunciones particulares de H, que des- 
ignaremos por ip^'^ (r) . La forma asintotica de estas ultimas son de la forma: 






n 



(10) 



V'[+^(r)~e'^"- + f(r)^ 



(11) 



donde, como es usual, p = 



fikyE 



2m 
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Esto corresponde a una onda plana para el haz incidente, y otra onda 
esferica divergente, como resultado de la interaccion entre el haz y el bianco. 
Estas soluciones de la ecuacion de Schrddinger existen en realidad, y pode- 
mos expandir V'(r, 0) en terminos de ondas planas o en terminos de V'fc(r)-' 



Esto puede verse entonces como que la onda esferica divergente no tiene 
contribucion alguna al paquete de ondas inicial. 

Dispersion 

Al viajar la onda, necesariamente sufre una dispersion, de manera que aqm 
no podemos ya pasar por alto el efecto de la onda divergente, pero podemos 
usar lo siguiente: 



^ = i-k' = i-[l^o + (k - ko)]^ = i-[2ko • k - kg + (k - ko)^] , (13) 



para poder despreciar el ultimo termino del parentesis, al sustituir uj en la 
relacion para if), requerimos que: ^(k — ko)^r ^ 1, donde T 2± con 
lo que obtenemos: 




(12) 



donde hu) 



2m 



(Afe)Vo 
ko 



< 1 



(14) 
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Esta condicion implica que el paquete de onda no se dispersa de manera 
apreciable al desplazarse una distancia macroscopica tq. 

Escogiendo la direccion del vector k de la onda incidente aUneado con el 

sistema de coordenadas esfericas, podemos escribir: 

Ya que H, el operador hamiltoniano (que hemos considerado hasta ahora 
no como operador, pero cuyos resultados son los mismos) es invariante ante 
rotaciones en el eje z, podemos escoger condiciones de frontera que tambien 
sean invariantes, de manera que: 

estas se conocen como funciones de onda de choque. El coe&ciente f{9) se 
conoce como amplitud de choque. 

Amplitud de probabilidad 

Recordemos ahora la ecuacion de Schrodinger a resolver: 
Y, tenemos 

Pir,t) = rir,tmr,t) = \i;ir,t)\^ (16) 
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se interpreta, de acuerdo a Max Born, como una densidad de prohabilidad. 
Esta funcion de onda debe estar normalizada de manera tal que: 

J \tjj{r,t)\'^d^r = l . (17) 

El coeficiente de normalizacion de ip, y la integral de normalizacion 
deben ser independientes del tiempo, si ha de cumplir con la ecuacion de 
Schrddinger. Esto lo podemos notar de la siguiente manera: 

,.||^p(,.,,,3,.|^(,.^^^^,„V (18) 

de la ecuacion de Schrddinger: 
entonces: 

= ^ ^^[V'^VV' - {Vr )^]ndA (20) 

donde se ha usado el teorema de Green para evaluar la integral sobre el 
volumen. A es la superhcie que limita la region de integracion y [ denota 
la componente en la direccion normal a la superhcie dA. 
De&niendo: 

s{v,t) = ^[rvi^-{vr)i^] (21) 

obtenemos: 

1 = ^1 P(r, t)d^r = - / V • SdV = - / SndA (22) 
ot Jq Jn J A 
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para paquetes de onda en los que se hace cero a grandes distancias y la 
integral de normalizacion converge, la integral de superGcie es cero cuando 
U es todo el espacio. Se puede demostrar (vease P. Dennery y A. Krzywicki, 
Mathematical methods for physicists) que la integral de superticie es cero, 
de manera que la integral de normalizacion es constante en el tiempo, y se 
cumple el requerimiento inicial. De la misma ecuacion para S obtenemos: 

^^|^ + V-S(r,t)=0 (23) 

esto nos recuerda la ecuacion de continuidad dc clcctrodinamica, en este 
caso con un Hujo de densidad P y corriente de densidad S, sin fuentes o 
sumideros. Asi, es razonable interpretar S como una densidad de corriente 
de probabilidad. Por semejanza con la electrodinamica, es el operador 
velocidad y: 

S(r,t) = Rei^P*—V^) (24) 
im 

Funcion de Green en teorfa de dispersion 

Otra forma de escribir la ecuacion de Schrddinger a resolver es (— + 
F)V' = E^Jj o (V^ + = donde: F = 2pE^ ^ jj ^ 

Resulta mas conveniente transformar esta ecuacion a una de forma in- 
tegral. Esto podemos hacerlo si consideramos a Uijj del lado derecho de la 
ecuacion como una inhomogeneidad, y de esta manera una solucion de la 
ecuacion se construye con las funciones de Green, por medio de: 

(V^ + A;2)G(r, r') = S{r - r') (25) 
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entonces, una solucion a la ecuacion de Schrodinger, se obtendra de la suma 
de la solucion a la ecuacion liomogenea y la solucion a la parte inhomogenea: 

V'(r) = A(r) - J G(r, r')U(r')VXr')d^r (26) 

Buscanios una funcion G que sea producto de funciones linealmente in- 
dependientes, como lo son las ondas planas: 

G(r,r' = y"A(q)e^'i-(^-'"')(ig (27) 



usando la ecuacion 21 : 



J A{ci){e - g2)e^^-('-'-')dg = 5{r - r') (28) 
lo cual se cuniple como identidad si: 

A{cD = i2n)-^ {k^-qT' (29) 

de lo que resulta: 

1 r piqR 
= (2^ J • 

con R = \r — r'\. Tras un calculo usando los metodos de variable compleja 

^ Uegamos a: 

1 pikr 

G{r) = - — — (31) 
47r r 

£^sta funcion no esta determinada de manera univoca; la funcion de 
Green puede ser cualquier solucion de la ecuacion p3 y existen una intinidad 



Vease el problema 7.1 
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de ellas; por lo tanto, la eleccion de una en particular impone condiciones a 
la frontera sobre las eigenfunciones ipk (r) ■ 

La funcion de Green obtenida es por lo tanto de la forma: 

De esta forma, Uegamos a la ecuacion integral para la funcion de onda 
de choque: 

^Pik, r) = <f{k, r) - — 2 / r^ir'Mk, r)dr , (33) 

zTTn J r — r' 

donde es una solucion de la ecuacion de Helmholtz. Haciendo: |r — r'| = 
R, 

(V^ + k^)^ = (V^ + k'^)[ip + J G(r,r')C/(r')^(r')dV'] (34) 

asumiendo que podemos intercambiar el orden y poner el operador V dentro 
de la integral: 

(V^ + A;2)V' = /(V^ + k'^)G{r, r')U{r')^{r')d^r' = U{v)i^{r) (35) 
y se veriHca entonces que G{R) = ^^-^ es solucion. 

Teorema optico 

La seccion diferencial total esta dada por: 

f da 



crtot 



(k)^j^<ia (36) 
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Exprcsamos ahora f{0) expresada en terminos del corrimiento de fase 
Si{k) = e2^'5*(^) de forma que: 

1 oo 

/(^) = T E(2^ + 1)^'^'^'^ sm6i{k)Piicose) (37) 



entonces 



[- ^(2/+l)e*^'('=) smdiik)Pi{cose)][ / [^(2r+l)e*^''« sin 6v{k)Pv {cos 

(38) 

Usando ahora J Pi {cos 6) Pi' {cos 6) = 21^^11' obtenemos 

(Ttot = ^J2(^l + I) sin 6i{kf . (39) 



1=0 

Lo que aqm nos interesa es que: 



-I -1 1 

Im/(0) = -5^(2/+l)Im[e^^'«sin(5Kfc)]PKl) = t E(2^+l) ^i^*^' = ^^^t 

(40) 

Esta relacion es conocida como el teorema optico. Su signiEcado fisico 
es que la interferencia de la onda incidente con la onda dispersada en angulo 
cero produce la salida de la particula de la onda incidente, lo que permite 
la conservacion de la probabilidad. 

Aproximacion de Born 

Consideremos la siguiente situacion (Fig. 7.2). 
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Fig. 7.2 



Nos situamos en un punto muy alejado de P, que corresponde a la region 
de inAuencia del potencial U, y r L r' <^ I. La longitud MP, que 
corresponde a |r — r'| es aproximadamente igual a la proyeccion de MP en 
MO: 

|r — r'l c± 9 — u • r' (41) 
donde u es el vector unitario en la direccion r. Entonces, para r grande: 

1 „jfc|r-r'| 1 „ikr 
G = J. ^r^oo 6-^^="- . (42) 

47r |r — r I An r 

Sustituimos abora en la expresion integral para la funcion de onda de cboque, 
y obtenemos: 

Mr) = e'^^ -—— [ e-'^'^-''Uir')i)iY')(fr' . (43) 
47r r J 

Esta ya no es una funcion de la distancia r = OM, sino solamente deO yip, 
y entonces: 

f{e,iP) = J e"^*^" •'■[/(r')V'(r')rfV' (44) 
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DeEnimos ahora el vector de onda incidente kj como un vector de modulo 
k dirigido a lo largo del eje polar del haz tal que: e^^^ = e*''"''"; de manera 
similar, kj, con modulo k y con direccion Bjada por 9 y (p se llama vector 
de onda desplazada en la direccion {6, tp): kj = ku 

El vector de onda transferido en la direccion (9, ip) se define como: K = 
kd - ki 




Fig. 7.3 



Con esto, podemos escribir la ecuacion integral de dispersion como: 

^(r) = e^'^'-'" + J G{r,r')U{r')i>{r')d^r' (45) 

Ahora podemos intentar resolver esta ecuacion por iteracion. Hacemos 
el camhio r — > r'; r' ^ r" y con esto escribimos: 

^{r') = e^'''-'"' + j G{r',r")U{r")i^{r")d^r" (46) 

Sustituyendo esta expresion en ^ obtenemos: 

^P(r) = 6^^^-'^+ J G{r,r')U{r')e^^'-'''d\'+ J J G{r,r')U{r')G{r' ,r")U{r")i;{r")d\"d^r' 

(47) 

Los dos primeros terminos de lado derecho de la ecuacion son conocidos y 
solamente el tercero contiene la funcion desconocida ip{r). Podemos repetir 
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este procedimiento: camhiando r por r" y r' por r'" nos da ipir") la cual 
podemos reinsertar en la ecuacion 41 con esto obtenemos: 

^P(r) = e*'^'-'+ J G{r,r')U{rY^'-'''+ J J G{r,r')U{r')G{r' ,r")U{r")e'^'-''" d^r'd^r" 

+ JJ |G(r,rO[/(r')G(r',r'')t/(r'')e*'^'-'-''G(r'',r''')[/(r''')^(r''') . (48) 

Los primeros tres terminos son conocidos; la funcion desconocida ip{v) se 
ha ido hasta el cuarto termino. De esta forma, por iteracion construimos 
la funcion de onda de dispersion estacionaria. Notese que cada termino de 
la expansion Ueva una potencia mayor del potencial que la que le precede. 
Podemos continuar de esta forma hasta obtener una expresion despreciable 
del lado derecho, y obtenemos ip{r) en terminos de cantidades conocidas. 

Sustituyendo la expresion de ip{v) en f{9, ip) obtenemos la expansion de 
Born de la amplitud de dispersion. Limitandonos al primer orden en U, 
todo lo que hay que hacer es sustituir ijj^r') por e*''' '" en el lado derecho de 
la ecuacion, con esto obtenemos: 

Att J 4tt J 

e-i'^-^'U{r')d^r' (49) 

K es el vector de onda dispersada definido anteriormente. La seccion de 
dispersion se relaciona entonces de manera muy sencilla a la transformada de 
Fourier del potencial, recordando: V{r) = ^U{r), ya{9,(p) = \ f{6,ip)\'^ 
tenemos que: 



-1 

47r 



^'""KO^V) = ^1 / e-^^-^V{r)dM' (50) 
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La direccion y modulo del vector de onda dispersada K depende del 
modulo A; de ki 7 kd 7 de la direccion de dispersion {0,(p). Para un 9 y 
if, la seccion eficaz varia con k, la energia del haz incidente, y de manera 
analoga, con una energia dada a^^^ varia con y ip. Con esta aproximacion 
de Born estudiando la variacion de la seccion diferencial eficaz en terminos 
de la direccion de dispersion y la energia incidente nos da informacion del 
potencial V{r). 

Notas: Uno de los primeros trabajos de dispersion cuantica es: 

M. Born, "Quantenmechanik der Stossvorgange" , Zf. f. Physik 37, 863-867 

(1926) 
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Problemas 

Problema 7.1 

Calculo por medio de variable compleja de la funcion de Green 

Recuerdese que obtuvimos: 

G{r, r') = j^^^ J -j^z^d^q , con i? = |r — r'| . Ya que <fq = sin OdqdOdcf), 
Uegamos, despues de la integracion angular, a que: 

G{r,v') = ^j^ ^-'';::p qdq. 

Hacemos: C = -^^z^; J dividimos la integral en dos partes: 

CiJZo f^<ldQ - IZo l^^dq) . 
Hacemos ahora q —q en la primera integral: 

I-oo W^^{-l)d{-q) = /~°° W^Qdq = - JZo l^Qdq 
por lo que: 

G{r,r') = -2CiJZ.g^dq) 
sustituyendo C, obtenemos: 

Gir,r') = ^JZo§^dq 

Ahora, esta integral podemos evaluarla por los residues que posee, us- 
ando los metodos de variable compleja. Notamos que existen polos simples 
en q =1 k. 



180 




-fc-ie ]c_ie 



Fig. 7.4: Reglas de camino alrededor de los poles para G^. y G_ • 

Usamos el contorno mostrado en la Ggura 7.4, el cual rodee a los polos 
de la manera senalada, lo cual dara el efecto fisico buscado, pues notemos 
que, de acuerdo al teorema del residuo, 

G(r) = -4^^ (lm^>0), 

G(r) = -^^ (lmfc<0) 

La solucion que nos interesa es la primera, porque da ondas dispersadas 
divergentes, mientras que la otra solucion representa ondas dispersadas con- 
vergentes, y aun mas, la combinacion lineal 

corresponde a ondas estacionarias. 

La evaluacion formal de la integral se puede hacer tomando k'^ — ^ 
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k"^ +ie-q^ , de tal manera que: J^^ w^dq (pf^p^^? ■ 

Esto es posihle porque R > 0, de manera que el contorno a calcular estara 
en el semiplano complejo superior. Asi, los polos del integrando estan en: 
q Vk'^ + ie ~i + proceso hmite cuando e — > debe hacerse 

despues de la evaluacion de la integral. 

Problema 7.2 
Forma asintotica de la expresion radial 

Como ya se vio antesQ la parte radial de la ecuacion de Schrddinger se 
puede escribir como: 

ii, + '^i)RnUr) - |?[^(r) + ^-^^]R,Ur) + ^-^Rnlrn{r) = 
n,l,m son los numeros cuanticos. De aqm en adelante se eliminaran por 
comodidad, y donde se sabe que R depende solo de r. Asumiremos que los 
potenciales decaen a cero mas rapido que 1/r, y: limr-^o r'^V{r) = 0. 

Usamos ahora u{r) = rR, y como: + = ^-^^ tenemos que 

|^.+ |?[^-y(r)-iig^]n = 0. 
Notese que el potencial ahora tiene un termino mas: 

V{r) ^ V{r) + ^(g^ 
que corresponde a una barrera repulsiva centrifuga. Para una particula libre 
V{r) = {) y la ecuacion toma la forma 

[S. + ',i)-'-^]R + k^ R = 0. 

^^Vease El dtomo de hidrogeno 
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Introduciendo la variable p = kr, ohtenemos 

^ + 2^-^fl + fl = . 
dp'^ p dp 

Las soluciones a esta ecuacion son las Uamadas funciones esfericas de 
Bessel. La solucion regular es: 

y la solucion irregular: 

Mp) = -{-pn\in^) 

Para p grande, las funciones de interes son las funciones esfericas de 
Hankel: 

hi^\p) =ji{p) + mi{p) y 
hf\p) = [h'^Xp)]* 

De especial interes es el comportamiento para p ^ I, cuyo compor- 
taniiento asintotico es: 

3i{p) ^ -sin(p- ^) 
P 2 

1 It: 

ni{p)^ — cos(p-— ) (51) 
P 2 

y entonces 

I p 

La solucion regular en el origen es: Ri{r) = ji{kr) 
La forma asintotica es, usando la ecuacion ^ 

Problema 7.3 
La aproximacion de Born para potencial de Yukawa 
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Consideramos un potencial de la forma: 



V{r) = (52) 

con Vo y a constantes reales y a positiva. El potencial es atractivo o re- 
pulsivo dependiendo de si Vq es negativo o positivo; entre mas grande \Vo\, 
mas intense el potencial. Asumimos que |Vb| es su&cientemente pequeno 
para que la aproximacion de Born sea valida. De acuerdo a la formula ya 
obtenida antes en el apartado de amplitud de probabilidad, la amplitud de 
dispersion esta dada por: 

Como este potencial depende solo de r las integraciones angulares pueden 
hacerse facilmente, y Uegamos a la forma: 

Con esto Uegamos a: 
/(-)(^,^) = -^^. 

De la figura se observa que: |K| = 2fcsin |; por lo tanto: 



^^[a2+4fe2sin|'']2 



La seccion de dispersion total se obtiene por integracion: 
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8. LAS ONDAS PARCIALES 



Explicamos brevemente en que consiste el metodo de ondas parciales en el 
estudio de problemas de dispersion. 

Introduccion. 

El problema de dispersion desde el punto de vista cuantico, consiste en tratar 
a una particula que interacciona con otra llamada dispersor (en la presente 
exposicion supondremos que el dispersor siempre se encuentra fijo) en una 
region muy pequena del espacio. Fuera de esta region, la interaccion entre 
ambas particulas es despreciable. De esta manera es posible describir a la 
particula dispersada por el siguiente Hamiltoniano: 



Donde Hq corresponde al hamiltoniano para la particula libre. Entonces 
nuestro problema consiste en resolver la siguiente ecuacion: 



Es evidente que el espectro sera contmuo (estamos tratando el caso de 
dispersion elastica). La solucion a la ecuacion anterior esta dada por: 



De un ligero analisis podemos ver que en el caso que V = obtenemos 
la solucion: \ 4>), es decir, la solucion correspondiente al caso de la particula 



H = Ho + V 



(1) 



{Ho + V)\iP) =E\ V) 



(2) 




(3) 
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libre. Hay que notar que el operador -^zrjj^ en cierto sentido es anomalo, 
pues tiene un continuo de polos en el eje real que coinciden con las valores 
propios de Hq, para "librarnos" de ese problema induzcamos un pequeno 
desplazamiento del carte que yace sobre el eje real, de esta manera tenemos: 

La ecuacion anterior es conocida como la ecuacion de Lippmann-Schwinger. 
Al final el desplazamiento de los polos sera en el sentido positivo de el eje 
imaginario (para que el principio de causalidad no se viole [segun Feynman]). 
Tomemos la x-representacion: 

(X I = (X I ^) + I d'x (x| ^_^^^.J x^ (x I F I i;^) (5) 

En el primer termino del lado derecho de la ecuacion anterior vemos que 
corresponde a una particula libre y el segundo termino se interpreta como 
una onda esferica que emerge del dispersor. El kernel de la integral anterior 
lo podemos asociar con una funcion de Green o propagador y es muy sencillo 
calcular: 

G±(x,x) = ^/x| I |x-\=-i-f (6) 

' ' 2m\'^-i^o±^£ / 47r | x - x' | ' 

donde E = Ti^k"^ /2m. Como veiamos anteriormente la funcion de onda 
la podemos escribir como una onda plana mas una onda esferica que emana 
del dispersor (salvo un factor constante): 

(x| V+) = e^-^ + ^/(k,k') (7) 
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A la cantidad /(k, k ) que aparece en la ecuacion 7 se le conoce como 
amplitud de dispersion y expUcitamente la podemos escribir como: 

/(k,k') = -i-(27r)3|f(k'|F|V'+) (8) 

Definamos ahora un operador T tal que: 

T\,p) = V\ij+) (9) 

Si multiplicamos la ecuacion de Lippman-Schwinger por V y a partir de 
la deRnicion anterior obtenemos: 

T\^) = V\4>) + V I T 1 4>) (10) 
III — riQ + le 

Asi iterando la ecuacion anterior (como en teoria de perturbaciones) 
podemos obtener la aproximacion de Born y sus correcciones de orden su- 
perior. 

El metodo de ondas parciales. 

Ahora considcrcmos el caso de un potencial central no nulo. Entonccs de 
la definicion (9) del operador T se deduce que conmuta con y con L de 
ahi que se diga que T es un operador escalar. De esta manera para facilitar 
los calculos es convcnicnte utilizar coordcnadas esfericas, puesto que dada 
la simetria, el operador T sera diagonal. Ahora veamos que forma adquiere 
la expresion (8) para la amplitud de dispersion: 

, 'n 1 2m ,„ , q 
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m2Y.J2 f f dE'{k! \ E'l'm'){E'l'm' \ T \ Elm) {Elm \ k) (11) 

I I' m' 

despues de algunos calculos podemos obtener: 

/(k,k') = -ij!^^rK^)yr(kyr*(k) (12) 

I m 

Escogiendo el sistema de coordenadas tal que k tenga la misma direccion 
que el eje orientado z, de esta manera a la amplitud de dispersion unicamente 
contribuiran los armonicos esfericos con m igual a cero; si deGnimos que 9 
sea el angulo entre k 7 k' tenemos: 



Yn^') = ^^Pi{cose) (13) 
hagamos la siguiente deGnicion: 

(14) 

as/ la ecn. (12) se puede escribir de la siguiente manera: 

00 

/(k,k') = f{e) = ^(2z + i)fi{k)Pi{cose) (15) 
1=0 

A la cantidad fi (k) le podemos dar una interpretacion sencilla a partir del 
desarrollo de una onda plana en ondas esfericas, veamos el comportamiento 
de la funcion (x | para grandes valores de r, que como ya habiamos 
establecido previamente debe tener la forma: 



(x I - ^ 



(27r)3/2 



pikr 



r 
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(27r)3/2 



1 



L I 



Y^{2l + l)Pi{cosB) 




gifcr ^i{kr—ln) 



[l + 2ikfi{k)] 



2ikr 



+ ^{2l + l)fi{k)Pi{cose) 



,ikr pi{kr—lir) 



r 



(16) 



La ecuacion anterior la podemos interpretar de la manera siguiente. Los 
dos terminos exponenciales corresponden a ondas esfericas, el primero a 
una onda emergente y el segundo a una onda convergente; y el efecto de 
la dispersion se ve en el coe&ciente de la onda emergente, y es igual a uno 
cuando no existe un dispersor. 

Corrimientos de fase 

Imaginemos ahora una superficie cerrada centrada en el dispersor, si asum- 
imos que no hay creacion ni aniquilacion de particulas se verifica: 



Donde la region de integracion es evidentemente la superficie anterior- 
mente deRnida y j es ia densidad de corriente de probabilidad. Aun mas, 
debido a la conservacion del momento angular la ecuacion anterior debe ver- 
i&carse para cada onda parcial (en otras palabras, todas las ondas parciales 
tienen diferentes valores de las proyecciones del momento angular, lo cual las 
hace en esencia diferentes y la formulacion ser/a equivalente si trataramos al 
paquete de ondas como un haz de particulas que no interactuan entre si; mas 
aun, debido a que el potencial de nuestro problema es central el momento 




(17) 
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angular de cada "part/cuia" se conservara de tal manera que podnamos de- 
cir que las particulas continuan siendo las mismas). Pot las consideraciones 
anteriores, podemos decir que tanto la onda divergente como la emergente 
difieren a lo mucho en un factor de fase, es decir, si de&nimos: 

Si{k) = l + 2ikfi{k) (18) 

debera suceder que 

I Si{k) \= 1 (19) 

Los resultados anteriores los podemos interpretar en vista de la conser- 
vacion de las probabilidades, y era de esperarse pues hemos supuesto que no 
existe creacion ni aniquilacion de particulas, asi que la influencia del disper- 
sor consiste en agregar simplemente un factor de fase en las componentes de 
la onda emergente, en virtud de la unitariedad del factor de fase lo podemos 
escribir como: 

Si = e^'^i (20) 

Donde 6i es real y es funcion de k. A partir de la definicion (18) podemos 
escribir: 

2ik k kcot{5i)-ik ^ ' 

La seccion total de dispersion adquiere la siguiente forma: 

ototai = j I S{e)?dn = 

^f"d<j) d(cos(e))^^(2/ + l){2l' + l)e'^' sin((5/)e*V sin((5;0-P/^i' 
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A'jr 

= ^E(2^ + l)«i^'('^r) (22) 

Determinacion de los corrimientos de fase. 

Consideremos ahora un potencial V tal que V se anula para r > R, el 
parametro R se le conoce como el "alcance del potencial", asi que en la 
region r > R evidentemente debe corresponder a una onda esferica libre. 
For otro lado la forma mas general de la expansion de una onda plana en 
ondas esfericas es de la forma: 

(x|^+) = ^^^E^'(2^ + l)AW^'Kcos^) ir>R) (23) 

Donde el coeficiente Ai esta definido como: 

Ai = c\^^ hf^ {kr) + cp^ hf (kr) (24) 

donde h\^^ y hf^ son las funciones de Hankel esfericas y sus formas 
asintoticas estan dadas por: 

i{kr-l7r/2} 

^ ^ ~ ^-I — 

ikr 

-i{kr~l7T/2) 

~ - —u 

ikr 

AI ver la forma asintotica de la expresion (23): 
1 



De esta manera vemos que: 



2ikr 2ikr 



(25) 



cp) = \e-^^ cp) = 1 (26) 
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Ahora podemos ver que la funcion de onda radial para r > R se escribe 
como: 

Ai = e^'^i [cos 6iji{kr) - sin 6ini{kr)] (27) 

A partir de ecuacion anterior podemos evaluar su derivada logaritmica 
en r=R, i.e., justo afuera del alcance del potencial: 

j'l cos 6i — n'l (kR) sin 5i 



(28) 



ji cos Si -ni{kR) sin Si 
Donde Ji es la derivada de ji con respecto akr y evaluada enr = R. Otro 
resultado importante que podemos obtener conociendo el resultado anterior 
es el corrimiento de fase: 

kRj'AkR) - pdiikR) 

tan Si = ^7 — { — — T (29) 

kRn[{kR) - pmiikR) ^ ^ 

Para encontrar la solucion completa a nuestro problema aun hay que 

hacer los calculos para cuando r < R, es decir, dentro del rango del alcance 

del potencial. Para el caso de un potencial central la solucion a la ecuacion 
de Scbrddinger en tres dimensiones se reduce a resolver la ecuacion: 

(fui /,2 2m^^ + , X 

-^^[^'~^y-^y^=^ (30) 

donde ui = rAi{r) y esta sujeto a la condicion de frontera ui \r=o = 
de esta manera ya podemos calcular la derivada logaritmica, la cual por la 
propiedad de continuidad de la derivada logaritmica (que es equivalente a 
la continuidad de la derivada en un punto de discontinuidad): 

A \dentro= A 1/ uera 

(31) 
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Un ejemplo: dispersion por una esfera solida. 

Ahora tratemos un caso especi&co. Consideremos un potencial definido 
como: 

^=1^ (32) 

] r > R ^ ^ 

Es sabido que una particula no puede penetrar en una region donde el 
potencial sea infinito, as/ que la funcion de onda se debera anular enr = R; 

de que la esfera es impenetrable tambien se deduce que: 

Mr) \r=R= (33) 
Asi de la ecuacion (27) tenemos: 

ni{kR) 

Vemos que se puede calcular facilmente el corrimiento de fase para 
cualquier 1. Consideremos ahora el caso I = (dispersion de una onda 
S) de esta forma tenemos: 

6i = -kR 

y de la ecuacion (27): 

. , , sin/cr . cos kr . . 1 . /, r- x 
Ai=o{r) ~ — - — cos do H r — sin do = — sin(A;r + So) (35) 

KV KV KT 

Vemos que la contribucion al movimiento libre es la adicion de una fase. 
Claro que en el caso mas general diferentes ondas tendran diferentes corrim- 
ientos de fase ocasionando una distorsion transitoria en el paquete de ondas 
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dispersado. Estudiemos ahora el caso de energias pequenas, i.e., kR « 1 
en este caso las expresiones para las funciones de Bessel (usadas para escribir 
las funciones de Hankel esfericas) son las siguientes: 



obteniendo as/ la expresion: 



= (a + i)[(a-i)!!p ■ (^«' 

De la formula anterior podemos ver que la contribucion apreciable al 
corrimiento de fase es principalmente de ondas con I = 0, pero como Sq = 
—kR tenemos para la seccion eficaz de dispersion: 



total 



^dn = AttR^ (39) 



Vemos que la seccion eficaz cuantica es cuatro veces mayor que la seccion 
eficaz clasica, y coincide con cl area total de la esfera. Para grandes valores 
de la energia del paquete incidente conjeturemos que todos los valores de I 
hasta un valor maximo Imax ~ kR contribuyen a la seccion eficaz total. 

(Ttotal = Tf E + • (40) 

1=0 

De esta forma a partir de la ecuacion (34) tenemos: 



sm^Si = ^ = , — „ ~ sin^ { kR . (41) 

l + tan2,5, \ji{kR)? + [ni{kR)]^ \ 2 > ^ ^ 
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Aqui hemos utilizado las expresiones: 

ni(kr) ~ — — cos (kr 

vemos que 6i decrece en | cada vez que I se incrementa en una unidad, 
asi es evidente que se cumple sin^Si + sin^(5;_|_i = 1, y aproximando sin^5/ 
por su valor promedio \, asi es sumamente sencillo obtener el resultado de 
la suma anterior, y la formula para la seccion e&caz total es: 

atotal = ^{kRfl = 2TrR^ (42) 

Una vez mas el resultado del calculo utilizando mecanica cuantica diBere 

del resultado clasico, pero veamos cual es el origen del factor 2; primero a 
la ecuacion (15) la vamos a separar en dos partes: 



l=kR ■ l=kR 

/(^) = ^ E {2l+l)e'^'' Pi cosid)+^ (2^+l)P/cos(^) = freflexidn+fsombra 

(43) 



evaluando J \f reflexionl^ '^^ 

27r ^ "'"f f ^ 

\frefiexi6nfd^ = ^Il {2l+lf[Picom?d{cose) 



^max p 1 7 2 

" " '^hnax 7->2 

, , —[J- = 

1=0- 



fe2 

(44) 



Ahora analizando f gQ^^iyj-Q para angulos pequeHos tenemos: 
f sombra ~ ^ E(2^ + l)'^o(i^) ~ ik j\Mkhe)dh = '3^1^ (45) 
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Esta formula es hastante conocida en optica, es la formula de difraccion 
de Fraunhofer; haciendo el cambio de variable z = kRO podemos evaluar la 
integral J IfsomhrXdCl: 

1 1/ sombrafd^ ~ ^^^Vq'^^^T^'^^ ~ ""^^ ^^^^ 

Finalmente despreciando la interferencia entre f reflexion ^ f sombra 
(porque la fase oscila entre 2Si+i = 2Si — tt). De esta manera obtenemos 
el resultado (42). Hemos etiquetado un termino con el nombre de sombra, 
el origen de esta contribucion se explica facilmente al apelar al compor- 
tamiento ondulatorio de la particula dispersada (en este punto da lo mismo 
[fisicamente] un paquete de ondas que una particula), tiene su origen en 
las componentes del paquete dispersadas hacia atras, entonces tendran una 
diferencia de fase con las ondas incidentes y habra una interferencia destruc- 
tiva. 

Dispersion en un campo de Coulomb 

Ahora consideremos un ejemplo clasico y algo complicado: la dispersion 
de particulas en un campo coulombiano, para este caso la ecuacion de 
Schrddinger es: 

(-|^V^-^)^W=^^(r), ^>0 (47) 

donde m es la masa reducida de las particulas que interaccionan y eviden- 
temente E > debido a que tratamos el caso de dispersion sin creacion de 
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estados ligados. La ecuacion anterior es equivalente a la siguiente expresion 
(para valores adecuados de las constantes k y j : 

V2 + A;2 + ^^^(r) = o (48) 

Si no consideramos la parte centrifuga del potencial efectivo, es decir, 
que nos quedamos unicamente con la interaccion de un campo de Coulomb 
puro, por eso podemos proponer una solucion de la forma: 

i^ir) = ^^^-^xiu) (49) 

con 

u = ikr{l — cos 9) = ik{r — z) = ikw 
\<i-v = kz 

'tp{r) es la solucion completa a la ecuacion de Schrddinger, y es de esper- 
arse que su comportamiento asintotico conste de dos partes: una tipo onda 
plana e**^'^; y otra tipo onda esferica r~^e'^^^ . Definiendo nuevas variables: 

z = z w = r — z X = (p 

y con ayuda de las relaciones anter lores, la ecuacion (48) adquiere la 
forma: 

r w2 w 1 

X{u) = (50) 

Para resolver la ecuacion anterior, primero hay que estudiar sus compor- 
tamientos asintoticos, pero como todo eso ya esta hecho, la funcion de onda 
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asintotica normalizada que se ohtiene al final de todos los calculos anteriores 
es: 

Mr) = (^e^[k-r-7^n(..-k.r)] + fc{k, e)e^^^+^^-'''^ ^ ^^^^ 

Como vemos, la funcion de onda anterior tiene algunos terminos que la 
hacen diferir de manera apreciable a nuestra ecuacion (7) y eso se debe a que 
una fuerza de tipo de Coulomb es de largo alcance. Hacer el calculo exacto 
para la amplitud de dispersion de Coulomb es bastante dificil de hacer (de 
hecho casi todos los calculos de este problema). Quien logre hacer todos los 
calculos para obtener la funcion de onda normalizada encontrara: 

1 / n*('M^ .-,i'kr+-l ii2kr]\ 

dondegii-f) = ■ 

El analisis de ondas paiciales lo reduciremos a mostrar resultados ya 
obtenidos de una manera lo mas clara posible. Primero escribamos la funcion 
de onda (49) V'(r) de la siguiente manera: 

^(r) = e'^-^xiu) = Ae^^ '' j e"*f^-^(l - t)-'Ut (53) 

J c 

donde A es una constante de normalizacion y toda la parte integral es la 
transformada de Laplace inversa de la transformada directa de la ecuacion 
(50). La ecuacion anterior se puede escribir de la siguiente forma: 

V'(r) = A f e*"'(l - t)e''''^\l - t)d{t, j)dt (54) 
J c 
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con 

d(t,7) =fT-i(l-i)-^T-i (55) 
Ya en el analisis de ondas parciales procedemos a escribir: 

OD 

^(r) = ^(2i + l)i^Pi{cose)Ai{kr) (56) 

1=0 

donde 

Mkr) = a[ e'^^^'jiikril - t)](l - t)d{t, 7) (57) 
Jc 

Y dada la relacion que hay entre las funciones de Bessel esfericas con las 
funciones esfericas de Hankel tenemos: 

Ai (kr) = a\^^ (kr) + ^ (kr) (58) 

La evaluacion de los coeBcientes anteriores no la vamos a hacer aqui 
(dada su extension ), nos basta con saber que: 

A'i^\kr) = (59) 

(60) 



donde 



"r(i + / + i7) ^^'^ 
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Calculo de la amplitud de dispersion de Coulomb 

Si tomamos la transformada de Laplace de la ecuacion (50) obtenemos: 

X{u) = A I e"*fT-^(l - t)-^^dt (62) 
Jc 

Donde el contorno C va de — oo a oo j se cierra por arriba del eje real, 
vemos que hay dos polos: cuando t = y t = 1. Haciendo el cambio de 
variables s = ut obtenemos: 

X{u) = A I e's^'^-^iu - s)-^^ (63) 

x('u) debe ser regular en cero y en efecto: 

X{0) = {-ly'-^A I -ds = {-ly'^Al-Ki (64) 

Ahora tomando el limite — ^ oo haciendo un desplazamiento infinitesi- 
mal (para quitarse el problema de que los polos estan sobre la trayectoria) y 
un cambio de variable tal que ^ = — pero como este termino tiende 

a cero cuando n — ^ oo entonces podemos expander (u — s) en series de po- 
tencias de ^ para el polo con s = 0. Pero el desarroUo anterior no es bueno 
en s = 1, porque aqm s = —sq + z(k ± e; y de ahi que ;^ = 1 — tiende 
a 1 cuando k — ^ oo; pero si hacemos el cambio de variable s = s — u ya no 
hay problema: 

X{u) = ^j^ {[e'^^-'^{u - s)-^^]ds + [e^'+"(-s')^^(^i + s'y-^-'^]ds^ (65) 
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DesarroUando las series de potencias es sencillo calcular las integrales 
anteriores, pero del resultado aun hay que tomar el limite f — para 
obtener las formas asintoticas correctas para la dispersion de Coulomb: 

27r5i(7) = i [ e's'^-^ds 
J C2 

27r52(7) = i [ e's-'^ds (66) 

JC2 

Y despues de tanto cambio de variables, regresamos a las originales y 
vemos: 

{u*y^ = {-iy^[k{r - z)Y^ = eT"^/2e^Ti°'=('-^) 

= {{)-''^[k{r - z)]-''^ = eW2g-i7lnfc(r-^) (67) 
Ya hemos calculado x, Jo cual equivale a tener V'k(r) a partir de (49). 

Aproximacion eikonal 

Hacemos una breve exposicion a cerca de la aproximacion eikonal. La 
aiosofia de la aproximacion eikonal en este caso es el mismo que cuando 
hacemos el paso de la optica ondulatoria a la optica geometrica, por eso es 
valida cuando cl potcncial varia mas lentamente que la longitud de onda 
del paquete de ondas dispersado, es decir, para el caso E » \V\. Asi esta 
aproximacion puede ser considerada como una aproximacion cuasiclasica. 
Primero propongamos que la funcion de onda cuasiclasica es: 
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Donde S satisface la ecuacion de Hamilton-Jacobi, y su solucion es: 
S 



1/2 

dz' + constante (69) 



La constante aditiva se escoje de tal forma que: 
g 

-^^kz cuando V ^0 (70) 

El termino que multiplica el potencial lo podemos interpretar como un 
cambio de fase del paquete de ondas, su forma espedfica es la siguiente 

La anterior aproximacion dentro del metodo de ondas parciales tiene la 
siguiente aplicacion. Sabemos que la aproximacion eikonal es valida para 
grandes energias, pero para alias energias hay muchas ondas parciales que 
contribuyen a la dispersion, asi que podemos tratar a I como una variable 
continuay por analogia con la mecanica clasica: I = bk. Ademas como ante- 
riormente ya habiamos mencionado Imax = kR, sustituyendo en la expresion 
(15) se obtiene: 

f{0) = -ik J hJo{kbe)[e^'^^^'^ - l]db (72) 
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